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Abstract 



We discuss exclusive charmless B decays within the Standard Model of particle 
physics. These decays play a central role in the on-going process to constrain the 
parameters of the CKM matrix and to clarify the nature of CP violation. In order to 
exploit the rich source of data that is currently being collected at the experiments, 
a systematic theoretical treatment of the complicated hadronic dynamics is strongly 
desired. QCD Factorization represents a model-independent framework to compute 
hadronic matrix elements from first principles. It is based on a power expansion in 
and allows for the systematic implementation of perturbative corrections. 

In particular, we consider hadronic two-body decays as -B — tttt and perform a con- 
ceptual analysis of heavy-to-light form factors which encode the strong interaction 
effects in semi-leptonic decays as B Triu. 

Concerning the hadronic decays we compute NNLO QCD corrections which are 
particularly important with respect to strong interaction phases and hence direct 
CP asymmetries. On the technical level, we perform a 2-loop calculation which is 
based on an automatized reduction algorithm and apply sophisticated techniques 
for the calculation of loop-integrals. We indeed find that the considered quantities 
are well-defined as predicted by QCD Factorization, which is the result of a highly 
complicated subtraction procedure. We present results for the imaginary part of 
the topological tree amplitudes and observe that the considered corrections are sub- 
stantial. The calculation of the real part of the amplitudes is far more complicated 
and we present a preliminary result which is based on certain simplifications. Our 
calculation is one part of the full NNLO analysis of nonleptonic B decays within 
QCD Factorization which is currently pursued by various groups. 
In our conceptual analysis of the QCD dynamics in heavy-to-light transitions we 
consider form factors between non-relativistic bound states which can be addressed 
in perturbation theory. We perform a NLO analysis of these form factors and discuss 
some open questions of the general factorization formula which is obtained from the 
heavy-quark expansion in QCD. These include the origin and resummation of large 
logarithms and the non-factorization of soft and collinear effects in the so-called 
soft-overlap contribution. We show that the latter can be calculated in our set-up 
and address the issue of endpoint singularities. As a byproduct of our analysis, we 
calculate leading-twist light-cone distribution amplitudes for non-relativistic bound 
states which can be applied for the description of Be and rjc mesons. 
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Introduction 



The analyses presented in this thesis rely on the Standard Model of particle physics 
P-[9j which reflects our current knowledge of three of the four known fundamental 
forces in nature. Electromagnetic, weak and strong interactions are described therein 
by a relativistic and renormalizable quantum field theory which is based on a gauge 
principle. The Standard Model represents an impressive theoretical achievement 
which successfully explains phenomena from everyday electricity to high-energetic 
quantum processes that are investigated at dedicated particle accelerator facilities. 
It is one of the best-tested theories of contemporary physics. This may be illustrated 
by comparing the experimentally measured value of the anomalous magnetic moment 
of the muoru 

a^^^P = 11 659 208.0 (5.4)(3.3) x 10"^° (1) 

with its theoretical prediction calculated within the Standard Model [lU] 

af^ = 11 659 185.8 (7.2)(3.5)(0.3) x 10"^°. (2) 

Despite its tremendous success the Standard Model has its insufficiencies and is 
commonly believed to be incomplete. Severe constraints from electroweak precision 
data may be interpreted as a hint that "something unknown" happens at the TeV- 
scale, an energy scale that has been out of the scope of todays collider experiments. 
The Large Hadron Collider (LHC), which is currently being built at CERN and is 
scheduled to start operation in mid-2007, has particularly been designed to explore 
the physics at the TeV-scale. Thousands of particle physicists from all around the 
world are looking forward to the first data taking of the LHC, in the hope that it 
will help us to reveal the limitations of the Standard Model and give a first clue 
about the theory that lies beyond it. 

Apart from these direct searches, the physical effects that are supposed to lie beyond 
the Standard Model (often referred to as New Physics) can be investigated indirectly 
in high-precision measurements of low-energy observables as the anomalous magnetic 
moment of the muon mentioned above. Every physical observable is in principle 
sensitive to arbitrarily high energies and thus to New Physics due to quantum effects. 

^The anomalous magnetic moment of the electron can be measured even more precisely. As it 
is less likely to be affected by physics beyond the Standard Model, its measurement is used for the 
determination of one of the Standard Model parameters, namely the fine structure constant [llj . 
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The quantitative investigation of these tiny effects represents a highly challenging 
task, both for experimental measurements and for theoretical calculations which 
have to match the experimental accuracy. 

Whereas the gauge sector of the Standard Model has been tested to remarkable 
precision in the era of the Large Electron Positron Collider (LEP), the flavour sector 
is experimentally less constrained. The flavour sector contains a large number of 
parameters as the quark and the lepton masses or the four parameters related to 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [T2IIT3] which describes the mixing 
of the quark mass eigenstates in weak interactions. The numerical values of these 
parameters are not predicted by the Standard Model but rather have to be extracted 
from experimental measurements before making any theoretical prediction. 

The ultimate goal of B physics is to precisely determine some of these parameters 
and to test the CKM sector of the Standard Model. The phenomenon of CP violation 
is of particular interest as it is related to the observed matter-antimatter asymmetry 
in our universe. Two CKM parameters are known to date at the percent level [13] 

A = 0.22717lHS}g?, A = 0.806l°:°J^ (3) 

and the remaining two parameters p and f] are conveniently discussed in the context 
of a unitarity triangle which reflects the unitarity of the CKM matrix. The current 
status of the unitarity triangle is shown in Figure [TJ The values of the parameters 
p and f] correspond to the upper tip of the triangle which is given by [T^ 

P = 0.195l°:°i, = 0.3261°:°?!. (4) 
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Figure 1: Global fit of the unitarity triangle (status September 2006) I14f - 
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Notice that these determinations are far less accurate than the ones of A and A 
and that various independent measurements of p and f/, which are indicated by the 
coloured bands in Figure [H, are all consistent with each other. The very fact that 
the area of the unitarity triangle is non-zero is a manifestation of CP violation. 

Two dedicated B factories, the BaBar experiment at SLAG and the Belle experiment 
at KeK, have contributed substantially to our current understanding of B physics 
and CP violation. One of the most important milestones in their physics programme 
was the observation of (indirect) CP violation in the neutral B system in 2001 [15.16j 
more than 35 years after the first observation of CP violation by James Cronin, Val 
Fitch and collaborators in the Kaon system jT7] . 

Since the start of BaBar and Belle in 1999, the B factories have produced ~ 10^ pairs 
of B mesons which correspond to ~ 1 ab~^ of integrated luminosity. This impressive 
wealth of data will continuously increase until the end of their physics programme 
in 2010. It will then be complemented by the LHC-b experiment and future plans 
concerning an upgrade of the B factories are already envisaged. With this ongoing 
experimental effort it will be possible to nicely overconstrain the unitarity triangle 
and to reduce the uncertainty of the parameters p and f/ to a few percent. 

The reader may wonder why some observables as the anomalous magnetic moment 
of the muon can be determined with an uncertainty of better than 1 part in 10^ and 
why it is apparently so difficult to reduce the uncertainty of the CKM parameters 
at the percent level. One of the reasons is that the CKM parameters are related to 
quarks which are always affected by strong interaction^. 

In the Standard Model the strong interactions are described by a non-abelian gauge 
theory called Quantum Chromo Dynamics (QCD) The most important prop- 

erty of QCD is asymptotic freedom, i.e. the fact that the coupling of the quarks to 
the gluons becomes weak at large energy scales (cf. Figure [2]). Strong interaction 
effects from high energies as e.g. the Z° resonance at Mz = 91.2 GeV can be pre- 
cisely calculated in perturbation theory as an expansion in the coupling constant 
as{Mz) ~ 0.12. At lower energies of order of the QCD scale Aqcd ~ 0.5 GeV, the 
perturbative expansion in the coupling constant breaks down and the quarks get 
confined into complicated colour-singlet hadrons. In our theoretical description of 
any hadronic process that is observed in experiment we thus have to deal with these 
non-perturbative effects. We may sometimes be lucky and find some observables 
which are almost free of these hadronic uncertainties as the "golden" sin 2/3 mea- 
surement from the B ^ J/ip Ks decay which corresponds to the thin dark blue ray 
in Figure [U As it is not always possible to find such clean observables, we have to 
look for sophisticated methods which allow to control the hadronic dynamics. 

In this thesis we mainly deal with strong interaction effects in B meson decays. The 
role of B mesons is special as they are (apart from the T-resonances) the heaviest 
mesons which show an extremely rich phenomenology including many interesting CP 



^Even though the anomalous magnetic moment of the muon is a leptonic quantity, the main 
hmitation in its determination also stems from small strong interaction effects. 
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Figure 2: Running of strong coupling constant llSf . 

violating observables. Furthermore, the hadronic dynamics in B meson decays turns 
out to be partly accessible within perturbation theory as the intrinsic mass of the 
6-quark mf, ~ 5 GeV corresponds to a perturbative scale with as{mb) ~ 0.22. More 
precisely, the weak decay of the fe-quark is accompanied by short-distance effects 
and unaffected by the hadronization of the quarks which occurs at much larger 
time scales. The technical procedure which disentangles perturbative from non- 
perturbative effects is called factorization. For more details concerning factorization 
we refer to the comprehensive introduction in the first chapter of this thesis. 

We in particular consider exclusive charmless decays of B mesons, as B ^ mi or 
B Tciu. These decays provide important information on the unitarity triangle, 
the former serve for the determination of the CKM angle a and the latter for the 
measurement of \Vub\- The impact of these decays on the unitarity triangle from 
Figure [H are reflected by the light blue band and the dark green circle, respectively. 
Rather than focussing on the phenomenological aspects of these decays, we examine 
the complicated strong interaction dynamics which is encoded in the hadronic matrix 
elements. 

The decays considered in this thesis share the feature that the B meson decays 
into very energetic light mesons (in the B meson rest frame). The factorization of 
short- and long-distance QCD effects for these heavy-to-light transitions has first 
been worked out in the framework of QCD Factorization [TUH^ and has later been 
translated into a field theoretical formulation which is called Soft-CoUinear Effective 
Theory (SCET) [22H21] . We are mainly concerned with the calculation of higher 
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order perturbative QCD corrections, but have also a general look at the factorization 
properties of heavy-to-light form factors which are among the "simplest" objects for 
studying the QCD dynamics in B ^ 71 transitions. 

The core of this thesis consists in a 2-loop calculation related to the hadronic two- 
body decay modes. Due to the complexity of the problem we split the calculation 
into two parts. We first compute the imaginary part of the hadronic matrix ele- 
ments which is technically simpler than the real part. It is in addition of particular 
interest in phenomenological applications as it is related to a strong phase shift 
between the final state mesons which "pollutes" the interesting information about 
the underlying weak CKM phases. Our calculation represents one part of the full 
next-to-next-to-leading order (NNLO) analysis of the topological tree amplitudes 
within QCD Factorization. Another part of this analysis concerning (1-loop) spec- 
tator scattering has been calculated recently by various groups [251427] . We remark 
that the phenomenological impact of our corrections is beyond the scope of this 
thesis as the full NNLO calculation (including topological penguin amphtudes) is 
still incomplete. 

Our conceptual analysis of heavy-to-light form factors is based on a particular sce- 
nario. We consider transition form factors between non-relativistic bound states 
which can be addressed in perturbation theory. We perform a next-to-leading order 
(NLO) analysis of these form factors using the same techniques that we have devel- 
oped for the aforementioned 2-loop calculation. We then compare our explicit results 
with the general factorization formula for heavy-to-light form factors and address 
some of its open questions concerning the so-called soft-overlap contribution, the 
issue of endpoint singularities and the resummation of large (Sudakov) logarithms. 
We emphasize that this analysis can be applied for the description of Be — > (cc) 
transitions, although we focus on the more conceptual aspects here. 
Large parts of this thesis deal with the calculation of loop-integrals. One important 
element in our calculation procedure is an automatized reduction algorithm which 
allows us to reduce the calculation of several thousands of loop-integrals to a much 
smaller set of so-called Master Integrals. The calculation of these Master Integrals 
represents the most difficult task of the entire calculation. We would like to point 
out that the technical difficulties that we encountered in the two considered calcu- 
lations are of rather different origins. Whereas the first calculation in the B —>■ irn 
context represents a highly challenging and complex 2-loop calculation, the other 
(1-loop) calculation related to the form factors is complicated due to the presence 
of various physical scales. In this case we restrict our attention to the leading power 
in a mass expansion which requires different techniques than for full loop-integrals. 
Consequently, our collection of loop-techniques that we present in the second chap- 
ter of this thesis, summarizes (almost) all of the most sophisticated techniques that 
have been developed so far: the method of differential equations the for- 

malism of harmonic polylogarithms ^U], the method of expansion by regions [HT] , 
Mellin-Barnes techniques [32)133] and the method of sector decomposition [3l] . Most 
of these techniques have rarely been applied in B physics so far. 
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The structure of this thesis can be outlined as follows: 

In the first chapter we present the theoretical background required for an analysis 
of exclusive charmless B decays. We give a comprehensive introduction to QCD 
Factorization and Soft-Collinear Effective Theory and comment briefly on several 
alternative approaches. The remainder of this chapter is devoted to a detailed 
analysis of hadronic two-body decays and heavy-to-light form factors which are of 
particular interest in this work. 

In Chapter [2] we collect the techniques that we have used in our calculations. We 
develop a systematic strategy which is based on an automatized reduction algorithm. 
As most parts of this thesis deal with the calculation of loop-integrals, we dedicate 
a sizeable part of Chapter 2 to the presentation of several sophisticated techniques. 

Due to this structure, the second part of this thesis is free of the technical issues 
related to the calculation of loop-diagrams. In Chapter [3] we consider the imaginary 
part of the topological tree amplitudes in hadronic two-body decays. Apart from the 
2-loop calculation, we address the issues of Fierz symmetry, evanescent operators, 
renormalization and IR subtractions. We finally obtain our results in an analytical 
form and conclude this chapter with a brief numerical analysis. 
Chapter H] is devoted to the real part of the topological tree amplitudes in hadronic 
two-body decays. The calculation follows the same lines as in Chapter O but is 
technically much more involved. So far, we have accomplished the technical part of 
this calculation and present some preliminary numerical results. 
In our final analysis in Chapter [5] we consider heavy-to-light form factors between 
non-relativistic bound states. We first address the non-relativistic approximation 
in this context and present subsequently the NLO (1-loop) calculation of the form 
factors. In our conceptual analysis we investigate the origin of endpoint singularities 
and comment on the resummation of logarithms. We show that we can isolate (and 
calculate) the so-called soft-overlap contribution in our set-up and calculate leading 
twist light-cone distribution amplitudes of non-relativistic bound states. 

We finally conclude and give an outlook on future developments. The results for all 
Master Integrals that appeared in our calculations from Chapter [3] - [5] as well as the 
explicit expressions for the NLO form factors are summarized in the appendix. 
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Chapter 1 

Exclusive charmless B decays 



In the first chapter we present the theoretical background for the description of 
exclusive charmless B decays. As the perturbative calculations in the second part of 
this thesis are based on QCD Factorization and Soft-CoUinear Effective Theory, we 
give a profound introduction to these two developments. We paid special attention 
to avoid unnecessary formulas in the introductory chapter in order to allow for a 
transparent presentation of the basic ideas behind these concepts. We comment 
briefly on alternative approaches to charmless B decays and have a closer look at 
those decays which are of particular interest in the work at hand. These include 
hadronic two-body decays and a conceptual analysis of heavy-to-light form factors 
which are important ingredients in semi-leptonic and radiative B decays. 

1.1 Basic concepts 
1.1.1 QCD Factorization 

The phrase QCD Factorization is closely related to the theory of hadronic two-body 
decays of B mesons. At the end of the 90s Beneke, Buchalla, Neubert and Sachrajda, 
to which we refer as BENS in the following, established this novel framework which 
allowed for the first time for a systematic treatment of these decays in QCD [T9ti2T] . 
However, QCD Factorization is a more general framework with applications covering 
a wide spectrum of semi-leptonic, radiative and hadronic B decays. 

QCD Factorization basically merged two different developments that were known 
at that time: the heavy quark expansion (for a review see [35]) and the theory of 
hard exclusive processes which is also known as coUinear factorization [36)137 ]. In 
the course of the 90s much progress has been made in the understanding of heavy 
mesons. It has been realized that the QCD dynamics of heavy mesons simplifies 
substantially when it is considered in the heavy quark expansioE0 (HQE), i.e. an 
expansion in the ratio Aqcd / "^q where tuq is the mass of the heavy quark. On the 

"'^The phrase heavy quark expansion is often used in the Hterature in the context of inclusive 
decays. Wc will refer to it here more generally whenever we speak about an expansion in Aqcd /niQ. 
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other hand the theory of hard exclusive processes can be seen as the counterpart of 
deep-inelastic scattering for inclusive processes. It was developed for the description 
of exclusive processes with a large momentum transfer ^ ^qcd- -Due to the large 
momentum transfer, the particles in these processes are very energetic and assumed 
to move coUinear to light-cone directions which leads to important simplifications. 
Charmless B decays naturally incorporate both of these aspects. The B meson in 
the initial state implies a systematic description in terms of a HQE in AqcB/^b- 
The final state being charmless, which means that it consists of light hadrons only, 
further implies that the particles in the final state are very energetic in the B meson 
rest frame and can be described to move almost on the light-cone. 

The basic idea of factorization is the attempt to disentangle physical effects from 
different length or momentum scales. This is a very general idea that can be ap- 
plied in many different fields of physics. Concerning the dynamics of the strong 
interactions this strategy is particularly suited due to the asymptotic freedom of 
QCD. Any decay or scattering process involving hadrons is sensitive to the scale 
Aqcd which is responsible for the confinement of the quarks into the hadrons. As 
the strong coupling constant at these scales is of 0{1), the respective effects cannot 
be addressed in perturbation theory and therefore we call them non-perturbative. 
In B physics we are confronted with an additional intrinsic scale in form of the 
mass of the b-quark rrib ^ Aqcd which is a perturbative scale with as{mb) ~ 0.22. 
We see that the idea of disentanghng the effects from the scales nib and Aqcd is 
equivalent to separating perturbative from non-perturbative effects in QCD. The 
predictive power of factorization lies in the fact that we can calculate the former 
systematically in perturbation theory whereas the latter typically give rise to univer- 
sal hadronic quantities which can be obtained from other methods as lattice gauge 
theory or QCD sum rules or they can even be extracted from experimental data. 
The essence of the QCD Factorization prediction is summarized in a factorization 
formula for a hadronic matrix element. The factorization formula illustrates how 
the perturbative and non-perturbative effects are disentangled {factorized) . In the 
remainder of this section we present several examples of factorization formulas for 
different classes of exclusive processes. In the first two examples we sketch the 
situation that was known before QCD Factorization was established. The first one 
deals with a heavy-to-heavy transition which can be described with the help of the 
HQE and the second one corresponds to coUinear factorization. In the last two 
examples we illustrate how QCD Factorization combines these two pictures. We 
give slightly simplified descriptions in order to concentrate on the main aspects 
concerning factorization. We hope that our presentation will help to understand the 
structure of the factorization formulas that we discuss in the following sections. 

Example 1: B ^ Dlv 

We start with the simplest example in form of exclusive B decays into final states 
that do not contain light hadrons as e.g. B Diu. The relevant scales in these 
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Figure 1.1: Factorization of short- and 
long distance effects in heavy-to-heavy 
transitions. The former are contained in a 
coefficient function H , the latter in a soft- 
overlap contribution S. The double lines 
denote heavy quarks with mq ^ Aqqd- 

processes are the mass(es) of the heavy quark(s) which we simply denote by mg and 
the hadronic scale Aqcd- The factorization formula for a generic hadronic matrix 
element of a current J' takes the schematic form 

{D\J\B) = Hifip) Sifip) + 0(AQCD/mQ). (1.1) 

The factorization formula is illustrated in Figure [TTTl We first notice that the factor- 
ization formula makes a statement about the leading power in the HQE. On the other 
hand it is predicted to be valid to all orders in perturbation theory. In writing (11.11) . 
hard effects from the scale mq and soft effects related to Aqcd have been disentan- 
gled, the former being contained in the coefficient function H{fip) and the latter in 
a remnant matrix element which we denoted by Slfip). Technically, factorization is 
achieved with the help of a factorization scale fip satisfying mq ^ fip ^ Aqcd- The 
effects form hard gluons with virtualities k"^ > n'j^ are contained in H{^f) and those 
from soft gluons with k"^ < are absorbed into S^fip)- As the factorization scale 
has been introduced artificially in the factorization formula (11.11) . the dependence 
of the functions H{jip) and S{hf) has to cancel in their product. 

As a side remark we mention that the concept presented above applies as well to 
the effective weak interactions in the description of low-energetic hadronic processes. 
These processes provide the hierarchy ^ where My/ is the mass of the W- 
boson and is a typical momentum scale in the process. In leading power in an 
expansion in q^ /M^r the hadronic matrix elements factorize similar to (II. ip into 
short- distance Wilson coefficients, which correspond to the H{hf) in our notation, 
and remnant matrix elements S{^f) which can be calculated in the Fermi theory 
of weak interactions. We come back to the effective weak interactions when we 
consider hadronic two-body decays in Section 11.21 

Example 2: nj* — > vr 

In a second example we consider hard exclusive processes which are not related 
to B physics, as e.g. tc'j* —>■ tt where the kinematics sets a perturbative scale in 
form of the large momentum transfer ^ ^qcd • Because of the large momen- 
tum transfer, light-cone dynamics comes into play and the energetic light mesons 
can approximately be described by their two-particle quark-antiquark Fock states. 
Similar to what is done in deep-inelastic scattering, the quark and the antiquark 
can be assumed to move coUinearly inside the meson and share its momentum with 
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Figure 1.2: CoUinear factorization of 
short- and long distance effects. The for- 
mer give rise to a hard- scattering ker- 
nel T, the latter to light-cone distribution 
amplitudes (j)-^. The dashed lines denote 
collinear quarks. 



fractions u and u = 1 — u, respectively. The corresponding factorization formula is 
sketched in Figure 11.21 and reads 



{n\J\n) = du dv (j)^{v; hf) T{v,u; hf) (pniu; fip) + 0{Aqcj)/Q^). (1.2) 



Again, perturbative and non-perturbative effects are systematically disentangled in 
leading power. In this case the hard effects from the scale give rise to a hard- 
scattering kernel T{v, u; ^f) which depends on the momentum fractions v and u of 
the quarks in the mesons. The collinear effects from the scale Aqqj-, are encoded 
in light-cone distribution amplitudes (p-niv; fip) and 4>n{u; iif) of the initial and final 
state particles. 

Let us stress two important differences between the factorization formulas fll.ip and 
fll.2p . First, in (11.21) there is no long-distance interaction between the initial and the 
final state particle at leading power. The interaction between both mesons is entirely 
described by the perturbative hard-scattering kernel. The non-perturbative input in 
form of the distribution amplitudes contains information about the structure of the 
participating mesons and is independent of the considered process. This is different 
in (II. ip where the soft matrix element depends on the overlap between the wave 
functions of the initial and the final state particle and therefore on the considered 
process. Second, in (II. ip all effects related to the large scale rrtb are contained in 
the hard coefficient function whereas the non-perturbative input corresponds to a 
soft matrix element of a local operator. The situation is more complicated in (II. 2p 
where all effects from large virtualities are encoded in the hard-scattering kernel, 
but the large scale still enters the non-perturbative matrix elements in form of the 
large energy of the mesons. The distribution amplitudes correspond to non-local 
matrix elements which are defined on the light-cone. The non-locality in position 
space translates into convolutions over the fractions u and v in momentum space as 
illustrated in (II. 2p . 

Example 3: B ^ Dir 

For the first example in QCD Factorization we choose exclusive B decays into heavy- 
light final states as e.g. B Dir. According to the two-particle kinematics in the 
final state, the pion is very energetic in these decays with E.^ = 0{mb) in the B 
meson rest frame. Similar to what we have seen in the last example, the pion can 
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Figure 1.3: QCD Factorization of short- 
and long distance effects in exclusive B 
decays into heavy-light final states. All 
related quantities have already been intro- 
duced in Figure and 



be described by its two particle Fock state in the collinear approximation. The 
factorization formula now becomes (cf. Figure 11.31) 



The soft matrix element Slfip) describes the long-distance dynamics in the B — D 
transition as in (11.11) . The energetic pion enters in form of its distribution ampli- 
tude 07r(^; /^f) and the perturbative effects are contained in a hard-scattering kernel 
T(m; Hp) similar to (11.21) . We see that there is no long-distance interaction between 
the pion and the B — D system at leading power. This corresponds to the famous 
argument of colour transparency which states that the soft gluons cannot resolve 
the fast moving colour- singlet pion [38|39] . 




Let us make one remark concerning strong phases which are important in phe- 
nomenological applications. Strong phases arise from final state interactions, but 
we have just seen that the pion decouples in our example from the B — D system 
in the heavy quark limit m?, oo. The final state interactions are entirely encoded 
in the hard-scattering kernel and therefore predicted to be perturbative. We will 
see in Section 11.21 that a similar argument holds for the case of 5 — tttt where the 
knowledge of strong phases is even more desirable. 

We finally point out that the factorization formula (II. 3p is restricted to the case 
where the D meson picks up the spectator antiquark from the B meson. In the 
opposite case when the spectator goes into the pion, factorization does not hold. 

Example 4: B ^ iriu 

In the last example we finally discuss exclusive charmless B decays. Here we take 
the decay B Tciu as an example which we reconsider in more detail in Section [L3l 
As we deal with a three-body decay in this case, the energy of the pion depends on 
the invariant mass of the lepton pair. We restrict our attention to the case where 
the pion is very energetic in the B meson rest frame with E^^ = 0{mi,). This is 
called the large recoil region which is similar to the situation in the last example. 
However, the factorization formula, which is illustrated in Figure 11.41 turns out to 
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Figure 1.4: QCD Factorization of short- and long distance effects in exclusive 
charmless B decays. New objects enter these decays. The first term contains a 
complicated overlap- contribution ^ whereas the B meson distribution amplitude 
(pB appears in the second one. 

be more complicated and reads 

{'k\J\B) c:^ H{hf) i{HF) + I duj du (f)Biuj; fip) T{uj,u; hf) (j)niu; fip)- (1-4) 

Jo Jo 

It is again restricted to the leading power in Aqcj^i/mf, which we illustrate from now 
on by the symbol "~" for brevity. The second term resembles the factorized form 
in (11.21) . It consists of a perturbative hard-scattering kernel T{uj,u; fip) convoluted 
with the light-cone distribution amplitudes of the pion (l)n{u; fip) and the B meson 
(j)B{u};nF)- Whereas the former already entered the factorization formulas in the 
last two examples, the appearance of the latter is new. At first sight, it might look 
unnatural to consider the B meson on the light-cone. Nevertheless its distribution 
amplitude is a well-defined object, although more complicated and less understood 
than the one of the pion. 

On the other hand the first term in (11. 4p looks like the factorization formula (11.11) . 
but again the situation is more complicated in this case. Whereas the hard effects are 
factorized in a coefficient functions H[ij,p), the remnant matrix element is not a soft 
matrix element as in (II. ip . We therefore wrote ^(/i^) for the overlap-contribution in 
(ll.4p which contains a highly complicated interplay of soft and coUinear interactions. 
A deeper understanding of this overlap-contribution is the main motivation for our 
analysis in Chapter [5l 

Concluding remarks 

Factorization formulas within QCD Factorization are complicated because of the 
presence of soft and coUinear effects at leading power. We have presented a simpler 
example in (ll.3p and a more complicated one in (II. 4p . In general QCD Factorization 
does not prove that a factorization formula is correct. Typically, the factorization 
formula is shown to be valid in the lowest (non-trivial) order of the perturbative 
expansion and then assumed to hold to all orders in a^. It was only the development 
of Soft- CoUinear Effective Theory which provided the necessary tools to formulate 
rigorous factorization proofs. 
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1.1.2 Soft-Collinear Effective Theory 

The origin of Soft-Collinear Effective Theory (SCET) goes back to the beginning of 
this millennium. Since then SCET has sparked the interest of many physicists which 
led to a fast development of this field. If we focus on the publications which have 
contributed to the formulation of the SCET-Lagrangian, we may refer to Bauer et 
al. [221123] . Beneke et al. [24j, Chay et al. |10| and Neubert et al. [H]. 
For our purposes there is no need to introduce the whole concept of SCET as only a 
small part of our calculations in Chapter O is directly related to it. However, SCET 
provides a deeper understanding of QCD Factorization and puts the factorization 
formulas onto robust grounds. We therefore give a brief introduction to SCET 
without going into the technical aspects. 

Preliminaries 

Effective theories in general deal with the strategy that we have already discussed in 
the last section: the idea to disentangle the physics from different length/momentum 
scales. They have restricted validity up to some momentum cut-off and are designed 
to correctly reproduce the IR behavior of a physical process. 

There are two different ways to construct an effective theory. In the bottom-up 
approach either the underlying theory, which is assumed to be valid at all momentum 
scales, is not known or it is not understood how to derive an effective theory from 
the underlying theory. It might nevertheless be possible to write down an effective 
theory motivated mainly by symmetry arguments and inputs from experimental 
observations. Chiral perturbation theory and the Standard Model of particle physics 
fall into this class, as we all hope that the Standard Model loses its validity above 
the TeV-scale. On the other hand in the top-down approach the effective theory can 
be derived directly from the underlying theory. This is case for SCET which stems 
from QCD. 

The same is true for the effective theory which describes heavy mesons, called heavy- 
quark effective theory (HQET) [12113! • We already introduced the phrases hard 
effects {k^ ~ mg) and soft effects {k'^ ~ Aqcd) in the last section, which is the 
appropriate characterization in this case. E.g. the gluon field is split into 

A'^{x)=A^,{x) + A^{x) (1.5) 

and the hard gluons A'^ and the soft gluons A'^ are treated as independent degrees 
of freedom in the effective theory. On the other hand the heavy quark field can be 
decomposed into its large [ip) and small (x) components schematically by 

Q{x) = ^p{x) + xix) with xi^) = C(Aqcd/?^q) ^(x). (1.6) 

This is similar to what is done in non-relativistic theories where this procedure is 
called the Foldy-Wouthuysen transformation [Hj. The small components of the 
heavy quark field and the hard gluons can then be eliminated as degrees of free- 
dom in the effective theory with their effects encoded in coupling constants of local 
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operators. The effective theory is identical to QCD but benefits from the HQE at 
the level of the Lagrangian. It is therefore better suited for the description of heavy 
mesons than QCD and gives rise to the famous heavy quark symmetry [^5l|^ in 
the heavy quark limit rriQ — > oo which had a strong impact on the phenomenology 
of B and D meson decays 

Complications in SCET 

SCET is the effective theory for the description of energetic particles and jets. It 
exploits the hierarchy E ^ Aqcd and follows in principle the same strategy as 
outlined for HQET: introduce power counting in Aqcb/E, identify relevant degrees 
of freedom and integrate out hard fluctuations. However, it turns out that SCET is 
far more complicated than HQET. 

The main reason for these complications lies in the fact that SCET is a non-local 
effective theory in contrast to HQET. We do not go into the details related to 
this point here, but mention that this was the real challenge in the formulation of 
SCET. Notice that we already encountered this non-local nature of SCET in the last 
section in form of the light-cone distribution amplitudes as described at the end of 
example 2. A second complication is related to the fact that there are many different 
degrees of freedom in SCET. During the development of SCET the correct scaling of 
the collinear modes was the origin of some confusion. Later it turned out that there 
are two different scalings depending on the physical process under consideration. 
This is the reason why there are two different versions of SCET usually referred to 
as SCETi and SCETn. Let us have a closer look at this point and at the same time 
introduce our notation and terminology that we use throughout this work. 

Notation and terminology 

The light-cone djTiamics is conveniently formulated with the help of two light- 
like vectors which we denote by n_ and n+ satisfying n\ = and normalized to 
■ n_ = 2. Any momentum can be decomposed according to its projections onto 
these light-cone directions and a two-dimensional transverse plane. We write 



Notice that in our notation k± = k ■ n^. For the purpose of power counting we 
introduce a dimensionless parameter A with = AQCD/f^b- Scaling relations are 
given in the form k ^ {k_, k±, k^) where we drop factors of which can be restored 
easily by a dimensional analysis. For example, hard momenta with ~ and 
soft momenta with ~ Aqcd now become ph ~ (1, 1, 1) and Ps ~ (A^, A^, A^). 

In charmless B decays, the energy of the fast moving light hadrons typically scales 
as ~ rrih in the B meson rest frame (cf. the examples 3-4 in the last section). 
There are two possible scalings of collinear momenta given by phc ~ (1, A, A^) and 
Pc ~ (1, A^, A^), the first entry reflecting the large energy in each case. These are the 




with 



k±- n± = 0. 



(1.7) 
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iermmoiogy 


Momentum scaling 


Virtuality 


hard 


(1,1,1) 


ml 


hard-collinear 


(1,A,A2) 




soft 


(A^A^A2) 


^^QCD 


coUinear 


(1,A2,A4) 


^^QCD 



Table 1.1: Terminology for the relevant momentum regions in charmless B 
decays. The momentum scaling corresponds to the light- cone decomposition 
k ~ {k-, k±,k+) with a dimensionless parameter = KqcDl^h, see text. We 
have restored the mass dimension in the respective virtualities for convenience. 



two different scalings of collinear fields that we mentioned above. The important 
difference lies in the respective virtualities ~ mfeAqcD ^ pI andp^ ~ ^qcd ~ pI- 
We see that the former, to which we refer as hard-collinear from now on, introduce a 
new scale fihc ~ ("^bAqco)^^^, whereas the latter, simply called coUinear, are related 
to the hadronic scale Aqcd as usual. The appearance of the hard-collinear scale is 
another complication in SCET compared to HQET. In the following we treat this 
intermediate scale as a perturbative scale with asifJ^hc) ~ 0.4 

We have summarized our terminology in Table II. 1[ Depending on the considered 
process some of these modes turn out to be irrelevant. Inclusive decays as e.g. the 
endpoint spectrum in B ^ Xg'y can be described with an effective theory called 
SCETi which contains hard, hard-collinear and soft modes. The invariant mass 
of the energetic jet is typically Mx ~ mb Aqqd which explains the relevance of 
the hard-collinear modes in this case. The effective theory for exclusive decays as 
e.g. B Dtt is SCETn containing hard, collinear and soft degrees of freedom. In 
this case, the virtuality of the collinear modes is related to the mass of the energetic 
pion with ~ Aqqj). Hard-collinear modes may appear in exclusive decays as well 
induced by soft-collinear interactions, which makes a two-step matching procedure 
QCD SCETi — s> SCETn necessary. This is the case e.g. in B ^ Tciu but not 
in B —> Dn where these effects turn out to be power suppressed. This explains 
why the factorization formula for the former is much more complex than that of the 
latter as we have seen in the examples 3-4 in the last section. 

QCD Factorization versus SCET 

What is the difference between QCD Factorization and SCET? Both frameworks 
provide a systematic description of QCD in the heavy quark limit mj, oo, i.e. there 
is none: they are equivalent! There has been some confusion about this point in the 
literature. It seems as if there were different predictions from QCD Factorization 
and SCET in hadronic two-body decays of B mesons. We stress that this is not 
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related to the underlying frameworks but to different treatments of power corrections 
and input parameters of two groups: BBNS and Bauer, Pirjol, Rothstein, Stewart 
(BPRS). Therefore, the predictions from BBNS and BPRS differ although the ones 
from QCD Factorization and SCET do not. More details concerning this issue can 
be found in ^7\M\- 

From the conceptual point of view, QCD Factorization and SCET are different. 
QCD Factorization relies on an explicit analysis of momentum regions of Feynman 
diagrams. In contrast to this, SCET is a rigorous theory in the sense that it is 
derived from an established theory (QCD) in a well-defined expansion (Aqcd/E) on 
the level of the Lagrangian. With the Lagrangian and the respective Feynman rules 
at hand, it is a convenient tool to address a large variety of inclusive and exclusive 
processes with applications going even beyond the domain of B physics. 

1.1.3 Alternative approaches 

For completeness we give an overview of alternative approaches to exclusive charm- 
less B decays. As this section is not directly related to the remainder of this thesis, 
our presentation will be brief and concentrate on the most prominent methods. 

Lattice Gauge Theory 

The first method which comes into mind may be lattice field theory. In recent 
years there has been considerable progress concerning the reduction of its systematic 
uncertainties. Most importantly, a steadily increasing number of unquenched lattice 
calculations becomes available reducing the magnitude of its uncertainties in some 
cases to the percent level (in particular in ratios of hadronic quantities). For recent 
reviews on lattice results for heavy quark systems we refer to [SUf^ . 
Unfortunately, lattice gauge theory can tell us very little about exclusive charmless 
B decays. The main difficulty is related to the implementation of energetic mesons in 
the lattice calculations (which are usually performed in the B meson rest frame). So 
far, there are reliable lattice results for B meson decay constants and heavy-to-light 
form factors at small recoil which can be extrapolated to the large recoil case using 
dispersion relations. There are some interesting considerations to directly access the 
large recoil region in an approach called moving NRQCD [52j. However, if it will be 
possible to address hadronic B decays on the lattice remains very challenging. 

QCD Sum Rules 

Although intrinsically limited in their accuracy, QCD sum rules have become a 
serious competitor to lattice gauge theory. In some sense it can even be considered 
as complementary to QCD Factorization since it provides important information 
about the non-perturbative input parameters of the latter in form of decay constants, 
heavy-to-light form factors (at large recoil) and light-cone distribution amplitudes. 
For a selection of recent Light-Cone Sum Rule (LCSR) results we refer to [53)151] . 
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Charmless non-leptonic B decays have been studied in the LCSR approach [551458] . 
It should be noted that this requires a substantial extension of the standard LCSR 
formalism [55]. The outcome of the LCSR analysis is in good overall agreement with 
the QCD Factorization prediction. This technique is sometimes used to address the 
importance of power corrections in the QCD Factorization framework. This issue 
should, however, be treated with care as the essential QCD dynamics is considered 
at finite rrih in this approach and the heavy-quark limit — > oo is only taken at the 
very end of the calculation after performing the continuum subtraction. It has been 
pointed out in ^\ that one should rather consider LCSRs within SCET in order to 
properly address distinct contributions of the QCD Factorization framework. 

Perturbative QCD 

Perturbative QCD (pQCD) [60f6T] . which is also known as fc^-factorization, is based 
on a hard-scattering approach and may at first sight look similar to QCD Factoriza- 
tion since decay amplitudes are expressed as convolutions of hard-scattering kernels 
with meson wave functions. However, pQCD does not share the same systematics 
as QCD Factorization and SCET and relies on the assumption that soft contri- 
butions to heavy-to-light form factors are suppressed by Sudakov effects (criticism 
concerning this point has been raised in [62]). The form factor is thus considered 
to be dominated by hard gluon exchange and the wave functions have to include a 
dependence on transverse momenta in order to avoid endpoint-singularities in the 
convolution integrals. 

As a consequence pQCD is very different from QCD Factorization. This is reflected 
by a different hierarchy of various contributions to the hadronic matrix elements, 
different non-perturbative input quantities and by the fact that pQCD does not 
recover naive factorization (cf. next section) in any limit. Only recently the authors 
of pQCD have included next-to-leading order corrections which have simply been 
taken over from the QCD Factorization and SCET analyses [63fM] . 

Phenomenological approaches 

Apart from these dynamical approaches there exist many ideas how to extract the 
interesting CKM information in exclusive charmless B decays without explicitly cal- 
culating the hadronic dynamics. The basic idea is to find a suitable parametrization 
of the decay amplitudes and to use symmetry arguments (e.g. isospin, SU{3)) to de- 
rive relations between different decay processes. A sufficiently small set of unknown 
amplitudes is finally fitted to experimental data. 

These approaches should be seen as complementary to the dynamical ones discussed 
above. On the one hand they require input from the dynamical approaches in order 
to control their own intrinsic uncertainties of e.g. SU (3)-breaking, on the other they 
can serve as a guide for the dynamical approaches and indicate where these fail to 
correctly reproduce the data. Rather than presenting some exemplary ideas here, 
we refer to a recent review about the phenomenology of B decays [03] . 



22 



1. EXCLUSIVE CHARMLESS B DECAYS 



1.2 Hadronic two-body decays 

After this general introduction we come to the first class of exclusive B decays which 
we consider in detail in this work: hadronic two-body decays. There is a very large 
variety of decay channels which fall into this class as e.g. B tttt, B Ktt, B pp 
or Bs —>■ TiTT. From the phenomenological viewpoint, we are mainly interested in the 
underlying weak interactions which differ among these decays due to the flavour 
contents of the mesons. In our analysis concerning the QCD dynamics, the only 
differences lie in the pseudoscalar (P) or vector {V) nature of the light mesons 
and small SU{3) breaking effects. As we do not focus on the phenomenological 
implications of these decays here, it will be sufficient to concentrate on the B —>■ nn 
channels. 



1.2.1 Preliminaries 

B meson decays are mediated by weak interactions. As the typical energy and mo- 
mentum scales in these decays are much smaller than the mass of the Vl^-boson, we 
may work with an effective weak Hamiltonian which consists of a sum of local op- 
erators Qi multiplied by short-distance coefficients Ci and products of CKM matrix 
elements Xp = VpbV*^. The effective Hamiltonian is given by 

^<^s = ^Yl ( ^1^1 + ^2^2 + J^^Qi + CsQsj + h.c, (1.8) 

p=u,c \ i=3 / 

where are the left-handed current-current operators, Qs-e are the QCD penguin 
operators and Qs is the chromomagnetic dipole operator. Their explicit form reads 

Ql = ipb)v-Aidp)v-A, Ql = {Pibj)v-A{djPi)v-A, 



Q3 = {dh)v-A Y.qiQQ)v-A, Qi = {dibj)v-A Y^qim 



i)V~A, 

Q5 = {db)v-A J2giQ(l)v+A, Qe = {dibj)v-A T.qiQjQi)v+A, 

Qs = -|^mfe(ia^,(l + 75)6-^^6, (1.9) 

where {qiq2)v±A = ^i7aj(1 =t 75)^2 and the sum runs over all active quark flavours 
in the effective theory, i.e. q = u,d,s,c,b. If no colour index i,j is given, the two 
operators are assumed to be in a colour singlet state. The definition of the dipole 
operator Qs corresponds to the sign convention iDf^ = id^ + gsT^A^. 
In principle, more operators have to be taken into account. The complete set of op- 
erators contains in addition electroweak penguin operators and the electromagnetic 
dipole operator which are important for B — > Kir decays. As we focus on the 
B ^ 7171 channels here, the effects from these operators can safely be neglected. 
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Figure 1.5: Topological amplitudes ai_4(MiM2) as introduced in 166^ . The 
quark to the right and the spectator antiquark with flavour Qs (not drawn) form 
the meson Mi with flavour content [qsu], [qsd], [qsd] and [qsq], respectively. 
The two upper lines form the meson M2 with respective flavour [ud] , [uu] , [qq] 
and [qd\ where q £ {u, d, s}. The flavour- singlet penguin amplitude a3(MiM2) 
does not contribute to the B vrvr decay amplitudes. 

For the parametrization of the B —>■ nTt decay amphtudes we follow the notation 
of j66]. Neglecting some smaller amplitudes related to weak annihilation, they can 
be written as 

(7r+7r-|Heff|5°) = |A„[ai(7r7r) + a^(7r7r)] + Aca^(7r7r)} A^^, 

- (7rV|Kff|fi°) = {XuIMtiti) - atinn)] - A,a^(7r7r)} A^^. (1.10) 

The aj(7r7r) are called topological amplitudes and are related to the flavour flows in 
the decays. More precisely, ai(7r7r) is the colour- allowed tree amphtude, a2(7rvr) is 
the colour-suppressed tree amplitude and aKim) is the QCD penguin amplitude as 
illustrated in Figure II. 5[ We see that B~ 7r~7r° is a pure tree decay within the 
approximations mentioned above. It is free of interference effects which typically 
introduce large uncertainties due to our poor knowledge of the weak phases and 
therefore particularly suited to test the QCD Factorization predictions. On the 
other hand tt^tt^ is a colour-suppressed decay as can be seen in (11.101) by 

the absence of ai(7r7r). Any prediction concerning this decay is expected to be 
accompanied by large uncertainties. 

For later convenience it is useful to write the normalization of the amplitudes as 

A.. = z^m|F^--(0)A, (1.11) 

where F^^^(O) is a transition form factor at maximum recoil and f.^ is the pion 
decay constant. Their precise deflnitions read 

{n^{p')\u{^-f)b\B\p))\^,^_^,^,^^ ^ m|F+^-(0), 

{n-iq)\d^,^,u\0) = -tf^q^. (1.12) 

As in the examples that we discussed at the beginning of this chapter, the light-cone 
distribution amplitudes of the pion and the B meson enter the factorization formula 
for hadronic two-body decays. Their explicit deflnitions will be given in Chapter [5] 
(for light mesons the deflnition can be found in (15.391) and for the B meson in (I5.46p ). 
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1.2.2 Factorization Formula 

According to the QCD Factorization framework [T9H2T], the hadronic matrix ele- 
ments of the operators in the effective weak Hamiltonian take the form 



where Mi denotes the hght meson which picks up the spectator quark of the B 
meson and M2 is a second hght meson. The factorization formula is illustrated in 
Figure II. 61 Let us compare this factorization formula with the examples that we 
discussed in the beginning of this chapter. First, we notice that the factorization 
formula is again restricted to the leading power in the HQE in AqcB/^b- Further, 
we have suppressed the dependence on the factorization scale fip for simplicity here. 
The perturbative information in fll.l3p is encoded in the hard-scattering kernels 
and T^^ which also appeared in the examples 2-4. The same is true for the 
light-cone distribution amplitudes of the light mesons 0Mi , 4>M2 whereas the one of 
the B meson (pB entered the factorization formula in example 4. The distribution 
amplitudes always come in combination with decay constants /b, /a/i, fA-h which 
have been suppressed in Section II. 1.11 for simplicity. 

The last piece in fll.131) is the heavy-to-light form factor p^~'^'^^ at maximum recoil. 
The form factor still contains perturbative and non-perturbative effects and there- 
fore does not fit into the pattern that we developed in Section fl.l. 11 The form factor 
itself can be factorized as we discuss in Section [L3l but this is not needed here. We 
follow the BBNS analysis which treats the form factor as an external input with its 
numerical value taken from a Light-Cone QCD Sum Rule calculation. 

In leading order (LO) of the perturbative expansion, the factorization formula sim- 
plifies tremendously since = const -|- 0{as) and T^^ = 0{as). The second term 




(1.13) 



+ Ib fhh fM2 / dudvdu Ti {uj,v,u) 0b (t^) 4>Mi{v) (pM^zi^) 





/ 



/ 




+ 





Figure 1.6: QCD Factorization of short- and long distance effects in hadronic 
two-body decays. The former are contained in perturbative hard- scattering 
kernels , T^^ , the latter in light-cone distribution amplitudes (f>B, <t>Mi, <t'M2- 
The heavy-to-light form factor F is not factorized in the BBNS approach. 
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in fll.lSp is thus absent in this case and the convolution in the first term simply gives 
the normalization of the distribution amplitude. This yields 

{M,MM\B) ~ F^-^^nO) /a/,. (1.14) 

This approximation corresponds to naive factorization. With naive factorization we 
mean that a hadronic matrix element of a local four-quark operator is split into two 
matrix elements of bilinear quark currents, as e.g. in 

{MiM^mB) ~ {M^\{ub)v^A\B) {M2\{du)v-A\0) (1.15) 

with the first matrix element on the right-hand side giving the form factor and the 
second one the decay constant as in fll.l4p . Naive factorization was used before QCD 
Factorization was established and gave fairly good predictions which was surprising 
at that time [BTfUS] . We now understand why this is the case: Naive factorization 
corresponds to the leading term in the combined expansion in as and Aqco/'^b- 
The fact that QCD Factorization reproduces naive factorization in this expansion is 
non-trivial. 

One obvious problem with naive factorization is that it does not give rise to strong 
rescattering phases because of the lack of final state interactions. In fll.l3p we see 
how this problem is cured in QCD Factorization. At leading power all final state 
interactions are encoded in the hard-scattering kernels and T^^. The strong 
phases are thus predicted to be perturbative in the QCD Factorization framework. 
Moreover, we just have seen that strong phases are absent in LO of the perturba- 
tive expansion. To summarize, strong phases are of 0{as) and 0(Aqcd/"^6) i-e. 
generically small in QCD Factorization. 

We conclude this section with a comment on the factorization proof of fll.lSp . The 
explicit calculation in P^] - BT| showed that the factorization formula holds to 0{as)- 
A recent analysis [23H2Z] in combination with our calculations in Chapter [S] and H] 
extends this explicit proof to 0{al). So far a rigorous proof to all orders in as is 
still missing. For the similar (but simpler) case of i? — Dtt, which we introduced 
in example 3 of Section 11.1. 11 a factorization proof has been formulated in SCET 
and can be found in [69] . A first step towards an all-order proof of (11.131) has been 
undertaken in |70j . 

1.2.3 Perturbative corrections 

The power of the factorization formula f ll.l3p lies in the fact that it allows for a 
systematic calculation of perturbative corrections. They are contained in the hard- 
scattering kernels and T^^ which describe two different mechanisms to which we 
refer as vertex corrections and spectator interactions, respectively. Whenever the 
spectator antiquark in the B meson enters the perturbative subgraph, this contribu- 
tion is assigned to T^^ otherwise to T^. We now present the status of the perturbative 
calculation and identify the contributions that we address in Chapter [3] and HI 
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Figure 1.7: Tree level diagram. The line to the 
left (right) of the vertex denotes the b quark (light 
quark which goes into Mi), the upper lines the light 
quark/ antiquark which form M2. When the spectator 
does not participate in the scattering, it is not drawn. 

Tree Level 

We have already anticipated in the last section that QCD Factorization reproduces 
naive factorization in LO of the perturbative expansion. There is only one diagram 
that contributes at this order which is depicted in Figure II. 7[ As the spectator 
antiquark does not participate in this case, this corresponds to a contribution to T^. 
We quote the result in form of the topological amplitudes introduced in f ll.lOp 

af(MiM2) = a + ^ + 0(a,), ^ = l,2,4 (1.16) 

where the upper (lower) signs apply when i is odd (even) and the superscript p is 
to be omitted for i = 1,2. To illustrate the phrases colour- allowed (ai) and colour- 
suppressed (0:2), we have a look at the numerical values of the Wilson coefficients 
Ci(mb) ~ 1.1 and C2{mb) ~ —0.2 giving ai ~ 1.0 and 02 ~ 0.2. Notice that 0:2 is 
particularly small due to a cancellation in the two terms. 




Next-to-leading order 

The next-to-leading order (NLO) corrections have been calculated by BBNS [T9] ~ 
[21] . They consist in the calculation of the (naively) non-factorizable diagrams in 
Figure 11.81 whereas the (naively) factorizable diagrams in Figure 11.91 turn out to be 
irrelevant. The terminology is such that the former contain interactions between the 
B — Ml system and M2 whereas the latter do not. An explanation for the fact that 
only non-factorizable diagrams have to be considered in this context is relegated to 




Figure 1.8: Non-factorizable diagrams in NLO: Vertex corrections (a), hard 
spectator interactions (b) and penguin contractions (c). The last diagram in- 
volves an insertion of the dipole operator Qs- 
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Section I3.3.2[ In addition to the vertex corrections in Figure 11.8b and the spectator 
interactions in Figure flTSb . the penguin contractions in Figure [TTHb have to be taken 
into account for the calculation of the penguin amplitude a^. As the spectator is 
not involved in these diagrams, the respective contribution is again assigned to T^. 

At NLO the general form of the topological amplitudes ai, a2 and 04 becomes 



2\ 



;i.i7) 



where the functions V^^^if ^^^Pp^'' stem from the diagrams in Figure ll.8b /b/c. 
respectively. They contain convolutions of the hard-scattering kernels with light- 
cone distribution amplitudes according to the factorization formula (11.131) . The 
respective arguments V^^^{M2), H^^^{MiM2) and Pp^\M2) have been suppressed in 
(11.171) for simplicity. The normalization of the hard spectator interactions reads 

NH = 47r^ ^-^^t^ . (1.18) 

In contrast to the remnant NLO diagrams, the spectator interactions correspond to 
tree level diagrams at this order which is the origin of the factor Att"^. The ratio of 
hadronic quantities involves the quantity = O(Aqcd) which is related to the B 
meson distribution amplitude (the definition can be found in (15.561) ). Notice that 
the normalization scales as ~ 1 in the HQE and in the large A^^c-expansion. 
One might therefore expect the spectator interactions to dominate over the vertex 
corrections and the penguin contributions because of the factor 47r^ ~ 40. However, 
for realistic values of the hadronic quantities this is not the Nh ~ 0.5-1.5. 



In the following we quote the explicit NLO results for the vertex corrections V^^^ 
and the spectator interactions H^^\ As our calculations in Chapter [3] and H] are 
related to the topological tree amplitudes, the penguin contributions Pp^^ are not 
needed for our purposes (they can be found in [BBJ). The NLO vertex corrections 
can be written in the form 

y(i)(M2) = -61n^ - 18+ / dug{u) 4>M,{u). (1.19) 

"^6 Jo 
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The scale dependence exactly matches the one of the Wilson coefficients in the LO 
terms in f ll.l7p as desired. The function g{u) is found to be 



, , 3(1 -2m), 

g[U) = z m M — Stir + 

u 



2 1 — 3n 



2Li2 (m) — In u z — In M — 2i7r In M — (m ^ m) 

u 



;l20) 



where we wrote u = 1 — u. The appearance of an imaginary part in g{u) is particu- 
larly interesting since it is the origin of a strong phase shift between the final state 
mesons. The penguin contributions Pp^'* are another source of an imaginary part 
whereas the tree-level spectator scattering gives a real contribution which reads 

if«(M,M2) = Hdu Cdu fdv ^ ^ (L21) 

Jo Jo Jo UJ V u 

Let us make one remark concerning the typical scales of the perturbative correc- 
tions. Following our terminology in Table II. ![ the perturbative effects in stem 
from hard modes {k"^ ~ ml) whereas T^^ contains hard as well as hard-collinear 
effects (fc^ ~ fi\^ ~ mbAqcD)- The tree level spectator scattering in Figure [TTSb is 
in particular most naturally associated to the hard-collinear scale which enhances 
this contribution by a factor of as{f^hc) / dsimb) ~ 2 with respect to the other contri- 
butions. This issue is hidden in our notation in f ll.l7p . A correct treatment of the 
perturbative scales favours a SCET analysis of the spectator interactions in order 
to resum large logarithms In /i/m^ into short-distance coefficient functions. 



Next-to-next-to-leading order 

The next-to-next-to-leading order (NNLO) calculation is to date incomplete. It has 
ffist been addressed in the so-called large Po-limit in [7T|7^ . In this approximation, 
only a very small subset of NNLO diagrams has to be calculated, namely the ones 
with massless fermion- loops which give rise to a factor Uf. At the end of the calcu- 
lation, the fairly motivated substitution rif —>■ — 3/3o/2 is made in the hope to catch 
the main contribution of the full calculation in this way. It must be admitted that 
the large /3o-limit sometimes gives a good approximation. However, whether this is 
the case or not can only be judged a posteriori when the full result is known. Con- 
cerning the parts of our NNLO calculation, we will answer this question in Chapter [3] 
andU 




Figure LIO: Sample of NNLO diagrams: Vertex corrections (a) and spectator 
interactions (b). 
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Up to the work in [73], which also considered a small and incomplete subset of NNLO 
diagrams, the penguin amplitude has not yet been addressed. We therefore focus 
on the tree amplitudes in the following and disregard all diagrams with penguin 
contractions as well as insertions of the operators Qs-q and Qs- We extend the 
general form of the tree amplitudes ai and 0^2 in NNLO to 



ai{MiM2) = Ci + — — + 



;i.22) 



+ 



QV^ + C,±iVr + Nh ( C,H'^> + Ciiijf ) + 0{at). 



A first systematic study of the NNLO corrections addressed the hard spectator in- 
teractions H{ 2 which were calculated within SCET in [25]. A subsequent SCET 
calculation can be found in [26j. There is a tiny discrepancjQ between both calcula- 
tions which is claimed to be of minor numerical importance [26]. A third independent 
calculation of the spectator interactions in full QCD may help to resolve the ori- 
gin of this discrepancy [27]. The QCD calculation consists in the computation of 
~ 50 1-loop diagrams as those depicted in Figure II. 10b with its main complication 
due to the presence of two perturbative scales Hhc and m^. This problem is circum- 
vented in SCET where the two scales are systematically disentangled in the two-step 
matching procedure QCD — > SCETi SCETn. Another advantage of the SCET 
calculation lies in the fact that the hard-coUinear effects in T^^ are exactly the same 
as those entering the B —>■ Mi form factor [17]. As they were already known at the 
1-loop level from a different analysis [7^1475] . Beneke/Jager and Kivel calculated the 
missing hard contributions [25f26] . 

(2) 

No such complication arises in the calculation of the vertex corrections which 
we address in detail in this work. However, the vertex corrections represent a chal- 
lenging 2-loop calculation including ~ 80 diagrams as the ones in Figure II. 10b . We 
divide the calculation into two parts: We first focus on the imaginary part of 2 
in Chapter O and calculate the real part subsequently in Chapter HI This is done 
for several reasons. First of all, each part is independently well-defined and repre- 
sents a lengthy calculation on its own. It is natural to focus on a simpler part of 
this highly demanding calculation in a first step. Second, the imaginary part enters 
the topological amplitudes at 0{as) which means that the imaginary part of V}^2 
has the complexity of an effective NLO calculation. Conceptual complications that 
arise in NNLO can therefore be relegated to our analysis in Chapter |H Finally, 
the knowledge of the imaginary part is particularly important in phenomenological 
applications as it is the origin of a strong phase shift between the final state mesons. 



One important remark is in order concerning the general form fll.22p of the tree 
amplitudes in NNLO. The reason for this compact form which implies ai ^ 0:2 by 
the exchange Ci ^ C2 is related to Fierz symmetry arguments. A Fierz reorder- 
ing of a four-quark operator reshuffles the contracted fields and thus the related 



^The corrected version of |2S] is in agreement with 
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flavours. Notice that ai and 0:2 interchange their roles under the exchange of the 
quark flavours u ^ d as shown in Figure ll.Si However, it is no longer obvious if 
Fierz symmetry is preserved in the effective theory at the level of the Wilson co- 
efficients. The NLO calculation showed the Fierz symmetric form fl 1.171) to hold 

(2) 

and we therefore expect the same for the imaginary part of 2 which has NLO 
complexity. In their SCET analysis of the spectator interactions, Beneke and Jager 
paid special attention to this point and found that their result can indeed be written 

in the Fierz symmetric form fll.22p . We have not yet solved this question for the 

(2) 

real part of which requires a Fierz symmetric definition of evanescent operators 
at the 2-loop level. As a consequence, we can only give a result for the real part 
of the tree amplitude ai in Chapter H] but so far we cannot relate it to that of 02 • 
More details concerning Fierz symmetry and evanescent operators can be found in 
Section 13. 1[ 

1.2.4 Power corrections 

The factorization formula f ll.l3p is exact in the formal heavy quark limit nib 00. 
Unfortunately, power corrections do not factorize in general and cannot be calcu- 
lated systematically. They represent the main limitation of QCD Factorization with 
a generic size of (9(Aqcd/'"^6) ~ 5 — 15%. The reader may wonder why we perform 
highly non-trivial NNLO calculations that are accompanied by substantial system- 
atic uncertainties. Our ultimate goal is to look for deviations between experimental 
data and the QCD Factorization prediction and to address them to effects from New 
Physics rather than to presumably sizeable power corrections. 

From our point of view QCD Factorization is nevertheless the most convincing theory 
to study non-leptonic B decays and we should therefore reduce all uncertainties as 
much as possible. Calculable higher order perturbative corrections can also turn 
out to be sizeable and it would be erroneous to associate them with unknown power 
corrections. The rich phenomenology and the interplay with alternative approaches 
to non-leptonic B decays may help us to test QCD Factorization and to get an 
estimate of the size of power corrections on the one hand and to focus on the 
interesting CP violating observables on the other. Rather than aiming at a perfect 
global description of all decay channels, we are convinced that we can achieve the 
necessary precision in selected observables. We also refer to our discussion after 
fll.lOp in this context, concerning the question why some decay channels are better 
suited for a QCD Factorization prediction than others. 

Obviously, it would be a significant improvement if we could get some handle about 
power corrections. In the remainder of this section we briefiy recapitulate how 
BBNS implement certain classes of power corrections. The price to pay for this 
procedure is that their predictions become model-dependent in this way. For a 
model-independent analysis we refer to [77] which contains a formal classification of 
power corrections according to factorizable and non-factorizable operators in SCET. 
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Chirally enhanced contributions 

One class of power corrections, which is related to the projection on higher twist dis- 
tribution amplitudes of the light mesons, appears to be enhanced by large numerical 
coefficients. These corrections are typically proportional to 

^^^^^ mb{fi){mu + md){fi)' (1-23) 

which is formally of (9(Aqcd/'"^6) but numerically close to unity. Fortunately, an 
important part of this contribution, the so-called scalar penguins, turn out to be 
factorizable (at least to NLO) and can be calculated systematically within QCD 
Factorization. This results in an additional contribution to the penguin amplitude 
^4 which is currently known to NLO and can be found in |66j . 



In contrast to this the twist-3 projections related to the spectator interactions in Fig- 
ure [L8b do not factorize which can be seen as follows. In this case, the contributions 
to «!, a2 and can be written as 

O^sCf ^ M,, f' r. r 0bM <Pm,{v) (PaM 



/'OO pL pL 

5a^,{M,M2) = -^C,^,NHr';;'XB / dcu du dv 
47riVc Jo Jo Jo 



UJ V u 

(1.24) 



which involves the twist-3 distribution amplitude of the meson Mi with its asymp- 
totic form (prniiv) = 1. As it does not vanish at f = 1, the contribution in (11.241) 
is divergent which states that factorization does not hold for this power-suppressed 
contribution. BBNS propose a parametrization of the logarithmically divergent in- 
tegral in the form 

' dv ^Xh - XH=il + PH e^'") In ^ (1.25) 

as the integral is expected to be regulated by a hadronic scale = (9(Aqcd)- In 
addition, they allow for a complex coefficient with pn = 0{1) and an arbitrary 
phase (pH which might be generated due to soft rescattering effects. Notice that Xh 
is treated to be universal, i.e it does not depend on the meson Mi and is assumed 
to be the same for all topological amplitudes in (11.241) . In our analysis of the tree 
amplitudes in Chapter [S] and HI we adopt the BBNS treatment of non-factorizable 
chirally enhanced contributions. 



Weak annihilation 

We briefly comment on a second class of power-suppressed effects: weak annihilation. 
The annihilation diagrams from Figure [T.lll turn out to be suppressed by one power 
in Aqcd/'"^6- They exhibit similar endpoint divergences as the chirally enhanced 
contributions and are therefore non-factorizable. In the BBNS approach, they are 
parameterized in analogy to (I1.25P by 

XA={l+PAe'^^) In^. (1.26) 
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Beneke and Neubert make a more relaxed assumption on the universality of anni- 
hilation contributions in [66] and examine a specific scenario with three different 
phases (pA for the final states PP, PV and VV. 

Our final remark is related to recent investigations of the annihilation amplitudes. 
They have been considered in [78j using a new type of factorization formula which 



includes zero-bin subtractions in SCETh [7^]. These subtractions render the convo- 
lution integrals finite which implies that the annihilation amplitudes become calcu- 
lable. We will come back to the formal aspects of this approach in Section [1.3.21 and 
focus on the outcome for the annihilation amplitudes here. Strong phases from an- 
nihilation are found in [78] to be of C(a^(/i/ic)AQCD/"^6) and are therefore expected 
to be small in qualitative agreement with predictions from a Light-Cone QCD Sum 
Rule analysis [HH]- However, whether or not this new type of factorization formula 
represents a model-independent prediction for the annihilation amplitudes is still a 
matter of debate. 
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1.3 Heavy- to- light form factors 

Our calculation in Chapter [5] is related to an analysis of heavy-to-liglit form factors. 
These form factors encode the strong interaction effects in exclusive semi-leptonic 
decays (e.g. B niu), hadronic two-body decays (e.g. B — > tttt) and exclusive 
radiative decays (e.g. B — > K*'y). Rather than aiming at a quantitative description, 
our analysis focuses on the conceptual aspects concerning factorization. Heavy-to- 
light form factors are among the "simplest" objects for studying the complicated 
QCD dynamics in exclusive charmless B decays. 

1.3.1 Factorization Formula 

We concentrate on the three independent form factors in i? — vr transitions which 
can be defined as 

{nip') I qYb I B{p)) = F^q') {p^ + p' ^ + F^{q') q', 
W) I q^y'^'q.h I B{p)) = lll^ [q\p^ + p' ^ - - O ?1 (1-27) 

with q = p — p'. We are interested in the large recoil region ^ ml where the pion 
is very energetic with = 0{mb) in the B meson rest frame. The factorization 
formula has already been introduced in example 4 at the beginning of this chapter. 
In a slightly different notation, it reads 

/•oo /•! 

F,(g2) ~ iJ,(g2) ^(^2) + / du du 0B(^)T,(a;,n;g2) 0^(u), (1.28) 

Jo Jo 

where we made the kinematical dependence on the momentum transfer explicit and 
suppressed the one on the factorization scale. We recapitulate its interpretation 
in terms of our SCET terminology from Table 11.11 The perturbative information 
is contained in the hard coefficient function Hi{q'^) and the hard-scattering kernel 
Ti{u,uj; q"^); the former including hard effects (~ m-l), the latter hard and hard- 
collinear effects (~ mbAQco)- The light-cone distribution amplitudes 4>b{^) and 
(j)Tr{u) encode soft and collinear effects, respectively, from the scale Aqcd- Finally, 
the overlap-contribution ^(g^) is poorly understood so far and the main subject of 
our analysis. The important point to notice is that the same function ^(g^) enters the 
three B ^ n form factors (there are two other overlap functions C.wiq'^) and ^±(g^) for 
the seven B ^ V form factors). This is the basis of approximate symmetry relations 
between different form factors [80] which are broken by perturbative corrections 
and power corrections. Even if we consider the overlap-function ^(g^) as unknown, 
the factorization formula fll.28p represents a useful simplification as the number 
of unknown hadronic quantities has been reduced from three to one (counting the 
distribution amplitudes as known hadronic quantities in this context as they are 
universal quantities that can be analyzed in different processes). 
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The structure of the factorization formula has first been seen to emerge in [81]. A 
rigorous factorization proof in SCET has been formulated later in [821481] . The proof 
is usually performed in two steps, matching first QCD —>■ SCETi at the hard scale 
and subsequently SCETj SCETn at the hard-collinear scale. The two terms 
in (11.281) are related to two different SCETi currents which contribute in leading 
power to the form factors. What remains to be shown is that one of them exhibits 
the symmetry properties mentioned above and that the other gives rise to a finite 
convolution integral of leading twist distribution amplitudes as in the second term 
of the factorization formula (11.281) (more details will be given in Section [F!^ . 

We finally remark that the overlap-contribution ^(g^) still contains hard-collinear 
effects as only hard effects have been factorized in the first matching step into Hi{q^). 
As long as we treat the hard-collinear scale (mbAqcD)^^^ as a perturbative scale, the 
factorization of short- and long-distance effects is thus incomplete. 

1.3.2 Closer look at ^{q^) 

We follow the notation of [H3] by defining the overlap-contribution as a matrix 
element of a SCETj current in the form 

(7r(p ) I {ihcWhc)K I B{p)) = 2E^ e(g'), (1-29) 

where denote the large components of the heavy quark field satisfying i/h^ = hy 
and ^hc the ones of the hard-collinear field with fi-^hc = 0. The definition of the 
hard-collinear Wilson line Whc can be found in [53] . 

The attempt to factorize ^{q'^) in SCETn into light-cone distribution amplitudes 
leads to several complications: Sub- leading twist projections and three-particle Fock 
states of the B meson and the pion are found to contribute at leading power. Even 
worse, some convolutions turn out to be divergent at the endpoints which renders 
the overlap-contribution incalculable similar to what we have seen in our discussion 
about power corrections in the context of non-leptonic decays in Section ll.2.4[ 

The non-factorization of soft and collinear effects in SCETn and the related question 
about the correct treatment and interpretation of endpoint singularities is currently 
not completely understood. The problem does not only appear in energetic -B — >■ vr 
transitions but also in the perturbative description of the pion form factor. In this 
case the soft-overlap contribution {Feynman mechanism) appears at sub-leading 
power. A solution within the effective theory approach would be a significant break- 
through on the conceptual level. Moreover, if the soft-overlap contribution could 
be expressed in terms of a few fundamental hadronic parameters, this would imme- 
diately lead to many interesting applications. This is to some extent realized in a 
recent approach from Manohar and Stewart [79j which we will discuss below. In any 
case, the problem of endpoint-divergences in SCET deserves further study which is 
the main motivation for our analysis in Chapter [5l 
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Let us briefly comment on two analyses of tlie overlap-contribution ^(g^) in SCET. 
Tlie first one has been performed by Neubert et al. [H] and is based on an extended 
formulation of SCETn [H5] which includes additional degrees of freedom, so-called 
soft-collinear messenger modes, with momentum scaling kg^ ~ (A^,A^,A^). These 
modes allow for a cross talk between the soft and the coUinear sector of the theory 
as they can be exchanged between soft and coUinear fields without bringing them far 
off-shell. Technically, these modes serve as a regulator of the divergent convolution 
integrals which leads to the following conclusion: If matrix elements involving mes- 
senger fields contribute at leading power to a hadronic matrix element, factorization 
into soft and coUinear effects is spoilt. 

While messenger modes provide a technical solution to calculate otherwise ill-defined 
quantities in perturbation theory, their physical interpretation in non-perturbative 
matrix elements remains obscure as virtualities with k^^ ~ Aq(2;£,/mh are below the 
confinement scale. 

A second analysis of the overlap-contribution has been performed by Manohar and 
Stewart [70]. They propose a new type of factorization in SCETn which separates 
modes in their virtuality k"^ and in their rapidity k_. Their method is based on a 
subtraction of zero-bin modes which is claimed to avoid double counting of soft and 
coUinear modes in the effective theory. As a result a divergent convolution integral 



where /i_ appears as the factorization scale in rapidity space. In [79] Manohar 
and Stewart present a factorization formula for the overlap- contribution ^(g^) which 
involves leading and sub-leading twist as well as three-particle distribution ampli- 
tudes together with their derivatives at the endpoints (without the need to introduce 
soft-collinear messenger modes). 

A still unresolved issue in this approach concerns the cancellation of the depen- 
dence which requires a more careful analysis of higher-order perturbative corrections 
beyond fixed-order perturbation theory. It is also to be emphasized that the quan- 
tities 0^(0), which naturally appear in this framework due to the 0-distributions in 
(11.311) . have to be interpreted as independent hadronic parameters since they cannot 
be derived from a finite number of moments of the distribution amplitude (pniu), 
see the argument in [SH]- Unfortunately, the factorization formula loses some of its 
power in this way since it involves a larger number of new hadronic parameters. 
However, the approach from Manohar and Stewart has appeared very recently and 
deserves further investigations. 

In our analysis in Chapter [5] we start from a different viewpoint. We consider heavy- 
to-light form factors between non-relativistic bound states which can be addressed in 




(1.30) 



is replaced by 
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fixed-order perturbation tlieory. As a consequence, all quantities in the factorization 
formula (11.281) . in particular the overlap-contribution ^(g^), can be calculated explic- 
itly in our set-up. Our results shed further light on the physics of the soft-overlap 
contribution and also allow to address the proposals mentioned above. 



Chapter 2 

Perturbative corrections 



In this chapter we address the technical aspects of the calculations that we discuss in 
the second part of this thesis. Most of the techniques that we present in the following 
have been developed within the last ten years, reflecting the steadily increasing 
effort in performing high precision calculations. We do not intent to give a review 
on multi-loop techniques but rather give a presentation from the perspective of our 
calculations in Chapter [HI H] and E] in the hope that the reader can get an idea what 
these calculations looked like in detail. 



2.1 Strategy 

We start with an outline of the general strategy that we have used to tackle the 1- and 
2- loop calculations in Chapter [3l H] and [5l The calculations have been performed on 
the basis of the computer algebra system Mathematica. Unless otherwise stated, 
all routines have been written by ourselves. We emphasize that our algorithm is not 
restricted to calculations in the Standard Model. However, for calculations that are 
even more complex than our 2-loop calculation in Chapter [3] and HI the efficiency of 
our algorithm should be improved. Our strategy consists of the following four steps: 

Step 1: Set-up for loop calculation 

We deal with up to 80 Feynman diagrams in our calculations. We first examine the 
kinematics of the process and calculate the colour factor of each diagram. In a 1-loop 
(2-loop) calculation we denote the loop momentum by k {k,l). The diagrams are 
expressed with the help of a minimal set of denominators Vi of propagators which 
we specify in each of our calculations in the second part of this thesis. 

We use a general tensor decomposition to express all tensor integrals as linear com- 
binations of scalar integrals Si multiplied by linearly independent tensors that can 
be formed out of the external momenta pi and the metric tensor g. For example in 
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the case of two external momenta pi and p2, a 1-loop integral of rank r = 2 becomes 

d'k — ^ rn'gf^'^ s^+p'tp^,S, + p^,p^,Ss + {ptpl^ + p^y,)S,, (2.1) 

/ 1 ■ ■ ■ I p 

reflecting the symmetry in of the tensor integral which reduces the number 

of independent tensor structures from five to four in this example. In contrast to 
this a 2-loop tensor integral of rank (r^, r^) = (1, 1) is decomposed into 

d'k dH — ^ = m'g^'' S[ + p^p^ S', + p^p^ S'^ + p^p^ S', + p\p^^ S'„ 

(2.2) 

since the tensor integral has no manifest symmetry in /i ^ z/ in this case. As the 
tensor decomposition becomes less trivial for higher ranks of the loop momenta, we 
develop a systematic decomposition in Section 12.21 

On the other hand we will see below that the reduction algorithm requires scalar 
projections Pi of the tensor integrals as an input. In our first example from (12.11) . 
we can form four different projections given by 

J J I ... I p 

whereas in the second example (12.21) . there are five projections P/ which read 

rp' p' jD' p' pn - f ^d,,dj {(fcpi)(/pi), {kpi){lp2), (fcp2)(/pi), (fcp2)(/p2), jkl)} 
x^ii ^z-> = a I ^j^r^ . 

(2.4) 

The scalar integrals Si and the projections Pi are related by a linear system of 
equations which can be solved easily. This allows us to express each tensor integral 
in terms of projections Pj which will be the subject of the second step. 

Step 2: Reduction to Master Integrals 

Our calculations are performed with the help of an automatized reduction algorithm 
which allows us to express the projections Pj of the general form 

d'^kd'^l \ " (2.5) 

in terms of a small set of so-called Master Integrals (Mis). The Si in (12. 5p denote 
scalar products of a loop momentum with an external momentum or of two loop 
momenta as in the examples in (12.31) and (12. 4p . The reduction algorithm will be 
described in detail in Section [^751 We emphasize that it makes the use of Dimensional 
Regularization (DR) mandatory for the calculation as anticipated by our notation 
writing d-dimensional integration measures. 
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Our interest to implement a reduction algorithm is threefold: First, multi-loop cal- 
culations as the one that we present in Chapter [3] and H] are extremely complex 
including typically several thousands of integrals making an automatization indis- 
pensable. Second, even in less complicated 1-loop calculations as the one in Chap- 
ter the algorithm is very helpful to avoid errors in the calculation as it is basically 
reduced to the computation of some Mis. Finally, the algorithm provides a powerful 
tool that can be used in the future for any kind of loop calculation as it is not based 
on certain Feynman rules but on particle propagators. 

Step 3: Manipulation of Dirac structures 

With the loop integrations reduced to a minimal form, we consider the Dirac struc- 
tures of the diagrams in the third step. We use the programme Tracer [86] for all 
manipulations concerning the Dirac algebra. We express all diagrams in terms of a 
minimal set of irreducible Dirac structures Di. A generic diagram "Dj takes the form 



where we indicated that the coefficients Aijk depend on the dimension d. We do not 
go into the details concerning the reduction of the Dirac structures here as it depends 
on the considered calculation and represents in general only a minor problem. 

The treatment of 75 in 7^ 4 dimensions requires special care. The matrix 75 
enters our QCD calculations because of the weak vertices and/or the projection 
on pseudoscalar meson states. We apply the Naive Dimensional Regularization 
(NDR) scheme which treats 75 as a completely anticommutating object. Despite 
the fact that this scheme is algebraically inconsistent [871489] . it leads to correct 
results provided that we can avoid traces as Tr(7^7'^7''7'^75) [90] . 

Step 4: Calculation of Master Integrals 

The most difficult part finally consists in the calculation of the Mis. We have devoted 
Section 12.41 to this subject where we present several advanced techniques that we 
have found useful in our calculation. We are looking for a solution of the Mis in 
form of an expansion in e = [4 — d)/2 



Expanding fl2.6p then determines the maximal order in the expansion of the MI that 
is required for the calculation. In our calculations, the expansion starts at most 
with double (quartic) poles due to soft and coUinear IR singularities at the 1-loop 
(2-loop) level. 




(2.6) 




(2.7) 
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2.2 Decomposition of Tensor Integrals 

The tensor decomposition has aheady been illustrated by means of a simple example 
in the last section. The procedure for more complicated tensor integrals is straight- 
forward but starts to become involved for higher ranks {rk,ri) in the 2-loop case. 
We therefore find it useful to present a systematic algorithm for the decomposition 
of 1- and 2-loop tensor integrals which can easily be extended to the multi-loop case. 



1-Ioop integrals 

We consider the decomposition of a 1-loop tensor integral of rank r 

d^k -V-^^ (2.8) 

/ 1 ■ ■ ■ I p 

according to n external momenta pi,...,pn- Notice that the integral is totally 
symmetric under the exchange of any pair of indices fii ^ fij. 

We start with all tensors that can be formed out of r external momenta pi. In total 
there are n'^ tensors of this type with ("^^~ ^) totally symmetric combinations. E.g. 
for n = 2 and r = 3 this results in (3) =4 totally symmetric tensors which read 

pfpf + + Prpfpf , Pfpf • (2.9) 

If r > 2, we continue with all tensors that can be formed out of r — 2 external 
momenta Pi and one insertion of the metric tensor g. From the external momenta, 
we obtain n^~'^ combinations. Further, there are {r-2) possible insertions of the 
metric tensor giving 1 2)72''"^ tensors of this type. From these, ("^-2^) are totally 
symmetric. In our example with n = 2 and r = 3 we thus find (1) =2 totally 
symmetric tensors given by 

^Mi^M2^M3 ^ g^'^^"'p'^' + ^^^^^pf + g^^^'p'^\ 

^MlM2pM3 ^ ^M3Mlp^2 ^ ^A'2M3pW _ 10) 

The procedure continues if r > 4 with r — 4 external momenta Pi and two insertions 
of the metric tensor g. We then find ("^^4 ^) totally symmetric combinations etc. 

We conclude that the tensor integral can be decomposed in this way according to 

tm^j: ("^::;-'^) (2.11) 

totally symmetric tensors where [x] denotes the greatest integer less than or equal to 
X. The proof follows by complete induction, but we refrain from presenting it here. 
We refer to Table 12.11 for examples that are relevant in typical 1-loop calculations. 
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n 

U 


1 
1 


9 


Q 
O 




1 


1 


1 


2 


2 


3 


2 


1 


2 


4 


6 


9 


3 


1 


3 


7 


13 


22 



Table 2.1: Number T(n;r) of totally symmetric tensors of rank r that can 
be formed out of n external momenta and the metric tensor ( relevant for the 
decomposition of 1-loop tensor integrals). The cases T{2;2) and T(2;3) have 
been discussed explicitly in Section [2A\ and \2.2\. respectively. 



2-loop integrals 

We come to the decomposition of a 2-loop tensor integral of rank (r^, r^) 

d'^kdH ^ 2.12 

according to n external momenta In this case, the integral is totally 

symmetric under the exchange of any pair of indices /ij ^ fij and z/j ^ uj. 

We start with all products that can be formed out of a rfc-dimensional tensor with 
indices /ii, . . . , /x,.^ and r^-dimensional tensor with ui, . . . ,1/^^. From our analysis 
of 1-loop tensor integrals we conclude that we can form T[n]rk)T{n]ri) tensors 
in this way which have the desired symmetry properties. To illustrate what we 
have obtained so far, we reconsider the case n = 2 and {rk,ri) = (1,1) which we 
introduced in Section 12. 1[ In this case, there are five independent tensors given by 



k'^l" - Ptpl. P^K, KK, pX, (2.13) 

whereas T(2; 1) T(2; 1) = 4 which corresponds to the first four tensors, only. We see 
that we have to add all contributions with metric tensors connecting the /i^ with the 
Vj. There are T(n; — 1) T(n; — 1) combinations with one metric tensor respecting 
the symmetry constraint, T(n; — 2) T{n] ri — 2) with two metric tensors etc. 

We conclude that the tensor integral in (12.121) can be decomposed according to 

mm(rfc,r;) 

r(n;rfc,n)= T{n;rk - j)T{n;ri - j), (2.14) 

i=o 



tensors respecting the symmetry in /ij ^ jjij and z/j ^ Uj. We again refrain from 
presenting a proof and give an explicit example for illustration. Let us consider the 
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case n = 2 and {rk,ri) = (2, 1). We find 10 independent tensors given by 

Pi P2 P2 + P2 Pi P2 P2 Pl> P2 ^2 P2 



^m/..^^ (2.15) 

according to T(2; 2, 1) = 4 ■ 2 + 2 ■ 1 = 10. We refer to Table [2^2] for further examples 
that are relevant in typical 2-loop calculations. 



rk\ri 





1 


2 


3 


4 





1 


1 


2 


2 


3 


1 


1 


2 


3 


4 


5 


2 


2 


3 


6 


7 


10 



rk\ri 





1 


2 


3 


4 





1 


2 


4 


6 


9 


1 


2 


5 


10 


16 


24 


2 


4 


10 


21 


34 


52 



Table 2.2: Number T{n;rk,ri) of tensors that can be formed out of n external 
momenta and the metric tensor which are totally symmetric in the first r^ 
and in the last ri indices. The tables correspond to n = 1 and n = 2 and are 
relevant for the decomposition of 2-loop tensor integrals. The cases T{2; 1, 1) 
andT{2; 2, 1) have been discussed explicitly in Section\2A\and\KM respectively. 
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2.3 Reduction to Master Integrals 



The reduction algorithm is a very important element in our strategy for the cal- 
culation of higher order perturbative corrections. In principle it can be used for 
the reduction of any (scalar) loop integral but it reveals its full power only in the 
multi-loop case. To give an illustration, it served in our 2-loop calculation from 
Chapter [3] and H] to express about 6.000 integrals in terms of 36 Mis which are in 
addition much simpler than the original integrals. Needless to say that this 2-loop 
calculation would have been (almost) impossible without the implementation of the 
reduction algorithm. 

The starting point is the representation of a scalar loop integral in the form 



where the Si are scalar products of a loop momentum with an external momentum 
or of two loop momenta. The Vi denote the denominators of propagators and the 
exponents fulfil ni,mi > 0. We focus on 2-loop integrals in the following with obvious 
simplifications (extensions) in the 1-loop (multi-loop) case. Notice that an integral 
can have different representations in terms of {S, V, n, m} because of the freedom to 
shift loop momenta in DR. In fact the underlying topology, i.e. the interconnection 
of propagators and external momenta, uniquely defines the integral. In the following 
we (loosely) use the word topology in order to classify the integrals. An integral 
with t different propagators Vi with nj > is called a t-topology. 

Preliminaries 

We first classify the scalar products Si into reducible and irreducible scalar products. 
The reducible scalar products can be written as linear combinations of denomina- 
tors Vi which leads to the first simplification. This is sometimes called Passarino- 
Veltman Reduction [91] which represents a standard tool for the computation of 
loop integrals. Since it appears as a preliminary step in our reduction algorithm, we 
briefly illustrate the idea by means of a simple example. 

We consider the following 5-topology 



where the right-hand side shows the underlying topology of the integral. We use 
dashed internal lines for massless propagators. The solid external line corresponds 
to the virtuality of the incoming momentum p in this case. As the figure on the 
right-hand side reflects the denominator of the integral only, the scalar products in 
the numerator appear explicitly in our notation. 






kp — ,' 



(2.17) 
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The scalar product turns out to be reducible in this 



kp 





kp = - {k + pY — k"^ — p^ 




(2.1 


ives rise to 










1 /'7^\ 1 /"'^\ 

2 2 


2 





(2.19) 



Notice that the first two integrals on the right-hand side correspond to simpler 4- 
topologies and that the last one has a trivial numerator. Moreover, the advantage 
of our notation becomes obvious in fl2.19p as it allows us to recognize easily that the 
first two integrals are topologically equivalent. This brings us to 



kp 



p" 



(2.20) 



All reducible scalar products can be rewritten along these lines. We may therefore 
concentrate on integrals of the form fl2.16p with irreducible scalar products in the 
numerator. In the 1-loop case there are n + 1 scalar products Si which can be 
formed out of n external momenta and the loop momentum: {fcpi, . . . , kpn-, k'^}- For 
a 1-loop t-topology we thus find n+l — t irreducible scalar products. In the 2-loop 
case the irreducible scalar products amount to 2r?, + 3 — t. This implies that in our 
example discussed above (n = 1, t = 5) all scalar products Si are in fact reducible. 



Integration-by-parts identities 

The core of the reduction algorithm consists in the use of integration-by-parts (IBP) 
identities |92]|93] which follow from the fact that surface terms vanish in DR 



(fk (fl 



d ...SI 



0, 



V e {k,l}. 



(2.21) 



In order to obtain scalar identities one may contract 02.211) with any loop or external 
momentum under the integral before performing the derivative. In the 1-loop (2- 
loop) case this generates n + 1 {2n + 4) identities from each integral. 

We illustrate the use of IBP identities by means of the famous example from 
Chetyrkin and Tkachov [92f93] . Let us consider the integral 



d'^kd'^l 



1 



k^{k + py{k-iY{i + pYP 

which can be reduced with the help of the following IBP identity 

d 1 



d^kdH 



^k^' 



{k-iy 



k^{k + pf{k - + p)H^ 



0. 



(2.22) 



(2.23) 
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We perform the derivatives explicitly on the integrand which leads to 



2k -{k- I) 



-2{k + p) -{k-l) 



-2{k-lf 



(2.24) 



where the dotted propagators are taken to be squared. The scalar products can be 
reduced with the help of the Passarino-Veltman Reduction to give 



-2k -{k-l) 



-2{k + p) -{k-l) 



~2{k - If : = -2 

Combining f l2.24p and ( I2.25P yields the famous result 



(2.25) 



d-A 



(2.26) 



The original 5-topology has thus been reduced to two 4-topologies, the first of them 
being particulary simple as it represents a product of two 1-loop integrals. 



Lorentz-invariance identities 



It has been pointed out in [91] that similar Lorentz-invariance (LI) identities can 
be formulated which become important in problems with many external momenta 
(multi-leg integrals). The LI identities exploit the fact that the integrals in 02.161) 
transform as scalars under a Lorentz-transformation of the external momenta. In 
this way one may generate up to six identities from each integral depending on the 
number of external momenta. 

As the LI identities are of minor importance in our calculations, we refer to [9l] for 
a detailed description. In our applications in Chapter [3l H] and [5l we deal with two 
external momenta p, q which leads to one LI identity per integral given by 



d'^k dH 



p ■ q \ p' 



dp^' 



^q^' 



+ q^p^ 



_d_ 

dq'' 



_d_ 



0. 

^2.27) 
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Reduction algorithm 

As we have seen in our exphcit example, the identities relate various integrals of 
the form fl2.16p with different exponents {?T,,m}. In general each of the identities 
contains the so-called seed integral which has been used to generate the identity, 
simpler integrals with smaller {n, m} and more complicated integrals with larger 
{n,m}. The identity is then used to reduce the most complicated integral rather 
than the seed integral itself. 

We do not want to inspect single identities but are looking for an automatized 
reduction algorithm which exploits the full power of IBP and LI identities. This 
can be achieved in a bottom-up approach generating systematically all identities 
from all sub-topologies of a given integral. In the 2-loop case, we thus start with 
the identities from all 2-topologies, 3-topologies etc. which leads to a large set of 
identities which contains the integral that we are interested in beneath simpler and 
more complicated ones. If we were to use all unknown integrals that appear in 
these identities again as seed integrals, the procedure would be without end as we 
generate identities with more and more complicated integrals in each step. It was an 
important observation from Laporta that the number of identities grows faster than 
the number of unknown integrals in this procedure. At some point we may therefore 
stop the generation of new identities and solve the (apparently) over-constrained 
system of equations by expressing more complicated integrals in terms of simpler 
ones. Not all of the identities being linearly independent, some integrals finally turn 
out to be irreducible to which we refer as Mis. 

The choice of the point where to stop the outlined procedure requires thorough 
experimentation. On the one hand we want the system of equations to be large 
enough to be sure that the remnant integrals in the reduction procedure are indeed 
irreducible. On the other hand the system should be as small as possible to assure 
an optimized realization of the reduction algorithm. Based on his own experiences 
Laporta proposed a "golden rule" for the choice of such a cutoff in [95] which we 
have found very useful. 

In our calculation we typically deal with systems of equations made of several thou- 
sands equations. The solution being straight-forward, the runtime of the reduction 
algorithm depends strongly on the order in which the equations are solved. As a 
guideline for an efficient implementation we have followed the algorithm described 
in [Sg. 
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2.4 Calculation of Master Integrals 

We now present a collection of techniques that we have used for the calculation of 
the Mis. This last step in our strategy from Section [2?T] represents the most difficult 
part of the perturbative calculation. We always intended to compute all Mis with 
two independent methods in order to become sure of the correctness of our results. 

We illustrate the calculation techniques with two explicit examples: The first one 
corresponds to a 1-loop integral that appears in our calculation from Chapter [5l 
the second one to a 2-loop integral from Chapter [3] and HI Albeit this section has 
become somewhat lengthy in this way, we think that it is instructive to present the 
calculation of these integrals with different methods. In both cases we have chosen 
simple (but non-trivial) examples which allow for a transparent presentation of the 
central aspects of the techniques. Despite being simple, the examples cover most of 
the conceptually interesting features. 




1-Ioop example: Definition 

In the first example we consider a 3-topology 

^ ^ [{mw' — mw + ky — m'^][k'^ — 'm'^]{k — mw')'^ " ''" ^ 

with w"^ = w''^ = 1 and w ■ w' = > 1. In our analysis in Chapter [5l we focus on the 
leading power in an expansion in m/M where M is the mass of a heavy quark and 
m the mass of a light (non-relativistic) quark. The parameter 7 corresponds to the 
large boost between the rest frames of the bound states and scales as 7 = 0{M/m). 
In f l2.28p we used dashed/solid internal lines for massless/massive propagators. 
Dashed/solid external lines denote virtualities of 0{m^) / 0{mM) which are writ- 
ten explicitly here as multiples of m?. The normalization of the integral reads 



[dk] ^ d'^k. (2.29) 



2-Ioop example: Definition 

In the 2-loop case we consider a massless 4-topology ^ 

f'^^^f'^^^ {p-q- lY{uq + -q-k-l)\uq + kf = -r-----^ (2.30) 

with 0<u<l, «=! — M, p^ = 2p-g = m? and q^ = 0. The usual ie-prescription of 
propagators is understood. Internal dashed lines correspond to massless propagators 
and dashed/sohd/double external lines to virtualities O/wm^/m^, respectivelju. 



"'^Notice that m denotes the mass of the 6-quark in this case, whereas the hght quarks are 
considered massless. 
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2.4.1 Feynman Parameters 

The standard method to compute loop integrals introduces Feynman parameters. 
As this approach is not practicable for most of our Mis, we do not describe this 
method here in detail. However, it turns out that our simple examples can be 
solved with this technique. We comment briefly on the calculations and quote their 
results which will serve as a reference for the other calculations. 

1-Ioop example: Feynman Parameters 

As we deal with a 3-topology, we introduce two Feynman parameters and obtain 

-{m'^)-^-' r(i -£)r(i + £) 

[(5 - 47)a;2 + y"^ + 2(2 - -f)xy + 4(7 - l)x]^^' ' 

(2.31) 

The integral is ultraviolet (UV) and infrared (IR) finite. We focus on the leading 
term in the e-expansion which allows us to set e = 0. The integral can then be 
solved in a closed form in m/M, but we refrain from presenting the exact result 
since it looks rather complicated. We are only interested in the leading power of the 
m/M expansion which takes a simple form 





IT 



21n21n7 + ln^2 + — + 0(e) 
3 



^2.32) 



Notice that we cannot expand the integrand in 02.311) in m/M before performing 
the parameter integrals as this would invoke IR singularities for x — 0. In other 
words, the two expansions in e and m/M do not commute with each other which 
complicates the extraction of the leading power. 



2-Ioop example: Feynman Parameters 

In this example, we introduce three Feynman parameters but one of these integra- 
tions can be done trivially. We obtain 



-2, „2,„ r(2E) r(i - B) 



1 nl-x 

dx I dy x^~'^{l - X - yy'^^x + uyy'^\ 







The integral exhibits a double pole in e reflected by the factor r(2£:) and a second 
singularity from x —>■ 0. The solution can be found in a closed form in e in terms of 
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an hypergeometric function and reads 




^2,., 2r(-25)r(i-e)^ 

r(2 - 2e) r(2 - 3£) 
r(2£) 2^1(1, 2e; l + e-u)- eT{ef u'^u'' 



(2.34) 



The integral enters our calculations from Chapter [3] and H] up to O^e^), but we 
restrict our attention to the first three coefficients here. Notice that the calculation 
of the imaginary part essentially requires one coefficient less in the e-expansion. 
The hypergeometric function can be expanded with the help of the Mathematica 
package HypExp [HS] which yields 



2.4.2 Method of Differential Equations 

The method of differential equations [251^ in combination with the formalism of 
harmonic polylogarithms (HPLs) [30] turned out to be extremely useful for our 
2-loop calculation from Chapter [3] and HI The method is based on the reduction 
algorithm that we presented in Section 12.31 and can therefore only be applied with 
such an algorithm at hand. 

In the following we start with a general description of the method, introduce the 
HPLs and comment briefly on the calculation of the boundary conditions to the 
differential equations. Finally, we apply this method to our explicit 2-loop example. 
We do not use this technique in our 1-loop calculation from Chapter [5] where we 
focus on the extraction of the leading power in m/M. We therefore refrain from 
calculating our 1-loop example with this method. 

Solving loop-integrals with differential equations 

The Mis are functions of the physical scales of the process which are given by scalar 
products of the external momenta and masses of the particles. In our calculations 
most of our Mis depend on two scales which give rise to one dimensionless ratio: 7 
(m) in our 1-loop (2-loop) example. We therefore restrict our attention to the case 
that there are only two distinct scales in the process with obvious modifications for 
the general case. We denote the ratio of the two scales by u and recall the general 




(2.35) 



+ Li2(u)H — In^ u + Inwlnu — 51nu + 0(e) . 



form of a MI 




(2.36) 
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We perform the derivative with respect to u and interchange the order of integration 
and derivation 



The right-hand side being of the same type as the IBP and LI identities in fl2.2ip 
and fl2.27p . this procedure again leads to a sum of various integrals with different 
exponents {n, m}. With the help of the reduction algorithm, these integrals can be 
expressed in terms of Mis which yields a differential equation of the form 



where we indicated that the coefficients a and bj depend on the dimension d. The 
inhomogeneity of the differential equation typically contains Mis of sub-topologies 
which are supposed to be known in a bottom-up approach. In some cases few Mis 
in the inhomogeneous part are of the same topology as the MI on the left hand 
side of f l2.38p and thus unknown. Writing down the differential equations for these 
Mis, we find that we are left with a coupled system of linear, first order differential 
equations. 

We are looking for a solution of the differential equation in form of an expansion 



Expanding (12.380 then gives much simpler differential equations for the coefficients 
Cij which can be solved order by order in e. The solution of the homogeneous 
equations is in general straight-forward. The inhomogeneous equations can then 
be addressed with the method of the variation of the constant. This in turn leads 
to indefinite integrals over the inhomogeneities which typically contain products of 
rational functions with logarithms or related functions as dilogarithms. With the 
help of the formalism of HPLs these integrations simplify substantially. 

Harmonic Polylogarithms 

The HPLs have been introduced in ^30j and several extensions of the formalism have 
been considered in pTH^ . We briefly summarize their basic features here, focussing 
on the properties that are relevant for our 2-loop calculation in Chapter [3] and HI 

The HPLs, denoted by H{ described by a w-dimensional vector of 

parameters and by its argument x. In their simplest form, the parameters can take 
the values and ±L The basic definitions of the HPLs are for weight w = 1 




(2.37) 




(2.38) 




(2.39) 



H{0;x) 
H{l-x) 

H{-l;x) 



\nx, 



— ln(l — x), 
ln(l + x) 



(2.40) 
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and for weight w > 1 

H{a,mu)-i; x) = / dx' f{a;x') H{m^_i;x'), (2-41) 
Jo 

where the basic functions /(a; x) are given by 

f(l;x)^^Hil;x) = -^, 
ax 1 — X 

f{-l;x)^-^H{-l;x) = -^. (2.42) 
ax 1 + X 

In the case rhy^ = 0^, the definition in fl2.4ip does not apply and the HPLs read 

H(0,...,0]x) = ^liTx. (2.43) 
wl 

The HPLs form a closed, linearly independent set under integrations over the basic 
functions /(a; x) and fulfil an algebra such that a product of two HPLs of weight 
Wl and W2 gives a linear combination of HPLs of weight w = Wi + W2- 

As anticipated above, the solution of the differential equation typically leads to inte- 
grals over products of rational functions with transcendental functions as logarithms 
or dilogarithms. More precisely, we often encounter integrals of the type 



dx' (1 H{m^;x') (2.44) 

{x' 1 — x' 1 + x' ) 

which become trivial within the formalism of HPLs as they simply correspond to an 
HPL of weight w + 1 according to fl2.4ip . Further integrals take e.g. the form 



(2.45) 



The solution of these integrals is also straight-forward within the formalism of HPLs 
as an integration-by-parts relates them either to an HPL of weight w + 1 or gives rise 
to a simple recurrence relation. Not surprisingly, the pattern in fl2.44p and (12.451) 
does not apply to all integrals that we encounter in our 2-loop calculation from 
Chapter [3] and m However, a large part of this lengthy calculation can be performed 
along these lines. 

Furthermore, the formalism of HPLs helps us to avoid the problem of "hidden zeros" . 
Let us illustrate this point with a simple example. In our calculation we may obtain 
the result 



7r2 



Li2(M) + Li2(u) + Inulnw — 

6 



- (2.46) 

e 
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and be surprised that it contains a singularity. A closer look reveals that it is indeed 
free of any singularities since (12.461) vanishes identically. This can be seen easily 
when we rewrite the second dilogarithm in the form 

Lis (u) = In M In M - Lia (u) . (2.47) 

6 

We therefore call (12.461) a "hidden zero" . In complicated multi-loop calculations with 
HPLs of higher weight and various arguments, it might become difficult to identify 
"hidden zeros". In order to avoid this problem, it is important to use a formulation 
in terms of a minimal set of transcendental functions. The HPLs provide such a 
minimal set. To give an example, in our 2-loop calculation from Chapter [3] we were 
able to express our results in terms of the following HPLs 

H{0;u) = Inu, 

H{l;u) = -ln(l -m), 

H{0,0,l;u) = U3{u), 

/f(0,l,l;M) = Si,2(n). (2.48) 
Comment on boundary conditions 

A unique solution of a differential equation requires the knowledge of its boundary 
conditions. In the considered case the boundary conditions typically represent single- 
scale integrals corresponding to u = 0, 1. It is of crucial importance that the integral 
has a smooth limit at the chosen point such that setting u = or u = 1 does not 
modify the divergence structure introduced in (12.391) . 

The boundary conditions can be addressed with any method presented in this sec- 
tion. In some cases the calculation becomes trivial as they correspond to simpler 
topologies which may turn out to be reducible. If so, the integral can be expressed 
in terms of known Mis with the help of the reduction algorithm. For the more com- 
plicated integrals this is not the case and the calculation of the boundary conditions 
represents sometimes the most difficult part in the computation of the Mis. 

2-Ioop example: Method of Differential Equations 

With the help of the reduction algorithm we derive the following differential equation 

— ''"'"^^^-^ - ^ (d - 4) ''"^"^^^ + (2 49) 

du 2uu 2uum? \^ /' 

where the inhomogeneous part contains a 3-topology which is supposed to be known 

,M-2. ...sr(-l + 2e)r(l-e)5 



r(3 - Ze) 



(2.50) 
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Concerning the boundary conditions, we examine if the integral has a smooth limit 
for M 0,1. From (12.331) we see that setting u = modifies the divergence structure 
of the integral in the form — > x""^"^, whereas u = 1 does not. We therefore use 
the latter boundary condition which corresponds to a reducible topology 



3^ 8 4:'--:^ = fm2)-2-e2-^ 



^ ^ 19 

(c/-4)m2 V y ' ^ [2^2 + 2£ ^ 2 ^ 

(2.51) 

We solve the differential equation (I2.49P order by order in e with the ansatz 



^2y2e ^2i.e , X ^ ^^ ^ ^^^^ I _ ^2.52) 







Notice that the differential equation has a particular simple structure in our example 
because of the explicit factor [d — 4) = 0{6). However, the non-trivial element in 
the solution of the differential equation consists in the integration over the inhomo- 
geneity which can be illustrated here. In the lowest order, we simply find 

d 

— C-Ju) = — *• c_4(u) = const (2.53) 
ou 

and the comparison with (12.511) leads to c_4(-u) = 0. The computation of c_3 follows 
exactly the same lines. The first non-vanishing coefficient is c_2 which again obeys 

d 

— c_2(^i) = c_2(m) = const, (2-54) 
ou 

but the boundary condition (12.511) yields c^2{u) = |. From now on, the differential 
equations take a non-trivial form which perfectly fit into the formalism of HPLs 

d 1 

— c_i(m) = c_i(m) = -H(0; u) + const. (2.55) 

ou u 

The boundary condition results in c_i(m) = | — i/(0; u). In the next step, we have 

^c,{u)={---\ H{0-u)-- 
Ou \u u J u 

Co{u) = H{0,0;u)-H{l,0;u)~5H{0;u) + const. (2.56) 

The comparison with (ICTj) gives co{u) = f - f + H{0,0;u) - H{l,0;u) -5H{0;u) 
which can be expressed in terms of the minimal set (12.481) with the help of 

HiO,0;u) = ^HiO;u)', 

H{1, 0; u) = H{0; u)H{l; u) - H{0, 1; u). (2.57) 

Collecting all coefficients Cj and using (I2.52p with (I2.48p . we finally reproduce our 
result from (12.350 . It is evident that this MI has been particularly simply but it 
nevertheless illustrated the usefulness of the HPLs in this context. 
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2.4.3 Expansion by Momentum Regions 

The method of expansion by regions [31j is one of the most powerful techniques to 
determine loop integrals to leading power in some small expansion parameter. We 
have already pointed out that the extraction of the leading power would simplify 
substantially if we could interchange the order of loop integration and power expan- 
sion. However, we have seen in the analysis of our 1-loop example after (12.321) that 
this procedure may generate artificial IR singularities which invalidates this naive 
approach. 

We will see in this section that we can do so if we split the loop integration into 
individual momentum regions and expand the integrands in each region separately. 
This in turn generates artificial singularities in each region which cancel when we 
reconstruct the full integral as the sum over all regions. When the calculation is 
performed within DR, no further cutoffs are required to separate the momentum 
regions due to the fact that scaleless integrals vanish in DR. We now switch to our 
explicit 1-loop example to demonstrate how this technique works in detail. 



1-loop example: Expansion by Momentum Regions 

We first address a power-counting to the external momenta and particle masses 
in the problem. We adopt the SCET terminology from Table 11.11 in our example 
by identifying = m/M in the non-relativistic set-up. The momentum of the 
spectator antiquark in the initial bound state corresponds to a soft momentum with 
mw = (m, 0,m) ~ (A^,0,A^) whereas the quark and antiquark in the final bound 
state are in a collinear configuration with mw' = (27m, 0, m/27) ~ (1, 0, A^). Notice 
that this parametrization fulfils the conditions = w'"^ = l,iy-w'~7~l/A^ 
(transverse components are suppressed in the non-relativistic approach). 

We want to compute the leading power in A of the integral 

^ ^ [{mw' — mw + kY — m'^][k'^ — m'^]{k — mw'Y ^ ^ 
Potential leading contributions come from hard, hard-collinear, soft and collinear 
momentum configurations of the particles in the loop. We may also look for con- 
tributions from soft-collinear messenger modes as discussed in Section ll.3.2[ We 
perform a power counting of the integral in the individual momentum regions of the 
loop momentum k. This gives rise to the following pattern 




hard 


(1,1,1) 


0- 


-0- 


-0- 


-0 = 


= 0, 


hard-collinear 


(1,A,A2) 


4- 


- 2 - 


- 2 - 


- 2 = 


= -2 


soft 


(A^A^A2) 


8- 


- 2 - 


-4- 


- 2 = 


= 0, 


collinear 


(l,A^A^) 


8- 


- 2 - 


-4- 


- 4 = 


= -2 


soft-collinear 


(A^A^A^) 


12 - 


- 2 - 


-4- 


-4 = 


= 2. 
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The first term in each region comes from the scaling of the integration measure 
and the other terms from the three propagators. We find a leading contribution of 
(9(1/A^) which stems from the hard-collinear and the collinear momentum region. 

We first address the hard-collinear region. As anticipated above, we expand the inte- 
grand in this region counting k ~ (1, A, A^) and do not need to introduce additional 
cutoffs in DR. The hard-collinear contribution becomes 




r(i-^) 



[k^k_ + k\ — 27m A;+ + ie] 



[k^k_ + k\ + 27m(A;+ — m) — mk_ + ie] [A;+A;„ + k\ + ie] 



[2.m) 



where we have made the i5-prescription of the propagators explicit and translated 
the c?- dimensional integration measure from (12.291) into light-cone coordinates. We 
first perform the k_ integration by contour methods. As all poles lie in the upper 
(lower) half plane for < (A;+ > m), the only non- vanishing contribution comes 
from </;;+< m. In this case one of the poles lies in the upper half plane and the 
two others in the lower one which yields 




r(i-^) 



TT 



dk< 



k+ — m 



d'^-^kj 



[kl - A-fk+{m - k+)][kl - 2-fk+{m - k+)] ' 
The remnant integrations are straight-forward and give 



(2.61) 




r(l-e)r(l + e) 1-2- 



dk^ {2'~fkj, 



m — k^ 



— (27m 



2,_i„, 4- - T r(l - e) r(l + e) r(l/2) Ti-e) 



2 In 2 



e r(l/2-£) 

+ 2 In 2 ln(7m2) + 3 In^ 2 + 0{e) 



Now we consider the collinear region counting k ~ (1, A^, A^'' 



(2.62) 




r(i-e) 

2i7r'^/2 



dk_ dk+ d'^-^ki 



[kj^k^ -|- /c^ — m? + ie] 



(2.63) 



[—mk- — 27771^ -|- ie][kj^k^ + k^ — 2'~fmkj^ — mk^/{2'~f) + m? + ie] 
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We perform the /c+-integration with contour methods and find a contribution for 
< A;_ < 27m. The other integrations are again straight-forward. 




rn-e) 







Til-E)TiE) 
27m2 

rfi-£)r(£) 



27m 

dk. 



1 



d'^-^k 



27m2(A;_ + 27m) [kl - (m - A;_/(27))2] 



27771 



dk^ {k_ + 27m)"' (m - A;_/(27)) 



-2e 



(m") 



2^1(1, 1; 2- 26; -1) 



1 -2£ 



47777.2 



21n21n777^ - 21n^2 + — + 0{e) 

£ 3 



(2.64) 



Combining the two contributions from fl2.62p and (12.641) . we see that the artificial 
divergences that appeared in each region cancel each other. We finally reproduce 
our result from (I2.32p 




+ 




47771^ 



21n21n7 + ln^2 + — + 0{e 



. (2.65) 



Notice that f l2.62p and (12.640 determine the leading power of the integral in a closed 
form in e whereas we were only able to determine the first coefficient in the e- 
expansion with the help of Feynman parameters in Section 12.4.11 



2.4.4 Mellin-Barnes Techniques 

In the last two sections we have introduced two sophisticated techniques for the 
calculation of loop integrals. The first one was particularly suited for our 2-loop cal- 
culation from Chapter [H] and m the second one for the 1-loop calculation in Chapter^ 
We now present a third method which we used in both calculations. In our 2-loop 
calculation we typically applied Mellin-Barnes techniques [32)133] for the calculation 
of the boundary conditions, i.e. for single-scale integrals. As the boundary condition 
in our 2-loop example is rather trivial, we present the computation of the full 2-loop 
integral in the following. Mellin-Barnes techniques are also a comfortable tool for 
the calculation of loop integrals in a power expansion as we will see in our 1-loop 
example below. 

The method is based on the following representation 



1 



(X + Y)^ r(A) 2m 



+too 



dz V{\ + z)V{-z) r"X-"-^ (2.66) 



where the contour separates left poles from right poles (we specify the meaning of 
this phrases in our examples) and z is called a Mellin-Barnes parameter. We see that 
(I2.66P basically converts a sum into a product at the cost of a contour integration. 



2.4. CALCULATION OF MASTER INTEGRALS 



57 



We typically introduce Feynman parameter first and use the Mellin-Barnes represen- 
tation f l2.66p to simplify the corresponding integrations. With all integrations over 
Feynman parameters done, we examine the analytical properties of the integrand 
and perform the Mellin-Barnes integrations with contour methods. This leads to 
infinite sums over residues which often represent the main obstacle in this approach. 



2-loop example: Mellin-Barnes Techniques 

We first address our 2-loop integral from fl2.30p . The starting point is equation 
(12.331) which we rewrite with the help of the substitution x = st, y = st with the 
usual bar-notation understood. This yields 



^0 



r(2 - 2e) 

(2.67) 

We introduce one Mellin-Barnes parameter to split the last factor into a product. 
The integrations over the Feynman parameters give rise to F-functions and we obtain 

...C-^^ = (m2)-2-e2-^^ ^ (2.68) 

— [ dz T{2e + z)T{l + z)T{-£-z)T{-z)u'. 
2-ni Jc 

The integrand has simple poles which we classify into left poles which stem from 
r(. . . + z) and right poles from r(. . . — z). We find 

left poles: z = -2e, -1 - 2e, -2 - 2e, . . . z = -1, -2, -3, . . . 

right poles: z = — e, 1 — e, 2 — e, . . . 2; = 0, 1, 2, . . . (2.69) 

The contour C separates these two types of poles. Notice that we cannot simply set 
e = as this would glue left and right poles together at 2; = 0. There is no such 
problem for < e < 1; let us take 6 = 1/2 for the moment. We see that we can 
choose the contour C as a straight line in this case (cf. Figure [?!Th ). For e — the 
contour is deformed into C as shown in Figure 12.1b . Alternatively, we may stick 
to the contour C and add the contributions from the poles that crossed the straight 
line for e — explicitly. The crucial point is that all singularities of the integral 
stem from such contour crossings. In our example we pick up the pole for z = —2e 
and may safely set e = in the remnant contour integral which results in 



(m^) 



2\-2e „2i-Ke 



Til-e 



i4 



'l-2e)V{2-?>e) 

'3/4+joo 



T{l-2e)T{e)T{2e)u 



-2e 



1 /•— 3/4+joo / 

— / dz [V{z)V{l + z)V{-zf + 0{e) 



(2.70) 
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(a) 



(b) 



-3-2-1 1 2 



-3 -2 -10 12 



Figure 2.1: Possible integration contours in \2. 6^) that separate left poles 
from right poles. The figures show the position of the poles for e = 1/2 (a) 
and e = (h). All poles lie on the real axis, hut we have assigned a positive 
imaginary part for illustrative reasons. 



We close the contour to the left and sum the infinite residues from z = — 1,— 2,— 3, 

'1 + n\\iu) 



-1 3/4+joo 

— / dz V{z)T{l + z)T{-zf u"" 



E 

n=l 



(2.71) 



Li2('u) + In-u Inw In u 



7r 



2 3 

The expansion of (12.701) in e yields together with (I2.7ip our result from (12.351) . 



1-loop example: Mellin-Barnes Techniques 



In our 1-loop example we start from (I2.3ip and perform the substitution x = st, 
y = st as in the last example 




— m 



2\-l-e 



T(l-e)T(l + e) r 



ds dt 

Jo [t + 2s{2t + ts){-f -1)] 



l+e 



(2.72) 



We introduce two Mellin-Barnes parameters in this case to obtain standard integra- 
tions over Feynman parameters in terms of F-functions. We find 




1 

27Ti 



dw 2" 



„ r{w - z) F(-w) T{w -2e - z)T{l + z - w) F(l + w) 



T{2 + z + w) 



[2.n) 



The contours separate left poles from right poles if the arguments of the F-functions 
are positive. This can be achieved by integrating along straight lines with wq = —1/2 
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and Zo = — 3/4 for £ = (we already know that the integral is finite and set e = 
in the following). The poles in w can be classified according to 



left poles: w 
right poles: w 

We close the contour to the right and find 



z, -1 + z, -2 + z, 
l + z,2 + z,3 + z 



w 
w 



-1, -2, -3, 
0, 1 2, ... 



^2.74) 



— I dw 2 

2m 



_^ r{w - zf T{-w) T{l + z- w) r(l + w) 



E (-1)" 



n=0 



T{2 + z + w) 
V{n- zfV{l-n + z) 



+ 2- 



V{2 + n + z) 
Vi-l- 



-n—l—z 



n - z)T{2 + n + z)T{n + 1) 
r{3 + n + 2z) 



2(1 + z) r(-l - zf 2Fi{l, 2 + 2Z-2 + Z- -1) 
- 2 



-'-'^^^ "^^^"^ W + z)-^{?,/2 + z)] 



r(2 + 2z) 



(2.75) 



where il){z) is the digamma function given by ip{z) = r'(2;)/r(2;). So far our compu- 
tation has been exact in m/M. The important point to notice is that the remnant 
contour integral has poles for z &7a. If we close the contour to the left, the leading 
contribution in m/M stems from a single pole at z = —1 because of the explicit 
factor (7 - 1)^ ~ 7Mn fl2T3D . This yields 




TT 



2 In 2 ln7 + In^ 2 + — + 0{e) 



(2.76) 



in agreement with f l2.32p . We see that the Mellin-Barnes technique is particularly 
suited to compute the sub leading contributions in m/M which stem from subsequent 
poles in z in our example. 



2.4.5 Sector Decomposition 

Let us summarize what we have obtained so far. Concerning our 1-loop calculation 
from Chapter [HI we presented two analytical methods for the calculation of the Mis: 
the method of expansion by regions and Mellin-Barnes techniques. The comparison 
between these two independent calculations is extremely helpful to avoid errors in 
this part of the calculation. So far, we have found only one method for our 2-loop 
calculation from Chapter [H] and H) the method of differential equations in combina- 
tion with Mellin-Barnes techniques for the calculation of the boundary conditions. 
A second independent method would also be very helpful in this case. 
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We use a numerical method, the method of sector decomposition [3lj, in order to 
check our analytical results of the 2-loop integrals from Chapter [3] and HI We refrain 
from presenting this technique here in detail as we will only use it as an (impor- 
tant) check. The basic idea corresponds to a decomposition procedure of Feynman 
parameter integrals which systematically disentangles overlapping IR divergences. 
These divergences can then be isolated with the help of an adequate subtraction 
procedure. At this step the MI is given in form of a Laurent expansion in e with 
coefficients that correspond to regular parameter integrals that can be evaluated 
numerically. 

The power of the considered algorithm lies in the fact that it can easily be automa- 
tized and be applied, at least in principle, for any multi-loop integral. Its numerical 
precision is mainly limited by the potential presence of thresholds in the parameter 
integrals. The obtained precision therefore varies strongly in our applications be- 
tween 1 part in 10^-10^ depending on the considered integral. This situation could, 
however, be improved by performing the numerical integrations with more efficient 
routines. 



Part III 
Applications 



Chapter 3 

Hadronic two-body decays I: 
Imaginary part 

In the second part of this thesis we consider several perturbative calculations in 
the framework of QCD Factorization. We first compute the imaginary part of the 
NNLO vertex corrections in charmless hadronic two-body decays as e.g. i? — > tttt 
(the real part will be considered in Chapter Hj). The results presented in this chapter 
will be pubhshed in |100j . 

The outline of this chapter is as follows: We first introduce a new operator basis 
of the effective Hamiltonian which is particularly suited for multi-loop calculations. 
We then present some details of the 2-loop calculation following our strategy from 
Section 12. 1[ As the factorization formula reveals a rather complicated divergence 
structure at NNLO, we elaborate the subsequent UV and IR subtractions in some 
detail. The imaginary part of the topological tree amplitudes is finally obtained in 
an analytic form. We conclude with a brief analysis of the numerical impact of the 
NNLO vertex corrections. 



3.1 Change of operator basis 

In view of the calculation of topological tree amplitudes, we restrict our attention 
to the current- current operators (5i,2 of the effective weak Hamiltonian (11.81) . Due 
to the fact that we work within Dimensional Regularization (DR), we also have to 
consider evanescent operators [90]ll01|ll02j . These non-physical operators vanish in 
d = 4 dimensions but contribute at intermediate steps of the calculation in d ^ 4 
dimensions. We have emphasized in Section [l.2.3l that the imaginary part considered 
here has effectively NLO complexity. We will indeed see that the calculation of 
the imaginary part only requires 1-loop evanescent operators (2-loop evanescent 
operators will contribute to the real part which we consider in Chapter H]). For our 
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purposes the complete operator basis is thus given by 

Q2 = [uiYL bj] [dj'j^ L Ui] , 

E, = [uYYYLb] [J7^7,7^Lm] - (16 - As) Qi, 

E2 = [u^n^YL b,] [J,7m7.7p L u.] - (16 - 4^) Q2, (3.1) 

with L = 1 — 75. We refer to this basis as the (standard) QCDF basis for convenience 
and denote the corresponding Wilson coefficients and operators with a tilde. 

It has been argued by Chetyrkin, Misiak and Miinz (CMM) that one should use a 
different operator basis in order to perform multi-loop calculations P~03| . Although 
the deeper reason is related to the penguin operators which we do not consider here, 
we prefer to introduce the CMM basis in view of future extensions of our work. 
The CMM basis allows to consistently use DR with a naive anticommuting 75 to all 
orders in perturbation theory. In the CMM basis the current-current operators and 
corresponding 1-loop evanescent operators read (denoted with a hat) 

Qi = [uYLT^b] [d-f^LT^u] , 
Q2 = [uYLb] [d-ff,Lu] , 

E, = [uYYYLTH] [d^^^^^pLT^u] - I6Q1, 

E2 = [uYYYL b] [J7m7.7p Lu]-1QQ2. (3.2) 



Comparing fl3.ip and (13.21) we observe two differences: First, the two bases use 
different colour structures which is a rather trivial point. More importantly, they 
contain slightly different definitions of evanescent operators which we will examine 
now in detail. 

The issue is related to Fierz symmetry. We will see below that the definition of 
evanescent operators in the CMM basis explicitly breaks Fierz symmetry which 
relates the two diagrams in Figure 13. 1[ This can be seen by considering the UV part 



7 (5 7 5 








Q, 



Figure 3.1: Generic 1-loop diagram with different contractions of fields in a 
four-quark operator Qi. The two insertions are related by a Fierz reordering, 
see text (a,(3,j,5 are spinor indices). 
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of the left diagram which involves the following combination of Dirac matrices 

(3.3) 

which we reshuffled into a more convenient form using an anticommuting 75. On 
the other hand the right diagram gives 

(YYrLi.i,)^, {l,L)^, = id- 2f irL)^, il,L)^, id - 2f irL)^, il,L)^s , 

(3.4) 

where we have performed a Fierz reordering in the second step. If we now impose 
Fierz symmetry, the expressions (13.31) and (13. 4p have to be equal and we arrive at 

irYYL)^^ (7m7.7p^),, = (16 - 4e - 4e') {YL)^^ in.L)^, (3.5) 

If we look at the remnant 1-loop diagrams we find that (13. 5p is in fact the only 
constraint from Fierz symmetry. In NLO the terms of O(e^) do not contribute as 
the loop-integrals have at most 1/e (UV) divergences. They can therefore simply 
be neglected here (we stress that this will be different in our analysis in Chapter Hj). 
We conclude that the QCDF basis from (13. ip is Fierz symmetric whereas the CMM 
basis from (13. 2p is not. In other words, the freedom in the definition of evanescent 
operators has been used in the QCDF basis to properly adjust the e-terms into a 
Fierz symmetric form. 

Why do we care about Fierz symmetry? In the considered calculation the contrac- 
tion depicted in the left diagram from Figure 13.11 is related to the colour-allowed 
tree amplitude (ai), whereas the one from the right diagram leads to the colour- 
suppressed tree amplitude (02) • On the technical level these two possible insertions 
of a four-quark operator correspond to two completely different calculations. It 
would be very tedious if we had to perform both calculations explicitly, in particu- 
lar in the considered 2-loop case. We have pointed out in Section [1.2.31 that ai and 
a2 can naturally be related by Fierz symmetry. For this, it is of crucial importance 
that we preserve Fierz symmetry when we work in the effective theory which fac- 
torizes the amplitudes into Wilson coefficients and matrix elements within DR. As 
we have argued above, this is indeed the case in the QCDF basis which allows us to 
derive from a\ by simply interchanging C\ ^ C^- 

We conclude that the CMM basis is the appropriate choice for a 2-loop calcula- 
tion whereas the QCDF basis provides a short-cut for the derivation of the colour- 
suppressed amplitude. We therefore propose the following strategy for the calcula- 
tion of the NNLO vertex corrections: We perform the explicit 2-loop calculation in 
the CMM basis using the first type of insertion in the left diagram from Figure 13. 1[ 
In this way we obtain a\[Ci). We then transform this expression into the QCDF 
basis which yields a\[Ci) and finally apply Fierz symmetry arguments to derive 
OL^iCi) from a\{Ci) under the exchange C\ ^ C2. 
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3.2 2-loop calculation 

In this section we present a brief overview of the technical aspects of the considered 
2-loop calculation. Herein, we follow the systematics of our strategy from Section [2?T1 

Step 1: Set-up for loop calculation 

We will see explicitly in Section [3.3.21 that we may restrict our attention to (naively) 
non-factorizable diagrams, similar to what we have seen in the NLO analysis from 
Section 11.2.31 These diagrams contain at least one gluon which connects the two 
currents in the left diagram of Figure 13.11 The full set of 2-loop diagrams to be 
considered here is depicted in Figure 13.21 but only about half of these diagrams give 
rise to an imaginary part. It is an easy task to identify this subset of diagrams as 
the generation of an imaginary part is always related to final state interactions. 

The colour factors of the 2-loop diagrams from Figure \32\ can be found in Table IXTl 
The diagrams can be written in terms of the following denominators of propagators 





= {p-q-kf , 


Vis = 
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where p denotes the momentum of the 6-quark (with mass m = rrib), p — q the one 
of the quark to the right of the weak vertex and the quark/antiquark of the emitted 
meson have uq/uq, respectively. The on-shell kinematics is reflected by = and 
p^ = 2p - q = m? . The variable Zf = mf/rrih is related to the diagrams with a closed 
fermion loop. For massless quarks in the loop we simply have Zq = 0, for an internal 
b-quark Zb = 1 and for the case of a charm quark we write z = Zc = rric/mb. 
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Figure 3.2: Full set of non-factorizable 2-loop diagrams. In each diagram the fermion 
line to the left of the four- quark vertex denotes the massive b quark, all other quarks are 
massless. The bubble in the last four diagrams represents the 1-loop gluon self-energy. 
Only diagrams with final state interactions, i.e. with at least one gluon connecting the line 
to the right of the vertex with one of the upper lines, give rise to an imaginary part. 



68 



3. HADRONIC TWO-BODY DECAYS I: IMAGINARY PART 



ODfirator 


^ 1 




line 1-3 (L) 
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Table 3.1: Colour factors of the diagrams in Figure WM , With "L" and "S" we refer to 
the large and small figures, with "A" and "NA" to abelian and non-abelian diagrams. The 
normalization is chosen such that the tree diagram from Figure [777| gives Nc for Q2. 




Figure 3.3: Scalar Master Integrals that appear in our calculation. We use dashed 
lines for massless propagators and double (wavy) lines for the ones with mass m (zfm). 
Dashed/solid/double external lines correspond to virtualities Q/um? /m? , respectively. Dot- 
ted propagators are taken to be squared. 
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Step 2: Reduction to Master Integrals 

The reduction algorithm represents an indispensable tool for the considered calcu- 
lation. It enables us to express all diagrams from Figure 13.21 as linear combinations 
of Mis which are multiplied by some Dirac structures. As the coefficients in these 
linear combinations are real, we may extract the imaginary part of a diagram at 
the level of the Mis which is a much simpler task than for the full diagrams. The 
calculation of the imaginary part involves 14 Mis which are depicted in Figure 13.31 

Step 3: Manipulation of Dirac structures 

We do not perform the bound state projections at this level of the calculation as 
this would yield unwanted traces with 75. We instead treat the two currents in- 
dependently and make use of the equations of motion for the on-shell quarks in 
order to simplify the Dirac structures. In this way we end up with three irreducible 
structures which are given by 



The second structure gives rise to 1-loop evanescent operators according to (13.21) . 
As the last structure only enters in the finite piece of the considered calculation, it 
can simply be evaluated in ci = 4 dimensions without the need to introduce 2-loop 
evanescent operators. 

Step 4: Calculation of Master Integrals 

Some Mis in Figure 13.31 can be solved easily with the help of Feynman parameters. 
This direct approach could be improved with the help of the Mathematica package 
HypExp [On] which allows to expand a special class of hyper-geometric functions 
to arbitrary order in e. However, as discussed in detail in Section 12.41 wide parts 
of the 2-loop calculation were performed with the help of the method of differential 
equations in combination with the formalism of Harmonic Polylogarithms (HPLs). 
We further applied Mellin-Barnes techniques for the calculation of the boundary 
conditions to the differential equations. Apart from the two Mis with an internal 
charm quark (wavy line) , we were able to express all Mis with the help of a minimal 
set of five HPLs given in (12.481) . 

The situation is more complicated for the Mis with an internal charm quark which 
introduces a new scale to the problem. However, a closer look reveals that these 
Mis depend on two physical scales only, namely uml and = z^ml. The Mis can 
then be solved within the formalism of HPLs in terms of the ratio ^ = /u if we 
allow for more complicated arguments of the HPLs as 77 = | (l — a/1 + 4^) . As an 
independent check of our results we evaluated the Mis numerically using the method 
of sector decomposition. The results of the Mis can be found in Appendix lA.ll 



[d'y^ L u] , 



(3.7) 



70 



3. HADRONIC TWO-BODY DECAYS I: IMAGINARY PART 



3.3 Renormalization and IR subtractions 

So far we have computed the (unrenormalized) matrix elements 

(Qi,2) = {MiM2\Qi,2\B) (3.8) 

to NNLO in perturbation theory (without spectator scattering). These matrix el- 
ements are UV and IR divergent. In this section we discuss adequate subtractions 
which will lead to a finite result for the NNLO vertex corrections. 



3.3.1 Renormalization 

The renormalization procedure involves standard QCD counterterms, which amount 
to the calculation of various 1-loop diagrams as the ones depicted in Figure 13. 4[ as 
well as counterterms from the effective Hamiltonian. We write the renormalized 
matrix elements as 



(3.9) 



where = zl^'^Zq^'^ contains the wave-function renormalization factors of the 
massive b-quark Zi, and the massless quarks Zq, whereas Z is the operator renor- 
malization matrix in the effective theory. We introduce the following notation for 
the perturbative expansions of these quantities 



(Qi) 



ren/bare 



fc=0 



^^^=^^^ + E(^) ^ (3-10) 



fc=l k=\ 

and rewrite fl3.9p in perturbation theory up to NNLO which yields 





Wj/bare' 




(^*)barc 




(^«)bare 



7« + 7«A.. 



7(1) , 7(1) r.. 



Wj/bare' 
(^i)bare + 



y{2) 7(1)7(1) , 7(2). 



(3.11) 



(4) 



(0) 



3 1 bare " 



The full calculation thus requires the operator renormalization matrices Z'^^'^-'. For 
the calculation of the imaginary part, the terms proportional to the tree level matrix 






Figure 3.4: Sample of counterterm diagrams. 
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elements do not contribute and Z^'^^ drops out in (13. lip as expected for an effective 
NLO calculation. 

Mass and wave function renormalization are found to be higher order effects. For 
the renormalization of the coupling constant we use 

Z, = l-^^(^^N,-ln,^+0{al). (3.12) 
The 1-loop renormalization matrix Z^^^ can be found e.g. in |lU4j and reads 

Z« = ^ 3 12 9 (3 13) 

V 6 1 y ^ 

where the two lines correspond to the basis of physical operators {Qi,Q2} and the 
four columns to the extended basis {Qi,Q2, Ei, E2} including the mixing of the 
non-physical evanescent operators into the physical ones. 

3.3.2 Factorization in NNLO 

In this section it will be convenient to introduce the following short-hand notation 
for the factorization formula (I1.13P 

{Qi)ren = F-Ti®<^ + ... (3.14) 

where F denotes the B Mi form factor, Tj the hard-scattering kernels and $ 
the product of the decay constant and the distribution amplitude 0m2- The 
convolution in f ll.l3p has been represented by the symbol ® and the ellipses contain 
the terms from spectator scattering which we disregard in the following. 

Formally, we may introduce the perturbative expansions 

OO , CO J CO , 

^-EO^'". ^-E(S)^'". *-E(£) (3.15) 

k=Q k=0 k=0 

Up to NNLO the expansion of (13.140 yields 





= ■ 


^(0) 
i 




{Qi)ien 


= • 


7^(1) 

i 


® ${o) + pW . r.(o) ^ $(0) + ^(0) . pm ^ ^(1) 


(Qi) rcn 


= ■ 


i 


^ ${o) + pw . T.d) ^ $(0) + ^(0) . pm ^ ^(1) 



+ . Tf ) ® $(0) + F(^) ■ Tf ) ® $(1) + ■ Tf ) ® $(2). (3.I6) 
In LO the comparison of (13. lip and (I3.16P gives the trivial relation 

(4)^°^ = (4)il = i^^°^-Tf (3.17) 
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which states that t/"'' can be computed from the tree level diagram in Figure ll.7[ 
In order to address higher order terms we split the bare matrix elements into con- 
tributions from (naively) factorizable (f) and non-factorizable (nf) diagrams 



(3.18) 



(4)£e - mf^ + (4)1?. 

In NLO the corresponding diagrams have been shown in Figure [L9] and Figure [L8 
respectively. To this order (13. lip and (I3.16P lead to 



nf 



(1); 



$(0) _^ ^(0) . 2^(0) 



$(1), 



(3.19) 



(3.20) 



which splits into 

(4)ii^-FZ,^)(4.)(o) = FW-7;«®$(°) 

for the calculation of the NLO kernels t/^"* and 

(g,)« + = i^^'^ ■ ® ^(o) + ■ ) ® <|.(^), (3.21) 

which shows that the factorizable diagrams and the wave-function renormalization 
are absorbed by the form factor and wave function corrections F^^^ and 

This suggests in NNLO the following structure 



= ■ T.^°^ ® + F(^) ■ r.^°^ (g) <l>(^) + ■ T.^°^ (g) $(2). 



(1) 



(3.22) 



(2) 

These terms are thus irrelevant for the calculation of the NNLO kernels which 
justifies that we could restrict our attention to the non-factorizable 2-loop diagrams 
from Figure 13.21 In NNLO the remaining terms from (13.1 ip and (13.160 contain 
non-trivial IR subtractions 



z^ 



(0) 



(3.23) 



This equation can be simplified further when we make the wave function renormal- 
ization factors in the form factor and the distribution amplitude explicit 



p 7I/2 p 

r — Zj^ Zj^ -Tamp; 



z $ 



(3.24) 



Notice that the resulting amputated form factor Famp and wave function $amp con- 

1 /2 3 /2 

tain UV divergences by construction. We recall that — ^ 



Z^^ Zq and find 
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Combining (13.231) with f l3.25p and fl3.20p . we arrive at the Master Formula for the 

(2) 

calculation of the hard-scattering kernels T- in NNLO 



(1) 



(1) 



amp 



amp" 



(3.26) 



The 2-loop matrix elements on the left-hand side of (13.261) have been considered in 
Section [32] and the 1-loop renormalization matrix has been given in (14.41) . The 1- 
loop matrix elements of the non-factorizable diagrams involve the calculation of the 
diagrams in Figure 11.8b . The tree level matrix elements as well as the 1-loop matrix 
elements of the factorizable diagrams from Figure 11.91 can be disregarded here as 
they do not give rise to an imaginary part. Hence, the only missing pieces for the 
calculation of the imaginary part of the NNLO kernels are the IR subtractions 
on the right-hand side of (I3.26p which we consider in the following section. 



3.3.3 IR subtractions 



Let us first address the NLO kernels T'f^ which can be determined from equation 
(13.201) . The renormalization in the evanescent sector implies that the left hand side 
of (I3.20p is free of contributions from evanescent operators up to the finite order 
e*^ . However, as the NLO kernels enter (I3.26P in combination with the form factor 
correction F'^^'> which contains double (soft and collinear) IR divergences, the NLO 
kernels are required here up to 0{e'^). Concerning the subleading terms of 0{e), 
the evanescent operators do not drop out on the left hand side of (I3.20p and we 
therefore have to extend the factorization formula on the right hand side to include 
these evanescent structures as well. Schematically, 



'(l)//^.\(0) 



(1) 



(0) 



(0) T.(l) 



i,E 



(0) 



(3.27) 



with a kernel T^^g = 0{1) and an evanescent matrix element F^'^^^' = 0{e). 
Similarly, the right hand side of (I3.26P has to be modified to include these evanescent 
structures. 

From the calculation of the 1-loop diagrams in Figure 11.8b . we find that the NLO 



^(0)^(0) 



kernels vanish in the colour- singlet case, T2 
part of the colour-octet kernels is given by 



(1) 



(1) 



2,E 



IT 



(1), 



u 



-3 -21nM + 21nM) 1 



0, whereas the imaginary 



+ 



'11 - 31nM 



37r2 



26 



In^ u + \y? u] 



2 + y) lnM+ (^9-y )lnu 



- - In^ u - - (in^ u - In^ u) 
2 3 ^ ^ 
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1 



sL 



1 



TT ■ AN^ y~ ' ' ■ V3 2 

where L = \iafj,'^/ml and we recall that u = l 



\nu- -\nu)e + 0{e^) 



(3.28) 



u. 



Form factor subtractions 

We now address the form factor corrections which require the calculation of the 
diagram in Figure 13.51 (for on-shell quarks) and its counterterm. According to the 
definition of -Famp in f l3.24p . we do not have to consider the wave function renormal- 
ization of the quark fields here. 

We again have to compute the corrections for physical and evanescent operators. 
Concerning the physical operators with Dirac structure [7^^-^] [y^L] the counterterm 
is found to vanish and we get 



amp 



mt 



e{l - 2e) 



pio) $(0) 



(3.29) 



refiecting the 1/e^-singularities mentioned at the beginning of this section. On the 
other hand, the evanescent operators with ['y^^y'^'y^'L] [j^j^'^pL] — 16[7'^L] [y^L] yield 
a contribution proportional to the evanescent and the physical operators. We find 
a non- vanishing contribution from the counterterm diagram in this case and obtain 



^(1) $(0) 
-^amp,E 



mt 



l-e) 



pio) $(0) 



Tie) 



mt 



e(l-2e)(l-£)2 



(3.30) 



The first subtraction term in fl3.26p then follows from combining the form factor 
corrections in f l3.29p and f l3.30p with the NLO kernels in f l3.28p . We emphasize that 

the corrections related to the evanescent operators do not induce a contribution to 

(2) 

the physical NNLO kernel in this case since 



- K 



(1) 



TT 



amp,E 



ImT(>g) 



0{e) 



i?(o)$(o)_ 



(3.31) 



Wave function subtractions 

Concerning the wave function corrections we are left with the calculation of the 
diagrams in Figure 13.61 for coUinear and on-shell partons with momenta uq and uq. 



Figure 3.5: 1-loop contribution to the form factor correction F[ 



(1) 



amp ■ 
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Figure 3.6: 1-loop contributions to the wave function correction $Lmp- The 
dashed line indicates the Wilson-line connecting the quark and antiquark fields. 



However, as in our set-up = all these diagrams vanish due to scaleless integrals 
in DR. We conclude that the wave function corrections are determined entirely by 
the counterterms. We compute these counterterms by calculating the diagrams 
from Figure 13.61 with an off-shell regularization prescription in order to isolate the 
UV- divergences. The counterterm for the physical operators is found to be 

,7/1 = — 

" amp V 



dw V{u,w) 



(3.32) 



with the familiar Efremov-Radyushkin-Brodsky-Lepage (ERBL) kernel [5B]I37] 



V{u, w) 



u 



e{w-u)- 1 



u 



e{u-w)- 1 



(3.33) 



w \ w — U J w \ w — u 

where the plus-distribution is defined as [f{u, w)]^ = f{u, w)—5{u—w) dv f{v, w). 
For the evanescent operators we obtain 

2C, 



dw 



2Ae Ve{u, w) <I>(°) {w) + V{u, w) ' 



w 



(3.34) 



where Ve{u,w) denotes the spin-dependent part of the ERBL kernel given by 



Vi U 

Ve(u, w) = 6(w — u) h 6iu — w) — 

w w 



(3.35) 



Notice that the evanescent operators do induce a finite contribution to the physical 

(2) 

kernel in this the convolution with the corresponding NLO kernel implies 



- Ff ) Im T- $-4,^ 



6C| 

Nr 



+ 0(6) 



ir(o)${o)_ 



(3.36) 



We finally quote the result for the convolution with the physical NLO kernel 



3 u 



\nu 



u 



+ In^ M — 2 In M In M — In^ m — 4Li2 ("U 



+ L 



TT^ 15 ^ 5^ — 4/lnM Inw 

— - — -2C3 + ^— — + 

2 4 2 \ u u 



^ 1 - 



3Li2(-u) In u + 

^ ' 2u 2u 



2 1 3li i2— ^1*^ i2i — 

" ■ ' m M — - m -u + m -umw 



-f ^ ln=^ M + 2 In M Li2(M) + 2Li3(M) + 2Si,2(m) + 0{e) | <I>(°) 



(3.37) 
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3.4 Tree amplitudes in NNLO 

The NNLO kernels can now be determined from the Master Formula fl3.26l) . We 
indeed observe that all UV and IR singularities cancel in the NNLO kernels which 
provides a very important and highly non-trivial cross-check of the calculatior0. In 
analogy to Section [1.2.31 we quote our results in terms of the tree amplitudes ai,2- 

3.4.1 ai in CMM basis 

The procedure outlined so far leads to the colour-allowed tree amplitude in the CMM 
operator basis from (13. 2p . We write 



ai(MiM2) = (72 + 



4^ m'r 



+ 



(3.38) 



where the ellipses denote the terms from spectator scattering which are irrelevant 
for our purposes. In the CMM basis, the imaginary part of the vertex corrections 
can be written in the form 



1 

- Im V^^^ = / du gi{u) 
Jo 



- Im 

TT 



(2) 



du 



/29 2 

(yiVc--r2/)(7i(?i) + Ci./ii(M) 



In At 

mi 



+ Cf h2{u) + Nc h-i{u) + (ri/ - 2) hi{u] 0) + hi{u] z) + hi{u] 1) \<Pm2{u) 



Im yi^^ = 



TT 



du 



Qgi{u) \n^ + ho{u) ^(pMiiu). 
mr 



(3.39) 



We stress that these quantities do not correspond to the \/(^'2) from (11.221) which 
will be given in the following section after the transformation to the QCDF basis. In 
writing (13.391) . we have made the dependence on the renormalization scale explicit 
and we have disentangled contributions that belong to different colour structures. 
The function h/^{u\ zj) stems from diagrams with a closed fermion loop and depends 
on the mass of the internal quark through zj = mf/nib- We write z = Zc = rric/mb 
for simplicity. 



In NLO we find 



g^{u) = -3 ~2\nu + 2\nu. (3.40) 

"'^As we do not distinguish between UV and IR singularities in our calculation, we cannot verify 
their cancellation independently. However, the aforementioned renormalization and IR subtraction 
procedure can be organized in a way that allows to control the cancellation of the leading poles in 
several intermediate steps of the calculation. 
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The NNLO kernels were so far unknown. They are found in this work to be 
ho{u 



155 4 
— + 4(3 + 4Li3(M) - 4Si,2(m) - 121ntiLi2(M) + - In^ m - 6 In^ m Inw 



+ 



u 



-U2{u) 



2u 



In^ M + 



2uu 



In M In M 



4 -11m + 2^2 47r2\, (5 + 6u^ - 12u^)7i^ , 
1 Inw V [u lAA 



u 



2Auu 



+ 



3-m + 7m\ 2 11-10m2_ , , l-14u' 



+ 



2u 
46 - 51m 



In u 



Inu 



u 



hi{u) = 36 + 
h2{u) 



21n^M-4Li2(M) + 



(41 - 42^2)77^ 
2M 

2(13 - 12m) 



Li2 (m) H — In M In M 



l-u 



Au 

- {u ^ u) 
In M — (m m) 



79 + 32C3 - 16Li3(M) - 32Si,2(m) + ^ In^ m + ^^^^Li2(M) 

3 u 

13 - 9m + 6^2 - i7_6^2_g^4 

in M H in M in M 



2m 

5-11m + 2m2 47r 
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4mm 

(1 + 14m2 - 8M^)7r2 



5MM 



+ (m m) 



+ 



2 9 — 14m2 
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3 MM 



+ 



+ 4 
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iuMiuM 



M 



M 



+ y ) in M - 



(26 - 21m 



6m 



— (m ^ m) 



1 379 4 — 
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M 



+ 
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4m 



in M 



4mm 



in M in M 



+ 



41-66m + 8m2 2 



4m 



TT 



(1 + 6m2 - 4M4)7r2 ^ 

iUM H ^ h (m ^ Ml 



mu 



+ 



1 15 — 26m2 

- 2Li3(M) + 4Si,2(m) + 4 in m Li2(M) + - in^ m + — — Li2(M) 

3 4mm 



+ 



11 - 14m - 42m2 



in M — 



11 - 14m2 
4m 



in M in M 
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+ 



/2165- 2156m tt^x (53 -42^2)712 ^ 

— m-u H (u ^ u) 

V 36m 3 / 24m ^ ' 

+ — - 2^2 ^ ^ 2 f (1 + Ai)xi + (1 + Qi)x^ In xi - 4^ xi In X2 
b u X2 \ / 

+ (I - 2(1 + 40x2) In + (^(1 + 20x1 + (1 + 6OX2) In M + (m ^ m) 

94^2 2(1-3^2), ^xi 2[(6 + 29e + 20e')xi + (29 + 380^X2] , 

— m mxi 

9^ 3 X2 94 

4(1 - 3^2) 2mm(6 + 29^ + 20^^)^^ + (1 - 2m) (6m - m^^) 

-Xi m X2 H :^ : m z 



3e 



In M In 2; 



9mm^ 

;i2 + 29e + 2e2)xi + (29 + 380^^2 

9e 



In M — (m ^ m) 



(3.41) 



The last kernel hi{u\ z) has been given in terms of 
1 



l + 4e-l), X2 = \[^l + ^+l 



M 



(3.42) 



In the massless limit hi{u\ z) simply becomes 



, / N 17 2,2 2,2 - 20, 38, „ 

hAu\\)) = m mH — m mH InM InM. 

^'^ 3 3 3 9 9 



(3.43) 



3.4.2 ai and a2 in QCDF basis 

We now perform the transformation of the colour- allowed tree amplitude ai into 
the QCDF operator basis from (13.11) . As discussed in Section [3.11 manifest Fierz 
symmetry in the QCDF basis allows us to derive the colour-suppressed amplitude 
a2 directly from ai under the exchange Ci ^ €2- 



The colour-allowed tree amplitude has been given in the CMM basis in f l3.38p and 
in the QCDF basis in f ll.22p . If we focus on the imaginary part and disregard 
contributions from spectator scattering, these relations become 



— Im «! 

TT 



V 



— Im 



47r A^c 



A-n Nc IT 



Im 



C,v^^) + !±[C,VI 



C2vr' 



(3.44) 



In order to compute V^^^'^), we need the relation between the Wilson coefficients in 
the CMM basis Ci 2 and the ones in the QCDF basis Ci,2 to NLL approximation. 
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The Wilson coefficients can be found e.g. in |103] . where the transformation between 
both bases has been studied in detail. From this, we derive 



Ci = 2C2 + 



as 



4Ci + y ^2 



0{a 



1 ^ 



Combining ^Ml, and fICTD we obtain 



(3.45) 



1 



TT 



Im V^^^ = - Im V(^) 



IT 



du gi{u) (pMiiu), 



IT 



Im v; 



(2) 



Im 



TT 



du 







1 

3 5(1 (m) In — + 2 5(1 (m) + - /lo (m) }> 0Af2 (w) 



- Im v. 
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(2) 



— Im 

1 

du 



3 



7 



6 



+ - h2{u) + 3 /i3(m) + {uf - 2) /i4(m; 0) + h4{u; z) + h^{u] 1)|0m2('")- 

(3.46) 

Notice that these expressions determine the vertex corrections for the colour-allowed 
amplitude ai and the colour-suppressed amplitude 02 according to fll.22p . The 
equations in (13.461) represent the central result of our analysis. The expression for 
is in agreement with (I1.19p . whereas the expressions for IZ/J are new. The 
kernels gi and /10-4 can be found in (13.401) and (I3.4ip . The terms proportional to 
Uf have already been considered in the analysis of the large /^o-limit in [7T|72] . Our 
results are in agreement with these findings. 



3.4.3 Convolution with distribution amplitude 



The NNLO vertex corrections have been given in (13.460 as convolutions of hard- 
scattering kernels with the light-cone distribution amplitude of the meson M2. We 
may explicitly perform the convolution integrals by expanding the distribution am- 
plitude into the eigenfunctions of the 1-loop evolution kernel 



)M2 {u) 



6uu 



l + J2at'^Ci'/^\2u 



n=l 



(3.47) 
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where g^Yid Cn^"^^ are the Gegenbauer moments and polynomials, respectively. 
It is convenient to truncate this expansion at n = 2. The convolution integrals with 
the kernels gi and /iq-s from (13.401) and (13.411) then give 



^0 
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The convolution with h^lu; z) from (13.411) can also be performed analytically 

Hi{z) = [ du h4{u; z) ^A/al"") 

+ 148^2 - QQz^F{z) - 36^^ In^ 



(3.48) 
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(3.49) 
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where we defined 

yi 



F{z) = Lis^-yi) - Si,2(-2/i) - \nyiU2{-yi) + ^ln?/iln^2/2 - ^2^^^ + Cs- 

(3.50) 

In the massless limit the function H/^{z) simply becomes 

i^4(0) = | + 7af + ^af . (3.51) 

The finiteness of all convolution integrals in fl3.48p and fl3.49p completes the explicit 
factorization proof of the imaginary part of the NNLO vertex corrections. 

We summarize our results for the vertex corrections in the considered representation 
of the light-cone distribution amplitude of the emitted meson M2 



-Im \/(l) = -3-3af^ 

TT 



1 ^ "(2) _ o ^1 , ,^ , 1189 , 477r2 /33 23vr2 



- Im V/'^ = 9 (1 + af^) In^ + ^ + ^-SCa- — 
TT ^ ^ ^ ^ m2 24 90 V 8 10 

173 97r 



60 35 



"2 5 



TT ^ V 3 ' y 72 135 3 



00793 166.^ ^j, _ /3155 _ 187^! ) ^ ^ 



V 72 15 y ' V 48 42 

(3.52) 

with H4,{z) given in fl3.49p . In order to illustrate the relative importance of the 
individual contributions, we set fi = rrib and z = rric/mb = 0.3 which yields 

Im = -9.425 - 9.43af^ 

Im Vl^^ = 141.621 + 58.36 af^ - 17.03 

Im V^'^^ = -317.940 - 115.62 af^ - 68.31 (3.53) 

We find large coefficients for the NNLO vertex corrections and expect only a minor 
impact of the higher Gegenbauer moments (in particular in the symmetric case with 
_ g-j_ Notice that all contributions add constructively in ai^2 due to the relative 
signs of the accompanying Wilson coefficients. In the case of ai the contribution 
from v}"^^ is found to exceed the formally leading contribution V^^^ due to the fact 
that the latter is multiplied by the small Wilson coefficient C*2. For a2 the impact 
of the NNLO vertex corrections is also substantial, roughly saying they amount to 
a 50% correction. A more detailed numerical analysis including the contributions 
from spectator scattering will be given in the following section. 
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3.5 Numerical analysis 

We conclude this chapter with a brief numerical analysis (an extended version will be 
given in flOO| ). We first consider the vertex corrections solely without the spectator 
scattering contributions which have been computed recently in [25)126] . These will be 
added in the second part of our analysis which will lead us to the full NNLO result 
for the imaginary part of the topological tree amplitudes in QCD Factorization. 

3.5.1 Vertex corrections 

We come back to the question if the large /3o-limit considered in [7T|7^ represents 
a good approximation for the imaginary part of the NNLO vertex corrections. For 
illustration, we introduce two functions Ci{u) defined by 



The Ci{u) correspond to a combination of Wilson coefficients and hard-scattering 
kernels multiplied by the distribution amplitude of the emitted meson M2. The 
imaginary part of the functions Ci{u) are shown in Figure 1X71 with the asymptotic 
form of the distribution amplitude, for simplicity. 

As we have stated at the end of the last section, we find that the NNLO corrections 
add constructively to the NLO results. They provide the dominant contribution 
to the imaginary part of ai and a substantial contribution to the one of 02- We 
further see that the large /3o-limit is a good approximation in the case of 02 but 
not for «!. This can be traced back to the fact that the imaginary part of V2 is 
reproduced well in this approximation whereas the one of Vi , which provides the 
most important contribution to ai, is missed completely. 

Concerning the scale dependence we recall that the imaginary part has only NLO 
complexity at the considered order in perturbation theory. We therefore use the 
2-loop expression for the running coupling constant (n/ = 5, A^ = 0.225 GeV) and 

O.ip . . . . ^ O.lp ' ' ' ' n 




(3.54) 
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Figure 3.7: Imaginary part of C,i{u) and C2{u) o-s introduced in \3.54-^ . The 
graphs show the 1-loop vertex corrections (blue), the large f3Q- approximation 
(gray) and our new results including the 2-loop vertex corrections (red) for 
/i = rrif, (with asymptotic distribution amplitude). 
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Figure 3.8: Renormalization scale dependence of the imaginary part of ai and 
a2 (without spectator scattering). Blue: NLO result. Red: NNLO result. 

consider the Wilson coefficients in NLL approximation which we take from jl05] . 
Contrary to the second Gegenbauer moment, the first moment would also be required 
in NLL approximation as it enters the expression for V^^^ in fl3.52p . However, in the 
following analysis we focus on the B ^ tttt decays where the first moment is absent 
and we only implement the LL evolution of the second moment which is given by 

= (^)'"'''°«^'^(/^o)' = (3.55) 

The scale dependence of the imaginary part of ai^2 is shown in Figure 13. 8[ We 
observe only a minor reduction of the scale dependence if we vary the scale between 
~ 2.4 GeV and 2mh ~ 9.6 GeV, in particular for ai where the NNLO correction 
dominates over the NLO result. 

3.5.2 Full NNLO result 

We finally combine our results with the spectator scattering contributions from [23] . 
One remark is in order concerning the scale dependence of the spectator term. The 
respective kernel T^^ receives hard and hard-collinear contributions which are en- 
coded in a hard coefficient function and a jet function, respectively (notice that 
the hard coefficient function represents the only source for an imaginary part). In 
the following discussion we simply evaluate all quantities related to the hard coef- 
ficient function at the hard scale and all other quantities at the hard-collinear 
scale fihc (this corresponds to equation (58) of |25] with U\\{fih, fihc) = !)• A more 
sophisticated treatment of the scale issues in the spectator term is relegated to [100] . 

Our input parameters for the B ^ tttt amplitudes are summarized in Table 13. 2[ 
They correspond to the values from previous analysis in QCD Factorization [25l|66] 
with updated values for the form factor and the Gegenbauer moment based on recent 
LCSR analyses In order to estimate unknown perturbative corrections we 

vary the hard scale in the range fj,h = 4:.8~^2a GeV and the hard-collinear scale 
independently between fihc = I-SIq' GeV. 
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Parameter 


Value 


Parameter 


Value 


■"■MS 

rUb 
rric 
f. 


0.225 
4.8 
1.3 ±0.2 
0.131 


fB 

Ff"*"(0) 
aJ(lGeV) 


0.2 ±0.03 
0.26 ±0.04 
0.35 ±0.15 
0.2 ±0.2 



Table 3.2: Theoretical input parameters (in units of GeV or dimensionless) . 



The complete NNLO result for the imaginary part of the topological tree amplitude 
is found to be 

Im ai(7r7r) = 0.012|^(,) ± 0.03l|^(2) - 0.019|^,2) 

= 0.025 ±0.021, 
Im a2(7r7r) = -0.077|^(i) - 0.052|^(2) ± 0.03l|^(2) 

= -0.098 ±0.035, (3.56) 

where we disentangled the contributions from V'^^\ V'^'^'> and if*^^-* according to (11.221) . 
In the case of ai the NNLO corrections exceed the NLO result which can be ex- 
plained by the fact that the latter is multiplied by the small Wilson coefficient C*2- 
In both cases the individual NNLO corrections are found to be sizeable, but we 
observe a large cancellation in their sum. The NNLO vertex corrections considered 
in this work turn out to dominate over the spectator terms resulting in a moderate 
additive contribution to the NLO (BBNS) result. 

The uncertainties quoted in (13.561) stem from the variation of the parameters shown 
in Table 13.31 As the dominant sources we identify the hadronic parameters A b and 
02 . Moreover, the sensitivity to the renormalization scale remains sizeable at NNLO 
as we have mentioned at the end of the last section. We finally emphasize that we 
have not yet assigned an error estimate to unknown power corrections which will be 
included in |100] . 







fJ'hc 


rUc 


fB 


pB. 


Xb 


^2 




0.011 


0.006 


0.000 


0.003 


0.003 


0.014 


0.008 


Ci2 


0.019 


0.009 


0.000 


0.005 


0.006 


0.023 


0.013 



Table 3.3: Uncertainties in our predictions of the imaginary part of ai{Tnr) 
and a2(7rvr) from the scale variation and the input parameters in Table COl 



Chapter 4 



Hadronic two-body decays II: 
Real part 

The calculation of the real part of the NNLO vertex corrections in hadronic two- 
body decays proceeds along the same lines as the one of the imaginary part which we 
presented in the previous chapter. However, we will see below that the calculation 
of the real part is far more complicated involving many additional Mis and the full 
NNLO complexity concerning e.g. the issue of renormalization and the treatment of 
evanescent operators. 

In this chapter we present a preliminary result for the real part of the colour-allowed 
tree amplitude ai (in the CMM operator basis). Similar to what we have seen in 
the previous chapter, the colour-suppressed tree amplitude a2 can then be derived 
from «! after the transformation into a Fierz-symmetric operator basis (which we 
called QCDF basis in Chapter E]). In NNLO we thus have to extend the basis from 
f l3.ip to include 2-loop evanescent operators with an appropriate definition of e- and 
e^-terms which guarantees manifest Fierz symmetry in this operator basis. As we 
have not yet worked out the details of this last step, we refer to |106j for the analysis 
of the colour-suppressed tree amplitude. 

4.1 2-loop calculation 

The 2-loop calculation will be performed in the CMM operator basis (cf. Section ISTTI) . 
Apart from the operators in (13.21) . we have to take into account 2-loop evanescent 
operators which are defined by 

E[ = [u^^YYYYL TH] [d^.iulplalr L T^u] - 20^1 - 256 Q,, 

E'^ = [u-f'^YYrYL b] [d-f^-fulplalr Lu]- 2OE2 - 256 Q2- (4.1) 

We turn to a brief characterization of the considered 2-loop calculation following 
our recipe from Section 12.11 
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Step 1: Set-up for loop calculation 

In contrast to the calculation of the imaginary part from Chapter [Sj we now have to 
consider the whole set of non-factorizable 2-loop diagrams from Figure I3.2[ Notice 
that the most complicated diagrams only enter the calculation of the real part. 
Whereas the diagrams that we considered for the calculation of the imaginary part 
contained at most one massive (6-quark) propagator, we now have to deal with up 
to three massive propagators. 

Step 2: Reduction to Master Integrals 

The fact that the diagrams involve between 0-3 massive propagators immediately 
leads to many distinct topologies and a large number of Mis. In addition to the 14 
Mis from Figure [331 we find 22 (real) Mis which are shown in Figure HUl We further 
remark that our Mathematica implementation of the reduction algorithm hardly 
succeeds to reduce the most complicated diagrams of the considered calculation 
within a reasonable amount of CPU time. 

Step 3: Manipulation of Dirac structures 

Concerning the irreducible Dirac structures we find the same set f l3.7p as in the 
calculation of the imaginary part. In contrast to our analysis from Chapter [3], the 
last structure now enters the divergent piece of the calculation giving rise to 2-loop 
evanescent operators according to (14. ip . 




Figure 4.1: Additional Master Integrals that have to be considered for the calculation of 
the real part of the NNLO vertex corrections. For details cf. Figure HOI 



4.2. RENORMALIZATION AND IR SUBTRACTIONS 
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Step 4: Calculation of Master Integrals 

In the last step of the calculation it becomes obvious that the analysis of the real 
part of the NNLO vertex corrections is much more complex than the one of the 
imaginary part. We now have to calculate the real parts of the Mis from Figure 
as well as the Mis from Figure 14.11 in general up to five orders in the e-expansion. 
In order to tackle this highly challenging task we applied the same techniques as in 
our calculation from Chapter [3l It turns out that almost all Mis can be expressed 
in terms of the set fl2.48l) supplemented by the following HPLs of weight w = 4 

H{0,0,0, l;u) = Luiu), 

i/(0,0,l,l;M) = S2,2(w), 

i7(0,l,l,l;n) = Si,3H. (4.2) 

Unfortunately, some Mis with two and three massive propagators do not fit into this 
pattern. We further have to include HPLs related to the parameter —1 and to the 
argument u — 1, cf. ( 1A.43I) - ( lA.SSj) . Still, we find two functions which we could not 
express in terms of Nielsen polylogarithms: H[0, —1, 0, 1; u) and A{u) from f[09p . 



Similar to what we have seen in Chapter [3], the situation is more complicated for the 
Mis which stem from the diagrams with a closed fermion loop. So far we have not 
yet calculated these Mis with a massive internal quark (wavy lines), i.e. we work in 
an approximation which treats all quarks in these loops as massless. We stress that 
this is certainly inconsistent for an internal 6-quark with respect to the remainder of 
the calculation. However, we expect that this will only have a minor impact on our 
(preliminary) result. A consistent treatment of these diagrams is relegated to [106] . 
The analytical results for all other Mis can be found in Appendix IA.2[ 



4.2 Renormalization and IR subtractions 

We have shown in detail in Section 13.31 how to perform appropriate UV and IR 
subtractions in order to extract the NNLO kernels from the hadronic matrix 
elements. The essence is summarized in the Master Formula f l3.26p . We briefiy 
address the issues that go beyond our analysis of the imaginary part from Chapter [3l 



4.2.1 Renormalization 

The renormalization of the 6-quark mass becomes relevant in the considered calcu- 
lation. We treat the 6-quark in the on-shell scheme according to 



"'^The appearance of the function A{u) seems to be an artefact of our calculation, i.e. of our 
special choice of the Mis. Though it enters two of our Mis from Figure it drops out in their 
sum in the respective Feynman diagrams and is therefore irrelevant for our purposes. 
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Due to the full NNLO complexity, we require the 2-loop renormalization matrix Z*^^) 
which we take from |104j 



41 
3 


58 
9 


125 
36 


73 
54 


19 5 
96 108 


-29 


4 


73 
12 





5 1 
24 9 


317 
36 


515 
54 


4493 
864 


49 
648 


1 

384 


349 
> 12 


3 


1031 
144 


8 
9 


35 
192 



e 



2 



_j_ 36 54 864 648 384 864 _ (4 4) 




where the two lines correspond to the physical operators {Qi, Q2] and the columns to 
the extended operator basis {Qi, Q2) Ei, E2, E[, E2} including evanescent operators. 



4.2.2 IR subtractions 

In order to perform the IR subtractions from the right hand side of (13.261) . we must 
compute the real parts of the NLO kernels T-''^^ to 0{e'^). The kernels vanish in the 
colour-singlet case T2^^ = tI^^, = 0, whereas the colour-octet kernels are found to be 

Re t}'\u) = ^^{{-QL + to{u)) {l + eL + ^s'L') 

+ (-L3 + t2(M))£2 + 0(£3)|, 
Re Tl%iu) = (2^ + iE,oiu)) (l + sL) 



(4.5) 
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M M / 12m 

In combination with the form factor corrections from (I3.29P and (13.301) . this deter- 
mines the first subtraction term in (13.261) . In the second subtraction we require the 
convolution of the NLO kernels with the wave function corrections from (13.321) and 
(Km . We find 



Fr Re n:> $1:^,^, - ^{12^ + t^A^) + o(^)} f^'^ $(°), (4.7) 
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1+m-3m2 2 -10m + 6^2 i_ 5^ + 3^2 

H z Si2(m) H r lnMLi2('u) z muLi2{u) 

uu ' uu uu 
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4.3 Tree amplitudes in NNLO 

The NNLO kernels Tf ^ follow from flX^ and are indeed found to be free of UV 
and IR singularities. We emphasize that this provides a very powerful check of 
our calculation which involves the cancellation of poles up to l/e"^ (1/^^) foi' the 
calculation of the real (imaginary) part. 



4.3.1 ai in CMM basis 



We now present preliminary results for the real parts of the NNLO vertex corrections. 
Our results are still preliminary in the sense that the calculation is not yet complete 
(massive fermion loops are still missing) and we have not yet performed numerical 
checks of all Mis. In analogy to fl3.39p we write 
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The NLO kernel is found to be 
3(1 -2m) 
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Concerning the NNLO kernels h^-g we do not quote the expressions for h^, hj and 
hs here, as they are extremely complicated and we have not yet expressed them in 
terms of a minimal set of HPLs. On the other hand, the expressions for and for 
hg (in the massless case) are much simpler and given by 
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(4.11) 



4.3.2 Convolution with distribution amplitude 

We now perform the convolution integrals by expressing the distribution amplitude 
of the emitted meson M2 in terms of its Gegenbauer expansion f l3.47p . We obtain 
analytical results for the convolutions with the kernels g2, and /ig 



du g2{u) (f)M2{u) 



45 , 11 M, 21 

«i 20 2 



du hRiu) 



^M2 



M = 348 



2 2^ 
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(4.12) 



and computed the remaining convolution integrals numerically 

M2 



[ du hiu) (pM2{u) = 322 - 213 af ^ + 3.8 
Jo 



du h7{u) (j)M2{u) = 731 — 348 flj^ 



M2 „M2 



du h^{u) 0m2(m) 



-409 + 412 af 2 - 32 a. 



,M2 



(4.13) 



The cancellation of all singularities and the finiteness of all convolution integrals 
completes the explicit factorization proof of the NNLO vertex corrections. 

We finally collect all contributions and illustrate the relative importance of the 
individual vertex corrections setting /i = 



Re = _22.5 + 5.5 af ^ - 1.1 > 



,M2 

'2 ' 



Re Vl^^ = -148 + 606 af' - 80 a 
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(4.14) 
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4.3.3 Preliminary numerical result 

We conclude this chapter with the presentation of a numerical result for the real 
part of the NNLO vertex corrections. We stress again that this corresponds to a 
preliminary result which treats the c-quark and the 6-quark in the closed fermion 
loops as massless quarks. If we reconsider our results for the imaginary part in this 
approximation, we find deviations of ~ 5% of the individual NNLO contributions. 
As the NNLO terms are subleading for the real part of the colour-allowed tree 
amplitude, we expect that this approximation will have only a minor impact here. 

Our preliminary result for the real part of the NNLO vertex corrections reads 

Re ai(7r7r) = 1.0l|^(o) + 0.03|^(i) + 0.03|^(2) 

= 1.06, (4.15) 

where we have used Wilson coefficients in NLL approximation for simplicity (they 
are indeed known to the required NNLL accuracy and can be found in |1U7] ). As 
expected, the contribution is of minor importance in absolute terms. However, 
the NNLO corrections are found to be as important as the NLO terms which are 
numerically suppressed by the small Wilson coefficient. Interestingly, the vertex 
corrections add again constructively and, as can be seen in comparison with the 
results from [23], come again with the opposite of the spectator interactions. 

The colour-suppressed tree amphtude a2(7r7r) is phenomenologically more interest- 
ing as the respective QCD Factorization prediction is rather low for a satisfactory 
description of the experimental data. In order to derive the NNLO result for ^2 
we still have to solve some conceptual aspects concerning a Fierz-symmetric defi- 
nition of (2-loop) evanescent operators. We therefore relegate the discussion of the 
colour-suppressed tree amplitude to |106] . 
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Chapter 5 



Heavy-to-light form factors for 
non-relativistic bound states 



In this chapter we investigate transition form factors between non-relativistic QCD 
bound states at large recoil energy. Assuming the decaying quark to be much heavier 
than its decay product, the relativistic dynamics can be treated according to the 
factorization formula (11.281) for heavy-to-light form factors obtained from the HQE 
in QCD. In contrast to the B ^ tt transition, the form factors can be calculated 
entirely in perturbation theory in the non-relativistic approximation which allows 
for an explicit analysis of the factorization formula. We perform a NLO calculation 
based on the methods that we developed in Chapter [2] and look for an interpretation 
of the results from the viewpoint of QCD Factorization. 

The basic idea of this work has already been presented in |108j . We emphasize that 
the formalism which we develop in this chapter can be applied for Be —>■ (cc) tran- 
sitions as all these particles can be considered approximatively as non-relativistic 
bound states. Notice that we treat the charm quark as a light quark in this case, 
rric <C mb- We will adopt this terminology throughout this chapter although we 
consider this analysis rather as a toy model for the B ^ n transition. A phenomeno- 
logical study of various Be decays in QCD Factorization will be given in |109j . 

5.1 Non-relativistic approximation 

The wave function for a non-relativistic (NR) bound state of a quark and an anti- 
quark with respective masses mi and m2 can be obtained from the resummation 





Figure 5.1: Resummation of potential gluons into a non-relativistic Coulomb 
wave-function (details of the resummation can be found e.g. in 11 lU}/ ). 
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of NR (potential) gluon exchange as sketched in Figure 15.11 The solution of the 
corresponding Schrodinger equation with Coulomb potential yields 

^5/2 

^c(p) (5.1) 

where k = mrCesCp and = 77111712/ [mi + 1712) is the reduced mass. The normal- 
ization of the wave function gives the (non-relativistic) meson decay constant 

JNR - 7 \ TTT^- (5.2) 

TT (mi + 7712) ' 

In this approximation, the meson is entirely dominated by the two-particle Fock 
state built from a bottom quark with mass M = rub and a charm antiquark with 
mass m = rric- Consequently to first approximation in the NR expansion, the Be 
meson consists of a quark with momentum Mw^ and an antiquark with momentum 
mWfj^, where is the four- velocity of the Be meson [w"^ = 1). The spinor degrees of 
freedom for the B^ meson in the initial state are represented by the Dirac projector 
PH = i(l+V^)75. 

Similarly, a pseudoscalar rjc meson is interpreted clS 3j CC bound state where both 
constituents have approximately equal momenta mw'^, where w'^ is the four-velocity 
of the rjc meson {w''^ = 1). The Dirac projector of the rjc meson in the final state is 
given by Vl = |(1 - V^Ots- 

In the following we consider heavy-to-light transitions at large recoil energy assuming 
M ^ m and working in leading power of the HQE in m/M. The QCD dynamics is 
then described by the SCET degrees of freedom from Table 11.11 where we identify 
the typical hadronic scale with the mass of the t]c meson ~ 2m = 0{m) in the 
NR approach. The relevant scales in the process are thus given by 



^ M > fihc ~ V Mm > fis,c ^ "m > /"NR ~ "mv, (5.3) 

where the NR scale refers to the virtuality of the potential gluons from Figure 15.11 
and f -C 1 is a NR velocity. Notice that all relativistic scales are perturbative in our 
set-up as M, m ^ Aqcd- The relativistic dynamics in the heavy-to-light transition 
can therefore be analyzed in perturbation theory. 

The momentum transfer can be approximated as 

g2 = {^Mbw - Mr^w'f ~ - AMui w ■ w' (5.4) 
which implies a large relativistic boost 

l^w-w' = = 0{M/m) (5.5) 



for large recoil energies with — = 0{M'^). 
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5.2 Perturbative calculation 

For simplicity, we focus on the — > rj^ transition and consider the form factors 
F^{q'^), F_{q^) and FT^q^) which can be defined in analogy to fll.271) . According to 
the general discussion for heavy-to-light decays at large recoil, we have to consider 
hard, hard-collinear, collinear and soft gluon exchange in order to describe the rela- 
tivistic dynamics of the transition form factors, while the non-relativistic modes are 
contained in the bound state wave functions of the initial and final state mesons. 



5.2.1 Tree level 

We have to require at least one relativistic gluon exchange in the large recoil case in 
order to rearrange the quark- ant iquark pair in the final state into a NR configuration. 
Consequently, the diagrams from Figure 15.21 contribute in LO of the perturbative 
expansion which imply the exchange of a hard-collinear gluon with virtually 

[mw — mw'Y — — 27771^ = 0{Mm). (5-6) 

The result for the form factors at LO becomes 



l + 2s 2s 2s 



(5.7) 



where and are the non-relativistic decay constants of the initial and final 
state mesons respectively and we defined 

0{1) (5.8) 



47m — g2 

with s = 1 at maximum recoil q^ = 0. 

From (15.21) . we read off that the decay constants scale as ~ {nivf^M-^'"^ and 
~ mv^l"^ in the combined NR and HQE. We find that the form factors scale as 



Ft'^iq') ^ a,v' i^-) (5.9) 

in agreement with the expected scaling Fi{q^) ~ M~^/'^ from the general discussion 
for heavy-to-light form factors at large recoil 



I 
I 




Figure 5.2: Tree level diagrams. The dot denotes the weak vertex mediating 
the heavy-to-light transition, the lower line the light spectator antiquark. NR 
gluons from the bound state wave functions are not drawn. 
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5.2.2 1-loop calculation 

The tree level results derived above do not give rise to endpoint singularities as 
the wave functions of the non-relativistic bound states have vanishing support at 
the endpoints. For a deeper understanding of the factorization formula we therefore 
have to consider the NLO contributions. We will see below that our NLO calculation 
indeed reveals the full complexity of the factorization formula. We now switch to 
the technical part of the calculation and evaluate the corresponding 1-loop diagrams 
following our strategy from Section 12.11 

Step 1: Set-up for loop calculation 

At NLO the form factors receive contributions from various 1-loop diagrams as 
shown in Figure The corresponding colour factors are summarized in Table EUl 
All diagrams from Figure can be expressed in terms of the following denominators 
of particle propagators 



Pi 


= {Mw + mw - kf - M^, 


Vs-- 


= (Mw + mw — 


mw' - 


-kf-M\ 




= {Mw - mw' + kf - 




= (Mw - kf - 


M^, 




^3 


= (2mw' — kf — m^, 


Vio 


= (2mw' — mw 


-kf 


— m^, 


V, 


= (mw' — mw + kf — m^, 


Vn 


= (mw' — k f - 


- m^, 




V, 


; 2 2 

= k — m , 


Vl2 


= k\ 








= (mw — kf , 


Vi3 


= (mw — kf — 


m^, 






= (mw' — kf , 


Pl4 


= (2mw' — kf 




(5 


step 2: 


Reduction to Master Intej 


^rals 









The reduction procedure is not as efficient in this calculation as it was in the 2-loop 
calculation from Chapter [S] and HI However, we were able to express about 200 scalar 
integrals in terms of 23 Mis which are summarized in Figure 15. 4[ The number of 
Mis is comparably high in this calculation due to the fact that we deal with three 
different types of propagators (with masses 0,m, M) and four external momenta 
(two of them linearly independent w,w'). 

Step 3: Manipulation of Dirac structures 

We project onto the NR bound states with the help of the projectors Vh and Vl 
which we specified in Section I5TT1 The fiavour singlet diagrams in the last two lines of 
Figure Em require special care as they involve traces like Tr(7'^7'^7^7°"75) which may 
invalidate the treatment of an anticommuting 75 within DR. However, the fiavour 
singlet diagrams turn out to give finite contributions and the traces can safely be 
calculated in d = 4 dimensions. 
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Figure 5.3: NLO diagrams. The bubble in the diagrams from the sixth line represents the 
1-loop gluon self-energy. The diagrams from the last two lines contribute to flavour singlet 
final states only. 
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Diagram 


line 1-2 


line 3-4 


line 5 


line 6 


line 7-8 


Colour 


^2 


(^F 2 


2 


Cf 


Cf 
2 



Table 5.1: Colour factors of the diagrams in Figure [531 The normalization 
is chosen such that the tree diagrams from Figure [5!^ give Cf- 




Figure 5.4: Scalar Master Integrals that appear in the 1-loop calculation. We use 
dashed/solid/double internal lines for propagators with masses 0/m/M, respectively. 
Double/solid/dashed/wavy external lines correspond to hard/hard- collinear/collinear/ soft 
external momenta. Notice that we refrain from associating all external scalar products to 
the Mis in this case. The figures therefore do not define the Mis unambiguously in our 
representation (for the explicit definitions of the Mis we refer to Appendix \A.3\) . 
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Step 4: Calculation of Master Integrals 

Though most of the Mis in Figure 15.41 correspond to apparently simple 1-loop 
3-topologies, the calculation is non-trivial due to the fact that the integrals involve 
several distinct scales and many massive propagators. The calculation simplifies as 
we only require the leading power in m/M of the Mis which we extract using the 
method of regions and Mellin-Barnes techniques as described in detail in Section [231 
The analytical results of the Mis can be found in Appendix IA.3[ 

5.2.3 Form factors in NLO 

In NLO we have to take into account the 1-loop diagrams from Figure 15.31 as well 
as standard (tree-level) counterterm diagrams like the ones depicted in Figure 15.51 
We use MS-scheme renormalization constants 

for the coupling constant and the gluon field, respectively. The (non-relativistic) 
quark fields are conveniently treated in the on-shell scheme with 



Zq{m) = 1 - 
Zmim) = 1 - 




(5.12) 



where we indicated that the wave function renormalization constant contains IR 
divergences. The renormalization of the weak vertex involves the Z-factor of the 
heavy-to- light current qV Q which is given by 

Zr = 1 + + 0(a2) for F = (5.13) 

whereas Z-p = 1 for the conserved vector current with F = 7^^. 

Adding up all diagrams and counterterms at NLO, all UV and IR divergences can- 
cel in the form factors as expected. Our explicit NLO results are summarized in 




Figure 5.5: Sample of counterterm diagrams. 
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M=5, m=1.5 



lOr 

8 
6 

4 ■ 
2 



M=5, m=0.3 




20 



15 



10 ■ 



1 2 



0.5 1 1.5 2 2.5 3 



Figure 5.6: Renormalization- scale dependence of the ratio R for different 
choices of quark masses. The solid (dashed) line denotes the NLO (LO) result. 
Left: Realistic case relevant for the — > ??c transition. Right: Toy example 
for the B ^ IT transition. 



Appendix [Bl Let us quote the result for the form factor at maximum recoil 

= here to illustrate the structure of the NLO contributions. We find 

^NLO 



(0) 



Mm 



Arc 



— Ac Hf \m Uf + 04 

3 3 V mM 9 ^ ^ 



Cf - In' — + 



M 35 - 201n2, M 



m 

- In^ h 
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6 

1 + ln2 



In — + - ln'2 + 31n2 

m 3 

. M 1,2 4 , 

In — + - In^ 2 - - In 2 - 

m 3 3 



77r2 
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6 


57r2 


73 \ 


36 





(5.14) 



where M and m are to be considered as pole masses. In (15.141) we isolated the con- 
tribution from the flavour singlet diagrams which is given by 5+ = 1/3 ln(M/4m). 

We now investigate the residual renormalization-scale dependence of Fj^^'-'(O). For 
illustration we will use two sets of quark masses with 



set 1 : M = 5 GeV, m = 1.5 GeV, 
set 2 : M = 5 GeV, m = 0.3 GeV, 



ViWm- 2.7 GeV, 



VMm ~ 1.2 GeV (5.15) 



which correspond to the physical decay rjc and a toy model for the B ^ tt 

transition, respectively (in the second case we thus drop the singlet contribution). 
We use 1-loop running of as with Uf = 4 and 0:^(5 GeV) = 0.2 and study the ratio 



R 



Mm 

f NR f NR 
J M Jm 



F40). 



(5.16) 



The result is illustrated in Figure ESI In both cases we observe a significant enhance- 
ment from the NLO contribution. For the B^ — >■ rjc case (set 1) the enhancement of 
the tree level result at the hard-collinear scale is about 35%. We find a substantial 
improvement of the scale dependence from about 40% at LO to 20% at NLO if 
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we vary /i between the soft scale m and the hard scale M. The hierarchy of scales 
M ^ m is not too large in this case, such that the formally large logarithms In M/m 
do not spoil the convergence of the perturbation series too badly. 

In the toy case for the B tt form factor (set 2), the soft scale m is close to the 
QCD scale Aqcd and consequently the convergence of the perturbative expansion 
breaks down at small scales. The theoretical error due to the renormalization scale 
is not under control in this case. From the conceptual viewpoint we are particularly 
interested in the structure and resummation of the formally large logarithms In M/ m 
in the B ^ n case where they should be counted as In M/m ~ 1. We therefore 
examine the origin of these logarithms in the subsequent section in detail. 
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5.3 Factorization Formula 

We now come to the interpretation of our explicit NLO calculation in terms of the 
factorization formula for heavy-to-light form factors at large recoil energy 

/•oo I* 1 

F,(g2) ^ ^(^2) ^ / ct>B{uj)Ti{u,uj-q^)<P^{u). (5.17) 

Jo Jo 

We start with a closer look at the factorization formula following [76]. The two 
terms in (15.1 7p are associated to the matrix elements of two distinct SCETj operator 
structures, the so-called A-type and i?-type operators (cf. (11.291) for notation) 

Jb = i^hcWhc) {WliIp±hcWhc) K. (5.18) 

The matrix element of the operator Ja defines the form factor ^(g^) and the Hi{q^) 
denote the corresponding matching coefficients from QCD to SCETi which include 
the contributions from hard momentum fluctuations. The operator Jb defines a 
non-local form factor S(r; g^) with short-distance coefficients Ci{q^]T). In the first 
matching step the factorization formula thus takes the form 

FM') ^ HM') i{q') + fdr CM': r) H(r; q'). (5.19) 

The non-local form factor S(r; q^) can be factorized further in a second matching 
step from SCETj to SCETn including the effects from hard-collinear fiuctuations. 
This results in a convolution of a perturbative jet function J(r; m, ui] q^) with leading 
twist distribution amplitudes given by 

/•oo /»1 

9^) = / duj I du (j)B{^) J{T]u,uj]q^) (j)r,{u). (5.20) 

Combining (15.191) and (I5.20p . we see that the hard-scattering kernel from (15.171) can 
be identified as 

Ti{u,uj;q'^) = / dr Ci{q'^;T) J{T;u,uj;q'^). (5.21) 
Jo 

In our analysis we calculate the leading twist distribution amplitudes of the Be and 
rjc meson at leading power in the NR approximation and including perturbative 
(relativistic) effects to NLO in QCD. We then use our explicit NLO results of the 
form factors together with the known expressions for Hi, Ci and J to extract the 
overlap-contribution ^(g^) in our NR set-up which may then be analyzed further 
concerning its factorization properties in SCETn. This immediately leads us to 
the problem of endpoint singularities which arise in our calculation at NLO of the 
perturbative expansion. We find it instructive to demonstrate how these endpoint 
singularities enter our calculation before looking at the factorization of the heavy- 
to-light form factors in detail. 
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5.3.1 Endpoint singularities 

The appearance of endpoint singularities is related to the issue of large logarithms 
which we already mentioned at the end of the last section. As these logarithms may 
spoil the convergence of the perturbation series, we are particularly interested in 
the resummation of such logarithms within the effective theory. The structure of 
formally large logarithms is particularly complicated in heavy-to-light decays due to 
the presence of various scales and degrees of freedom. 

We encountered these logarithms in our explicit NLO calculation in form of In 7 
where 7 = 0{M/m) is the large relativistic boost between the meson rest frames 
(cf. our results in Appendix [B]). We find that the origin of the leading double 
logarithms is related to two different mechanisms which we consider in detail in the 
remainder of this section: The first one gives rise to (standard) Sudakov logarithms, 
the second one is related to endpoint singularities. 



Sudakov logarithms 

Sudakov logarithms appear naturally in processes with (hard-)collinear and soft 
degrees of freedom.. In our calculation they arise for instance from the diagram in 
Figure ISTTl (in Feynman gauge). Rather than looking at the full diagram here, we 
profit from our reduction algorithm and consider an underlying (scalar) MI which 
allows for a more transparent presentation of the interesting aspects. 

The MI to be considered here is 



J ri H% 

(5.22) 

with the Vi given in (15.101) . Notice that the integral contains a double logarithm in 7 
and that we suppressed a prefactor — l/(87mM) for simplicity, cf. (lA.67p . The inte- 
gral may be calculated with Mellin-Barnes techniques as described in Section I2.4.4[ 

In order to understand the origin of the double logarithm, we disentangle the con- 
tributions from different momentum regions of the loop integration. Following the 
method of expansion by regions which we described in Section 1!^. 4. 3[ the MI is found 
to receive leading contributions from the hard, hard-collinear and soft region. 




Figure 5.7: Origin of Sudakov logarithms (in Feynman gauge). 
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The hard momentum region with k ~ (1, 1, 1) is found to give 




6 

(5.23) 



with double poles in e which arise in the limits k± ^ and A;+ — 0. The contribution 
from the hard-collinear region with k ~ (1, A, A^) becomes 



where the singularities stem from A;_|_ oo and A;_ — 0. Finally, the soft region 
with k ~ (A^, A^, A^) generates poles for k± ^ oo and /c+ — > cx) and reads 



Adding up these contributions, all singularities drop out and we reproduce f l5.22p . 

From the viewpoint of QCD Factorization, the hard contribution from fl5.23p is to 
be considered as part of a QCD SCETj matching calculation, the hard-collinear 
effects from (15.241) contribute in the SCETj — > SCETn matching procedure and 
the soft contribution from (15.251) to the (perturbative) calculation of a SCETn ma- 
trix element. Let us think of all regions to be MS-subtracted (the poles contain 
information about the anomalous dimensions of operators in the effective theory). 
We see that the first matching step in (15.231) is free of large logarithms if we choose 
the matching scale fi ~ M. We may then evolve the scale down to ~ fihc using 
RGEs in SCETi which implicitly resum logarithms of the type In M/fi. At the 
hard-collinear scale we perform the second matching step and we see that the choice 
/i ~ (Mm)^/^ guarantees the absence of large logarithms in (I5.24p . The procedure 
continues summing logarithms of the type In Hhc/f^ with the help of corresponding 
RGEs in SCETn and evolving the scale down to n ^ m. In the B —>■ n case the 
SCETn matrix elements have to be calculated with the help of a non-perturbative 
method. As the soft scale is still perturbative in our NR approach, the corresponding 
SCETn matrix elements are calculable giving (I5.25p . 

The aforementioned factorization and resummation procedure is of course oversim- 
plifying as we considered a single integral instead of operators in the effective theory. 
However, it allowed us to emphasize an important feature concerning the resumma- 
tion of Sudakov logarithms: We saw that we were able to choose the factorization 
scale /i such that there are no large logarithms in each matching calculation and we 
pointed out that the resummation of logarithms can be performed using standard 
RG-techniques varying between two well-separated scales which correspond to the 
virtualities of the respective modes. We will see in the following that part of the 
double logarithms in our explicit NLO results are of conceptually different origin. 




(5.24) 




(5.25) 
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Non-factorizable logarithms 

In a second example we consider the pentagon diagram from Figure 15^ (in Feynman 
gauge, see also the discussion in |108j ). For simplicity, let us again focus on a related 
MI here given by 

which may be calculated with Mellin-Barnes techniques. Notice that the MI gives 
rise to a double logarithm and that we dropped a prefactor —l/{2'-jm?), cf. ( 1A.76I) . 

At leading power the MI receives contributions from the hard-coUinear, coUinear 
and soft momentum region. The hard-collinear contribution with k ~ (1,A, A^) is 
found to be divergent for fc^ — > and and reads 



+ + (5.27) 



,'hc hc\ 1 1 , 1 o U 

---\ ) — ~ — + - In „ + - In^ — 

e 2 27m 

This looks very much like (I5.24p but all poles are of IR origin in this case. 

The coUinear region with k ~ (1, A^, A'^) is found to be divergent for — ^> 00 and 
k_ ^ and similarly the soft region with /c ~ (A^, A^, A^) for kj_ ^00 and k^ ^ 0. 
In contrast to what we have seen so far, it turns out that the longitudinal integrations 
along fc_ (fc+) in the coUinear (soft) region are not regularized in this case (notice 
that DR only influences the integration over transverse momenta in d± = 2 — 2e 
dimensions). However, we may introduce an additional regularization procedure in 
order to render these contributions finite. Following [5^108j . we apply an analytic 
continuation replacing V^^ — (— z/^)'^Pg"^~'^ such that endpoint singularities for 
fc_ — > (fc+ 0) show up as poles in 1/6. We thus find 

2 TW^ 

1 , 










m'^ J 


1 


- - In — 




e 











+ (- + ln£^) (- + ln— )+- + 0(.). (5.29) 




Figure 5.8: Origin of non-factorizable logarithms (in Feynman gauge). 
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Notice that the dependence on the ad-hoc parameters 6 and u disappears when we 
consider the sum of colhnear and soft momentum regions 

/ ^ hc^^ ' he 

- 2'"^+ U^'"^) + - + 

and that we reproduce the full result from f l5.26p by adding up fl5.27p and f l5.30p . 

To summarize, we have seen that we cannot disentangle colhnear and soft effects 
in the considered integral as the respective contributions suffer endpoint singular- 
ities related to the integrations over longitudinal light-cone momentum fractions 
which are not regularized in DR (at fixed transverse momenta/virtuality). How- 
ever, the sum of colhnear and soft momentum region turns out to be unambiguously 
well-defined. We may now ask the question if we can factorize the perturbative 
hard-collinear effects in (15.271) from the remnant SCETn matrix element in (I5.30p 
summing all formally large logarithms into short-distance coefficient functions. 

We repeat our analysis from the last section and think of the SCETi — SCETu 
matching calculation in (15.270 and the SCETn matrix element in (I5.30p to be MS- 
subtracted. Similar to what we have seen in the last section, (I5.27P is free of large 
logarithms choosing /i ~ (Mm)^/^. We then evolve the scale down to /i ~ m using 
RG-techniques which implicitly resum logarithms of the type In fi^c/ yU. However, 
in contrast to (15.250 the SCETn matrix element (I5.30p is found to contain a large 
logarithm ln(27) at the low scale /i ~ m (the term vanishes in our example for the 
specific choice ^0 = m but is present in general for any variation of the scale with 
fi ^ m). In other words, we do not resum all large logarithms with the help of 
standard RG-techniques in SCETn- 

The remnant logarithms are related to the appearance of endpoint singularities 
as can be seen from the artificial decomposition into individual colhnear and soft 
regions in (I5.28P and (15.290 . We see that they do not belong to the variation of 
the factorization scale n but rather to the soft-coUinear cross-talk at fixed virtuality 
n"^ ~ m^. The low-energy SCETn matrix element in (I5.3O0 still depends on the 
high scale due to the fact that the meson rest frames are related by a large boost 
7 = 0{M/m). We conclude that these logarithms cannot be factorized into short- 
distance coefficient functions and we therefore call them non-factorizable. 

We emphasize in this context that we do not find contributions from soft-collinear 
messenger modes in our set-up. By looking at the scaling of soft-collinear momentum 
regions in the loop integrals we find that these modes always give power-suppressed 
contributions and that they are in particular not needed to describe the physics of 
endpoint-singularities. The reason for this is that the light quark mass m provides 
a physical IR cut-off in our calculation such that messenger modes with virtualities 
smaller than cannot contribute to on-shell amplitudes. 
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5.3.2 Factorization of tree level result 

For the following discussion it will be convenient to introduce a short-hand notation 
for the factorization formula (I5.17P 

Fi c:^ H, ^ + (j)B Ti ® 0^, (5.31) 

where the symbol ® represents the convolution integrals. We write the perturbative 
expansion of the form factors in the form 

F, = (^) Ff (5.32) 

fc=l 

and similarly for all other quantities in (15.311) . Notice that the perturbative expan- 
sion of the form factor starts at 0{as) in our set-up according to our discussion in 
Section [5.2.11 Since the perturbative expansion of Hi starts at 0{1) and the one of 
Tj at 0{as), the form factors take in LO the form 

^(1) ^ ^(0) ^(1) ^ ^(0) ^ ^ ^(0) (5^33) 

with the hard-scattering kernels given by (Km as r/^^ = ^® J(i). The jet function 
J^^^ has first been computed by Beneke and Feldmann in [81] using a different 
definition of the soft-overlap contribution C,{q'^)- Our convention corresponds to the 
one of Beneke and Yang in [76] from which we read off 

(o)_ „m_„m_cf\T) _ d'\r) _CP{r) _^ 



r>{T;u,uj) = — , (5.34) 

JMc 'jm uju 

where /m (/m) is the decay constant of the heavy (light) meson defined in HQET 
(QCD). In LO of the perturbative expansion we simply have /m = ff^, fm = fm^- 

In order to extract the overlap-contribution ^^^^ from (15.331) we have to calculate the 
distribution amplitudes of the Be and rjc meson to LO in the perturbative expansion. 
As we treat the quarks in the static approximation at leading power in the NR 
expansion, we simply find 

(j)f{uj) = 6{u - m), (m) = 6{u - 1/2). (5.35) 

The second term in (15.331) can now be calculated giving 

/NR/NR o^2r< \ s — 2 i = + 

<t^f ® ® 0f = -ImJu^ ^JL±f )s a ^ = - (5.36) 

Ac 7m I .^^ 
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We may now isolate ^^^^ using our explicit LO results from (15.71) together with fl5.33p 
and (15 .34^ . We indeed find a universal contribution in this way which reads 

We have shown in some detail how to extract the soft-overlap contribution .^(g^) from 
our explicit perturbative calculation in Section I5.2[ As the conceptually interesting 
aspects related to the physics of endpoint singularities and the appearance of non- 
factorizable logarithms only enter our calculation at NLO, we proceed with a similar 
analysis in the following section. 

5.3.3 Factorization in NLO 

At NLO of the perturbative expansion the extraction of the overlap-contribution is 
much more involved. We now start from 

pf) ^ ^(0) ^(2) ^ ^(1) ^(1) ^5^3g) 

with Tf^ = Cf ) ® J(2) +cf ) ® J(i). Whereas the hard coefficient functions have 
already been computed in [SI], the calculation of C^^^ and J*^^-* has been performed 
recently |lllf74f75f76] . In order to extract from (I5.38p . we still have to consider 
the corrections to the distribution amplitudes 0^'' and (j)^^ and to perform the 
subsequent subtractions. In the following we present our results for the distribution 
amplitudes whereas we relegate the extraction of the overlap-contribution to |112j . 



Light-cone distribution amplitude of rjc meson 

We start with the calculation of the leading twist light-cone distribution amplitude 
(LCDA) of the "light" r]^ meson defined by [mfST] 

{Vc{q)\ c{yh,l5c{x) |0)|(^_^),^, ^ -zfm q, f du S^^-y^^^'^) (5.39) 

J 

where we omitted an appropriate path-ordered exponential of gluon fields which 
makes the definition gauge- invariant. Notice that there is an additional leading 
twist two-gluon LCDA for fiavour singlet mesons which we will not consider here. 

The NR bound states are described by parton configurations with fixed momenta 
which correspond to the light-cone momentum fraction m = 1/2 for the r\c meson, 
cf. (I5.35p . Relativistic effects from collinear gluon exchange lead to modifications: 
First, there is a correction from matching SCETn to the NR theory. Second, there 
is the usual evolution under the change of the renormalization scale [36ll37] . In 
particular, the support region for the parton momenta is extended to < m < 1. 
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Figure 5.9: Relativistic corrections to the light- cone distribution amplitudes. 
The dashes line indicates the Wilson-line. 



The first-order relativistic corrections arise from the diagrams in Figure 15.91 Apart 
from these 1-loop diagrams, we have to include the wave-function renormalization 
of the external quark lines. We first consider the local limit of the light-cone ma- 
trix element (15.391) which determines the relativistic corrections to the NR decay 
constant. In this case, the diagrams with the gluon attached to the Wilson-line are 
absent and we find 



fm /, 



NR 

m 



1-6 



s'^F 

47r 



(5.40) 



which coincides with the result in |114] for equal quark masses. 

The NLO corrections to the leading twist distribution amplitude are found to be 



0«(«;/^) = 4C^lnAr 



1 + 



1/2 -u 



u9{l/2 -u) + {u^u) 



-4Cf 



1/2 -u 



ln(l - 2uf 



1/2 



[1/2 -u)^ 



-\- (u ^ u] 



ud{l/2 - u) 



(5.41) 



where we have introduced plus-distributions defined by 

\u [. . .]+ f{u) ^ j\u [. . .] (^f{u) - /(1/2) - /'(1/2) (u - 1/2) ) . (5.42) 



Notice that our result for the distribution amplitude obeys the evolution equation 



d 



0r,(M;/i) 



dv V{u, v) (pniy; /i) + 0{a 



din jj ' TT jQ 

where V{u,v) denotes the familiar Brodsky-Lepage kernel [36ll37] 



') 



(5.43) 



V — U / V 



e{v-u) + (i + -^ 



V — U / V 



u 



[u — v) 



(5.44) 



While writing this thesis an independent calculation of the leading-twist LCDAs of 
the ?7c and J/tp meson appeared in |115j where Ma and Si perform a similar analysis 
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factorizing perturbative effects from NRQCD matrix elements at leading power in 
the NR expansion. Our result is not in agreement with these findings which can 
already be seen by the fact that the distribution amplitude in |115j is not properly 
normalized to unity. We will have a closer look at the origin of this discrepancy 
in [TT^ . 



The factorization formula (15.381) contains the distribution amplitude (f)^^ in form of 
a convolution with the jet function J^^^ from (15.341) . We thus require the knowledge 
of the 1/u moment which is found to be 



du 



u 



2(3 



21n2) In A7 + 41n2- -7r2 + 12 

3 



(5.45) 



Light-cone distribution amplitude of meson 

The calculation of the LCDA for the meson goes along the same lines as for the 
r]c meson. However, important differences arise because the heavy fe-quark is to be 
treated in HQET which modifies the divergence structure of the loop integrals. As a 
consequence, the evolution equations for the LCDA of heavy mesons |116fll7j differ 
from that of light mesons 



We define the two-particle LCDAs of the Be meson following |118|8T] 



{Q\cp{z)h^mB,{p))\ 



lfM{fi)M 



al3 

(5.46) 



with p ~ Mw, t = w ■ z and Jm the decay constant in HQET (we again omitted 
an appropriate Wilson-line in the definition). The Fourier-transformed expressions, 
which usually appear in factorization formulas, are given through 



du e"*^*-^^ 



(5.47) 



In the following we focus on the distribution amplitude 0b (cj) = <I)^{uj) which enters 
the factorization formula (15.171) . In the strict NR limit it simply corresponds to a 
delta- function, cf. (15.351) . Similar to what we have seen in the case of the rjc meson, we 
now have to take into account relativistic effects from soft gluon exchange according 
to analogous diagrams to those in Figure 15.91 which now extend the support region 
of the distribution amplitude to < < oo. 



In the local limit we derive the corrections to the decay constant. We find 



/m(/x) 



J M 



1 + 



31n^ -4) +0(a 
m 



(5.4J 



Notice that the decay constant of a heavy meson in HQET is scale- dependent. 
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The NLO corrections to the distribution amplitude read 



'^g)(u;;/i) = 2C^ln 



u — m) 



Cf S{lj - m) ( ^ In^ 



uj6{m — uj) — m9{uj — m)6{A — uj) 
m(m — uj) 

(A - m)2 m? 12 



+ 



(^-A) 



w — m 



+ 2C, 



{m + u) [uj6{m — a;) + m6{ijj — m)9{A — u)] 



—6{lu — m) I In 



m 



+ 



m{m — u)"^ 
2m 



+ 



{m + u)e{u- A) 



(w — m) 



, , — mo [00 — m) \ 2 In 1 1 — 

A — m A — m/ \ m m 2 



(5.49) 



with an analogous definition of plus-distributions as in fl5.42p . One remark is in order 
concerning the dependence of our results on an ad-hoc parameter A > m. In the 
computation of the LCDA we encounter IR singularities for uj = m (similar to the rjc 
case for M = 1/2) which we regularize with the help of adequate plus-distributions. 
We rewrite for instance 



{u — m) 



[w 



m 



|l+2e 



9{u 



m] 



w 



m 



(5.50) 



This is to be understood in the sense of a distribution which is convoluted with a 
kernel f{uj) that is assumed to be non- vanishing for uj = m and to scale as l/cu for 
large uj (as usual in application of QCD Factorization). The point to notice is that 
the (IR regular) plus-distribution gives rise to UV singularities in the convolutions 
as f{ijj) ~ l/oj is replaced by f{uj) — f{m) ~ 1 for large u. In order to avoid 
these artificial UV divergences we introduce an ad-hoc cutoff A > m in the plus- 
distributions such that f l5.50p is replaced by 

-2e 



9{uj — m) 
(w — m)^+2^ 



{u - m)9{A - u)' 



w — m 



+ 



e{uj - A) (A 



m] 



w — m 



2e, 



■6{u -m) +0{e). 

(5.51) 



This expression is thus free of UV singularities whereas the IR divergence has been 
isolated in the local term. Notice that the dependence on the cutoff parameter drops 
out in the convolution with an arbitrary kernel f{uj). This can also be seen below 
in our expression for the 1/uj moment of the LCDA in (15.561) . 



The distribution amplitude in (15.491) obeys the evolution equation 



0i?(^^;/^) 



TT 



du' 7+(u;, u'; fi) 4>b{uj'; fi) + C(aJ, 



where the anomalous dimension 7+ (07,0;'; /i) is given by [11611117] 
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u'{u' — u) u{u — u') 



(5.52) 



(5.53) 
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This can be verified most easily by integrating f l5.52p over a test function f{uj) sucli 
tliat tlie dependence of the LCDA on the cutoff parameter A drops out. 

In contrast to the rjc case, the normalization of the Be distribution amplitude is 
ill-defined. Imposing a hard cutoff Ayv ^ rri and expanding to first order in m/Auv, 
we derive 



Auv 



1 



An 



2 A2 



+ In 



+ 12 



A2 



+ 0{al 



(5.54) 



and similarly for the first moment 



Auv 

d(jj UJ (psi^', 



An 



2 In 



6 



(5.55) 



The last two expressions provide model-independent properties of the distribution 
amplitude which have been studied within the operator product expansion in [119] . 
Our results are in agreement with these general findings. Notice, however, that our 
subleading terms in m/Ayv cannot be compared with |119] as we deal with massive 
light quarks here. 

We finally quote our result for the l/oo moment which appears in the factorization 
formula f l5.38p in form of the convolution with the jet function J*^^^ from f l5.34p 



UJ 



m 



2 m2 



m'^ A 



(5.56) 



Notice the presence of a double Sudakov logarithm as in fl5.25p . 



Part IV 
Conclusion 



Conclusion 



The dedicated study of charmless B decays has become one of the most active 
and promising fields in particle physics. It aims at a precision determination of 
the flavour parameters in the Standard Model and may help to reveal the nature 
of CP violation. On the experimental side, the B factories BaBar and Belle are 
continuously accumulating larger data samples and provide us with measurements 
of unprecedented precision. On the theory side, the main challenge consists in the 
quantitative control of the complicated hadronic dynamics. 

QCD Factorization represents a model-independent framework to compute hadronic 
matrix elements from first principles. Based on a power expansion in Aqqd /rrib, it 
allows for a clear separation of short- and long-distance effects in certain classes of 
B decays. Whereas the development of Soft-Collinear Effective Theory provided 
the means to put factorization formulas onto a firmer basis, the implementation of 
higher-order QCD corrections has started only recently. 

In this thesis we have mainly addressed the most challenging class of charmless B 
decays: exclusive nonleptonic B decays with its most prominent example B ^ mi. 
This wide class of decays is of primary importance as most of the observables at 
the B factories are related to it. Our calculation considers one class of the NNLO 
corrections to hadronic two-body decays within QCD Factorization. Whereas the 
(1-loop) spectator scattering contributions have been considered recently by various 
groups , we have presented the first calculation of the 2-loop vertex corrections 

which constituted the missing piece for a full NNLO analysis of the topological tree 
amplitudes in QCD Factorization. 

The knowledge of the NNLO corrections is particularly important with respect to 
strong interaction phases and hence direct CP asymmetries. As strong phases arise 
first at 0{as) in QCD Factorization, they were known so far to LO accuracy only. 
As in any perturbative calculation the considered corrections are thus necessary to 
eliminate large scale ambiguities of the LO result. They may even drastically change 
the pattern of direct CP asymmetries in QCD Factorization. 

We indeed found that the considered corrections can exceed the formally leading 
contributions whenever the latter are accompanied by small numerical coefficients. 
This is the case for the imaginary part of the colour-allowed tree amplitude which 
receives its dominant contribution from the considered 2-loop vertex corrections. 
The respective contribution to the colour-suppressed tree amplitude also turned out 
to be sizeable. However, in absolute terms our corrections are small. In other words. 
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we did not encounter any source of sizeable strong phases from the perturbative 
calculation. We further remark that we have observed an accidental cancellation 
between the individual NNLO contributions which resulted in a moderate correction 
to the NLO result for the imaginary part of the tree amplitudes. 
The calculation of the real part of the topological tree amplitudes is more involved. 
We have presented a preliminary result for the colour- allowed tree amplitude which 
is based on certain technical simplifications. More importantly, we have already 
tackled the most challenging part of the calculation which consisted in the evaluation 
of a large set of 2-loop integrals. We could further verify that factorization holds 
at the considered order in perturbation theory which is a non-trivial statement as 
the formal all-order proof is still missing. We emphasize that this is the result of 
a highly complicated subtraction procedure which can also be considered as a very 
powerful cross-check of our calculation. 

From the technical point of view we have presented a calculation of 2-loop hadronic 
matrix elements. Even in the wider sense of perturbative corrections to weak decays, 
there exist only very few calculations of this type. In the context of inclusive decays, 
2-loop matrix elements have been considered for B — »• Xsi~^i~ |120lll21j and 3-loop 
matrix elements for B X^'y \1'2'2 \. Concerning the exclusive modes 2-loop matrix 
elements have been calculated for B K*'~f or i? — > pi^i~ |123] . However, our 
analysis represents the first calculation of 2-loop hadronic matrix elements for the 
more complicated nonleptonic decays which are mediated by a 6 — uud transition. 
For this class of decays 2-loop hadronic matrix elements have never been calculated 
so far, neither in the inclusive nor in the exclusive channel. 

In our final analysis we examined the formal factorization properties of heavy-to- 
light form factors. The development of Soft-CoUinear Effective Theory has led to 
a deeper understanding of the QCD dynamics in heavy-to-light transitions and the 
i? — > TT form factor (at large recoil) has emerged as the central object in these stud- 
ies. Rather than stating a complete separation of short- and long-distance effects, 
the factorization formula contains a classification into (non-calculable) symmetry 
conserving contributions, the so-called soft-overlap contribution, and perturbative 
symmetry breaking effects. Remnant short-distance effects in the soft-overlap cannot 
be factorized by standard means as this would lead to endpoint-divergent convolu- 
tion integrals. The question arises if one can understand this non-factorization of 
soft and coUinear effects from the viewpoint of the effective theory. 
In order to address this question we have considered a simplified scenario. We 
investigated heavy-to-light form factors between non-relativistic bound states which 
can be addressed in perturbation theory. We have performed a NLO analysis of 
these form factors and looked for an interpretation in terms of the factorization 
formula. We showed that the soft-overlap contribution can be calculated in our 
set-up. We have addressed the question of the origin of formally large (Sudakov) 
logarithms and argued that one class of logarithms related to endpoint-singularities 
cannot be resummed with standard RG techniques. We in particular did not need to 
introduce soft-coUinear messenger modes to describe the physics of the soft-overlap 
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contribution. As a byproduct of our analysis, we have calculated leading-twist light- 
cone distribution amplitudes for non-relativistic bound states which can be applied 
for the description of Be and rjc mesons. 

To summarize, the development of QCD Factorization and Soft-Collinear Effective 
Theory has substantially improved our understanding of QCD in B meson decays. 
To date, the QCD Factorization predictions are in good overall agreement with the 
experimental measurements. In order to perform a precision study of B decays, we 
have to reduce theoretical and experimental uncertainties as much as possible. This 
thesis has contributed to the former by calculating higher-order QCD corrections 
to exclusive charmless B decays. Apart from the perturbative corrections, further 
improvements on the theory side concern the control of power corrections and precise 
determinations of the hadronic input parameters. We think that it will be exiting 
to confront these improved theoretical predictions with updated experimental data 
in the near future. 
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Appendix A 
Master Integrals 



In the following we present the analytical results for all Mis that appeared in our 
calculations from Chapter [HI HI and O The Mis are normalized according to 



We give the results up to the order in e that was required in our calculations. For 
some simpler Mis, we obtain the results in a closed form in e. 

A.l Hadronic two-body decays I 

In the calculation from Chapter [S] we found 14 Mis which are shown in Figure I3.3[ 
The corresponding particle propagators Vi can be found in (13. 6p . Although the 
calculation in Chapter [3] requires the knowledge of the imaginary part only, we give 
the full results of the Mis here in view of our extension in Chapter |H 

1— loop integrals 

Some 2-loop integrals can be written as products of the following 1-loop integrals 



[dk] 



V{l-e) 



(A.l) 






(A.2) 
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3— topologies 
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The following Mis can be written in terms of the variable ^ = /u and 
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^ In^ ^ + ^ In^ r/) - (1 - 2r/) (^2Li2(r7) + 21n^lnr7 - In^r/) +r/ln^^ 



TT 



+ 51n^+ (4^1n^-8^1nr/ + 4r/)y + i7r 



- 2(1 - 2r]) Inr/ 



5-^nn^-21nr7j +2r/ln^ 



(A. 18) 



with 



^(47) 



.(47) 



.(47) 



2_2N-l-2e 



V1VIV22V2Z 



J=~2 



(A.19) 



-e(lne + ivr 

-^4- (1 -2r7)(^2Li2(r/) + 21n{lnf/-ln2r7j +r7ln^^ + 21n^ 



+ 2in(l + fj\n^ - (1 - 2r7)lnr7^ 



27r2 



(A.20) 
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5— topologies 



[dk][dl] 



1 



j=-i 



(A.21) 



with 



Inu — Inu] \- in 

6 



TC 



U2{u) 



In M — In M In M 



.(51) 



4l7r^ 



+ 13Li4(u) — 682,2 (w) — Si,3(u) — ^lllnw + 6 Inu^Li^lu) + ^ In"^ m 

f 1 1 In u — In Si 2 (u) In'^ u \nu -\ — In^ u u In m In'^ u 

V / 6 2 6 

- ^5 In^ u+ 12 Inw In-u — In^ M — 5Li2('u) jLi2('u) — Cs (^H Inw — Tln-uj 

2C3 + 7Li3(n)-2Si,2H 
— ^ In ti + 2 In u j Li2 (u) — - In^ m + - In^ u In m — In m In^ u 

(A.22) 



\26U2{u) -41n^M + leinwlnu + ln^ uj— + in 



5 In M — In M I — 



witli 



„(52) 



.(52) 



[dk][dl] 



1 



= (m 



TT 



C3 - Si,2(^i) + «vr ( —- U2{u) 



^-2S2,2(u)-5Si,3(«) + — Li2(n) 



i-^i 7(3 - 2Li3(M) - 5Si,2(m) 



j=-i 



(A.23) 



(A.24) 
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6— topologies 

The non-planar massless 6-denominator integral has been considered in |1 24111 26] . 
Recently it has been formulated in terms of hypergeometric functions which can 
be expanded to arbitrary order in e with the help of the Mathematica package 
HypExp [127]. Our result 





(A.25) 



is in agreement with the previous findings. 
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A. 2 Hadronic two-body decays II 

Our analysis in Chapter H] involves the calculation of the Mis from Figure 14.11 The 
corresponding particle propagators Vi can be found in fl3.6p . So far we have not 
yet calculated the Mis with an internal charm quark (wavy lines) which will be 
considered in |106] (these amount to four Mis out of 22). 



3— topologies 







] 


with 








1 

~ 2' 






3 

~ 2 ' 


u 

^4' 


"o 


= 34 


' y 




15 

~ T 


+ 7r 



(31) 



[dk][dl] 



(m 



2\l-2e 



J2dt'>e' + 0{e')\ (A.26) 



J=-2 



13u ^ , . l-u\ 
8 ^ ^ 2m 



115 vr^ 

— + y )m + Li3(M) + 4Si,2(m) + 



1 + Qu- 2u' 
2u 



Li2(M) 



+ 



In M — 



u 



13(1 -M^ 

4m 



21 



57r^ 
36" 



3C3- 



In u, 

865 137r2 C3 



32 



/ N 1 + 6m - 2^2^ , , 
+ 16Si,3(m) + Li3(m) 



2u 



12 2 
4 - 12m - 2m2 

M 



U + Li4(M) - 8S2,2(^i) 



Si,2(M)+3Li2(M) 



u 



1-T/N 4n_ 13, 2- vr , , _ 

3 In M Li2(m) + - In w —In u + [ -— + — ) In m 

3 2 V 8 3 



1 

+ — 

4m 



13 + 24 + 



n - 26n' 



Li2(u). 



(A.27) 



= /[^^H^^]^^^ = (m^)-^^ ||^df)e^ + 0(e3)| (A.28) 

with 
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^(32) 


9 

~ ~2 


3 




= C3- 


27 

T ~ 



TT 



-^Li2(n) + -f21n2n-31nn), 

Si,2(w) + Li3(M) + 3Li2(M: 



6 Inw Li2(M) + - In^ M - 6 In^ M + (9 + — ] In 

3 V 3 / 



,(32) OA 81 2 Stt^ qt-2/ X 12(2 -m) 

4 = SCa - 17 - 37r - — - 3Li2(m) Si,2(m) 

2 3o 



+ 



u 

+ - 



1 - 2u 



u 



U^{u) - 8S2,2(m) + 16Si,3(u) + 3Li3(u) + (9 + — ]U2{u) 



24 In M Si,2 iu) - 6 In M Lis («) + 12 In^ u U2 {u) + - In^ m - 18 In m Li2 (m) 



4 ^ 

8 In^ M + (|l8 + ^) In^ u - {21 + 2'n'^ - 2(3) In 



u 



(A.29) 



The massive sunrise integral can be found in |128] . Our result for the equal mass 
case is in agreement with these findings and reads 




4— topologies 



(A.31) 



with 
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Uq 



(41) 



5 M , _ 
- + - lllM, 

2 u 

19 Tl"^ _u 

Y ^ y ^ M 

65 Stt^ 



Li2 (m) + - In^ m — 5 In m 



(41) 



2+ 3 

^ 15 , 2 - 

+ — m u 



211 197r2 Svr^ n 
^ ^ ~3~ ^ "36" ~ ~ u 



Li3(M) — 3Si^2('") — 31nuLi2(M) — - In^w + 5Li2(u) 



19 + — Inu 



Li4(n) — 3S2,2(m) + 9Si,3(u) — 3 \\iu\Aj,{u) 



9 9 15 

+ 91nMSi,2(u) + -ln^uLi2(M) + -In^'u + 5Li3(u) - 15Si,2(u) In^w 

'2 8 '2 



15 Inw Li2('u) 



19 



27r2 



lOvr^ 



(Li2(M) + -ln2M) - (^65 + ^-2C3j Inw 

(A.32) 



/ 



with 



7(42) 



7(42) 



1 

2' 

5 M 1 _ 

- + - Inu, 
2 u 

19 vr^ M 

Y ^ y ~ M 



1 



[m 



2\-2e 



J=-2 



(A.33) 



(42) 



65 57r ^ ^ M 

^ h 4(3 

2 2 ^^ u 



2Li2(M) + 21n^M - 51nM , 

4Li3(M)-8Si,2(M) + 10Li2(M) 



\t? u 



- 81nMLi2(n) + lOln^u - (19 + 7r2)lnu 



4"' 



211 197r2 TT^ ^ M 

— + + — + 20C3 - - 

2 2 2 ^'^ u 



8Li4(M) - 16S2,2(m) + 32Si,3(m) + - In^ u 



161nMLi3(M) + 321nMSi,2(M) + 16 In^ m Li2(M) + 20Li3(M) -401nMLi2(M 
40 

40Si,2(m) -—\Yi^u + 2(19 + 7r2)(Li2(M) + In^ u) - (65 + 571^ + 8(3) In u 

(A. 34) 
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[dk][dl] 
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2,-2e J ^^(43) ^.^(^^^2^ 
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(A.35) 
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"o 



d\ 
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2 



2C3 + 2Li3(M) +MLi2(M) 



+ 4S2,2(m) + 3Li2(M) + 41nMLi3(M) + (5 - 3M)(Li3(M) - (3) 
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(A.36) 
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d^ 
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d'^f 


1 
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d^lf 


1 + u 
2u 


] 

+ - 


"0 


= -2 - 


6 


df'^ 


= 2C3 + 


2(1 



2^2 



2\-l-2e 



(A.37) 
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-Li2(M) 
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6m 



2m2 
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2(3 + m^)Si,2(m) + 2(1 + 2u^) lnMLi2( 



u 



l + hu\. . , 4(1 + ^2) 
Li3(m) H 



, 3_ , (1 + m)(1 + 3u)^. , , 
In"* u + ^ '-j- Li2(m) 



2 

- {I - Au - u^)\-n^ u + {2{l + u^) + (3 + «^)y) Inw 
27r2 TT^ 2(2 - u 
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+ 



-C3 + -T 



1 -3m^ 



Li4(u) -2(3-5n2)S2,2(M) 



+ 2(1 - 9m2)Si,3(m) - 2(1 - 2u^) lnMLi3(M) - 12(1 + u^) Inw Si,2(m) 
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3u^U2{uy - (5 + 7u') \n' u U2{u 



4(1 + ^2) (1 + 5^ 
In u Lia^Mj 



+ 2(3 - 6m - m^)Si,2(m) + 2(1 - 5m - 2u^) \nuU2{u) + 



4(1 -Au-u^ 



2 2 
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2 

+ (^2(1 - 4m - M^) + (3 - 10m - M=^)y + 2(2 - ^2)^3^ 



u 



(A.38) 
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T + ¥ 
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65 57r" 
2^6 



(A.39) 
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2u 
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5 + y )Li2(M) 
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(A.40) 



Jmidi]^^ 
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(A.42) 



136 



A. MASTER INTEGRALS 



The following Mis give rise to HPLs with parameter —1 
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(A.43) 
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. 27r2 ^ ^ TT^ 2(2 + u)(1-2m) , , 
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(A.44) 
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[dk][dl] 



1 



[m 



2,-2. J ^^(49) ^.^^(^2) 
J=-2 



with 



j(49) 
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- 7C3 In 2 - 8Li4(|) - 12Li4(M - 1) - (^4 In^ m - m In m + j Li2(M - 1) 

+ 101nMLi3(u - 1) - (^In^M - -In^M + — Inw - 7(3 —j ln(2 - u) 

X./ -,^ ^./N 4-5m + 2^2 6 - 12m + 

- uLi3(M - 1) - 2^(m) + Li3(u) Si,2H 



6 - 10m + 4^2 

u 

15 - 16m 



ln-uLi2('u) 
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(A.48) 



where we defined an auxiliary function A{u) which is given by the integral 

r 1 

A(u)= du if(l,0,l;M')- (A.49) 

Jo 2 - m' 
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[dk][dl] 
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The auxiliary function A{u) can be found in flA.49p . 
5— topologies 
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(A.61) 



where we have introduced a constant Cq, which arises in the calculation of the 
boundary condition for the 6-topology MI. We could not find an analytical expression 
for C(fl. With the help of our implementation of the method of sector decomposition, 
we obtain Cq = —60.2 ± 0.1. 



"'^We thank R. Bonciani and A. Ferroglia for pointing out an error in an earlier version of this 
thesis. 
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A. 3 Heavy-to-light form factors 

The Mis in our calculation from Chapter E] are shown in Figure The correspond- 
ing propagators Vi can be found in (15.101) which involve two external momenta w 
and w' with w'^ = w' 

of the Mis in m/M and write s = M/{A'ym) = 0{1). 



1 and 'f = w-w' = 0{M/m). We extract the leading power 
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4— topologies 
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[dk] 
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Appendix B 

Form factors in NLO 



We summarize our NLO results for the transition form factors F^{q^), F-{q^) and 
Friff) which can be defined in analogy to (11 .271) for the i?c i]c transition. All 
results are expressed in terms of pole masses M = rrih and m = and correspond 
to the leading power in the HQE in m/M. The relation between the momentum 
tramsfer and the relativistic boost 7 = 0{M/m) can be found in (15.51) . We 
further introduce s = M/{4'ym) = 0{1) and refer to (15. 7p for the LO expressions. 
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l + 2s 3(1 + 2s) 9(1 + 2s) 

where 5+ refers to the contribution from the flavour singlet diagrams which reads 
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with tlie flavour singiet contribution given by 
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(B.5) 



and the flavour singlet contribution is found to be the same as in (IB. 40 . i.e St = S^. 
Notice that the tensor form factor Fx^q^) has a residual scale dependence of 0{a1) 
in NLO due to the fact that the tensor current is not conserved. 
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